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SECTION  1 


INTRODUCTION 

The  conventional  propulsive  nozzle  is  an  axisynwetrlc  converging- 
diverging  nozzle  terminating  in  a  plane  perpendicular  to  the  axis  of 
the  nozzle.  If  the  flow  entering  the  nozzle  is  axi symmetric,  then  the 
flow  within  the  nozzle  is  also  axi symmetric*  and  the  resultant  thrust 
vector  lies  along  the  nozzle  axis. 

In  some  applications,  such  as  a  nozzle  exiting  through  the  side 
of  a  missile  where  the  nozzle  axis  is  not  aligned  with  the  missile 
axis,  the  nozzle  terminates  along  the  line  of  intersection  of  the 
nozzle  contour  and  the  missile  outer  skin.  Such  a  nozzle  is  called  a 
scarfed  nozzle.  Figure  1  illustrates  several  possible  scarfed  nozzle 
configurations. 

Unbalanced  side  forces  are  generated  in  scarfed  nozzles,  and  the 
resultant  thrust  vector  does  not  lie  along  the  nozzle  axis  or  the 
missile  axis.  When  two  or  more  identical  scarfed  nozzles  are  located 
symmetrically  around  the  missile  axis*  the  resultant  thrust  vector 
does  lie  along  the  missile  axis.  However,  even  in  that  case  the 
missile  axial  thrust  is  reduced  below  the  nozzle  axial  thrust  due  to 
the  misalignment  of  the  nozzle  axis  and  the  missile  axis. 

Expansion  waves  nr  shock  waves,  depending  on  the  pressure  ratio, 
emanate  from  the  nozzle  exit  lip  contour  as  the  internal  flow  leaves 
the  protection  of  the  solid  wall  and  flows  into  the  surrounding 


(c)  Conical  nozzle,  contoured 


(d)  Contoured  nozzle,  contoured 


extension.  extension. 


(g*  Conical  nozzle,  cylindrical  (h)  Contoured  nozzle,  cylindrical 
extension.  extension. 


Figure  1.  Several  possible  scarfed  nozzle  configurations. 


atmosphere.  When  scarfing  Is  not  too  severe,  the  waves  emanating  from 
the  top  of  the  no2zle  at  the  initial  point  of  scarfing  do  not  intersect 
the  opposite  wall,  and  the  portion  of  the  fiowfield  affecting  the 
pressure  on  the  nozzle  wall  remains  axi symmetric.  That  is  the  case 
considered  In  the  present  analysis.  When  scarfing  is  severe,  the 
aforementioned  waves  Intersect  the  opposite  wall,  and  the  flow  down¬ 
stream  of  that  intersection  is  fully  three-dimensional.  That  case  is 
not  considered  In  the  present  analysis. 

The  present  report  presents  an  analysis,  and  a  computer  program 
for  implementing  that  analysis,  for  determining  the  performance  of 
scarfed  axi  symmetric  propulsive  nozzles  when  the  scarfing  is  small 
enough  so  that  the  portion  of  the  flowfield  affecting  the  pressure  on 
the  nozzle  wall  remains  axi  symmetric. 
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SECTION  II 


NOZZLE  GEOMETRIC  MODEL 

1.  INTRODUCTION 

Several  scarfed  nozzle  configurations  are  illustrated  In  Figure  1. 
Six  specific  nozzle  configurations  are  considered  in  the  present  study. 
These  configurations  are  illustrated  in  Figure  2. 

Figure  2(a)  illustrates  a  scarfed  conical  nozzle,  where  the 
nozzle  contour  is  specified  analytically  as  a  straight  line. 

y  *  a  +  bx  (1) 

Figure  2(b)  illustrates  a  conical  nozzle,  specified  by  equation  (1), 
followed  by  a  scarfed  conical  extension  starting  at  the  end  of  the 
conical  nozzle,  where  the  nozzle  extension  is  specified  analytically 
as  a  straight  line. 

y  -  e  +  f x  (2) 

When  the  slope  f  Is  zero,  a  cylindrical  extension  is  obtained. 

Figure  2(c)  illustrates  a  scarfed  quadratic  nozzle,  where  the 
nozzle  contour  Is  specified  analytically  as  a  second-order  (l.e., 
quadratic)  polynomial. 


a 


(a)  Conical  nozzle. 


(b)  Conical  nozzle  followed  by  a 
conical  (or  cylindrical)  extension 


Figure  2.  Nozzle  geometric  models  considered. 


y  *  a  +  bx  +  o«2 


(3) 


figure  2(d)  Illustrates  a  quadratic  nozzle,  specified  by  equation  (3), 
followed  by  a  scarfed  conical  extension  starting  at  the  end  of  the 
quadratic  nozzle,  where  the  nozzle  extension  Is  specified  by  equation 
(2). 

‘  4 

Figure  2(e)  Illustrates  a  scarfed  tabular  nozzle,  where  the  basic 
nozzle  contour  is  specified  In  tabular  form. 

y*  «  f(xf),  (1  *  1,...,  n)  tabular  (4) 

Figure  2(f)  illustrates  a  tabular  nozzle,  specified  by  equation  (4), 
followed  by  a  scarfed  conical  extension  starting  at  the  end  of  the 
tabular  nozzle,  where  the  nozzle  extension  Is  specified  by  equation  (2). 

Figure  3  presents  a  meridional  plane  view  o?  the  geometric  model 
employee  for  the  scarfed  nozzk  flowfleld  analysis  considered  in  the 
present  Investigation.  The  nozzle  consists  of  a  conventional  exisym- 
metrle  throat  and  supersonic  expansion  contour  up  to  point  E  where 
scarfing  begins.  That  portion  of  the  scarfed  nozzle  is  called  the 
basic  nozzle.  The  scarfed  section  fror  point  E  to  point  F  Is  a  conical 
extension  to  the  basic  nozzle.  That  portion  of  the  scarfed  nozzle  is 
called  the  nozzle  extension.  In  the  present  analysis,  the  scarfed 
section  always  starts  at  point  E. 

The  specification  of  the  scarfed  nozzle  geometry  Is  discussed  in 
In  the  following  paragraphs. 
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2.  BASIC  NOZZLE  GEOMETRY 

The  basic  nozzle  geometry  consists  of  a  double  circular  arc  throat 
contour  joined  smoothly  to  a  supersonic  expansion  contour.  The  super¬ 
sonic  expansion  contour  can  be  conical*  quadratic*  or  tabular.  The 
specification  of  the  basic  nozzle  geometry  Is  discussed  in  the  following 
paragraphs. 

Throat  geometry.  The  throat  of  the  i-aslc  nozzle  consists  of  a 
double  circular  arc  contour  that  attaches  smoothly  to  the  supersonic 
expansion  contour  at  point  A.  The  throat  geometry  Is  completely 
specified  by  the  throat  radius  yt»  the  upstream  throat  radius  of 
curvature  ptd*  the  downstream  throat  radius  of  curvature  p^,  and  the 
throat  attachment  angle  e.  where  the  supersonic  expansion  contour 

O 

attaches  smoothly  to  the  circular  arc  throat  contour.  The  location  of 
point  A  Is  determined  from  the  following  two  equations. 

‘  ptd  s1"  es  (5) 

V  V  ptd (1  - 005  V  («> 

The  upstream  throat  radius  of  curvature  is  used  in  the  calculation 
of  the  supersonic  Initial-value  line  which  spans  the  nozzle  throat  from 
point  T  to  the  point  on  the  nozzle  centerline  where  the  Mach  number  Is 
1.0. 

Conical  nozzle.  Wher  the  supersonic  expansion  contour  is  conical  * 
the  cone  angle  a  must  equal  the  throat  attachment  angle  6A.  Five 

o 


options  exist  for  specifying  the  conical  nozzle  contour.  /  / 

■'■'I.  Specified  throat  attachment  angle  6a  and  nozzle  length  xe.  In  this 
case,  the  nozzle  exit  lip  radius  ye  Is  determined  from 

ye  *  ya  +  (Xe  -  xa)  tan  e&  (7) 

2.  Specified  throat  attachment  angle  9a  and  exit  lip.  radius  ye.  In 
this  case,  the  nozzle  length  xe  Is  determined  from 

xe  *  xa  +  (ye  -  ya)  /  tan  8a  (8) 

3.  Specified  threat  attachment  angle  Oa  and  nozzle  area  ratio  e  * 
(ye/yt)2-  In  this  case,  the  nozzle  exit  lip  radius  ye  is  determined 
from 

ye  *  (9) 

The  nozzle  length  xe  Is  then  determined  from  equation  (8). 

4.  Specified  nozzle  length  xe  and  exit  lip  radius  ye.  The  throat 
attachment  angle  ea  is  determined  from 


5.  Specified  ruszzle  length  xfi  and  nozzle  area  ratio  e.  The  nozzle  exit 
lip  radius  ye  is  determined  from  equation  (9)  and  the  throat  attachment 
angle  ea  is  determined  from  equation  (10). 


r^j'»  * » . ' 


>.v. 


v?  **. 


In  all  five  cases,  point  A  (xa,  ya),  point  E  (xe,  ye),  and  the 
throat  attachment  angle  0a  are  known.  The  equation  for  the  supersonic 
expansion  contour  is  given  by 

y(x)  3  a  +  bx  (11) 

where 


b  =  -  tan  a  and  a  =  ya  -  bxa  (12) 

Quadratic  nozzle.  When  the  supersonic  expansion  contour  is  quad¬ 
ratic,  the  supersonic  contour  must  attach  smoothly  to  the  circular  arc 
throat  contour,  and  the  nozzle  length  xe  must  be  specified.  The  equa¬ 
tion  for  the  supersonic  expansion  contour  is 

y(x)  =  a  +  bx  +  cx^  (13) 

Three  options  exist  for  specifying  the  supersonic  expansion  contour. 


1.  Specified  throat  attachment  angle  6a,  nozzle  exit  lip  angle  8e,  and 
nozzle  length  xe.  Substituting  the  three  known  values  (xa,ya),  (xa,ea), 
and  (xe,0e)  into  equation  (13)  and  its  derivative  (i.e.,  dy/dx  =  b  +  2cx 
3  tan  e)  yields 


tan  6  -  tan  6 
2(xe  -  xa) 


CHI 
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b  =  tan  6a  -  2cxa 


(15) 


a  =  ya  ■  bxa  -  “a2  05) 

2.  Specified  throat  attachment  angle  0a?  nozzle  length  xe,  and  exit  lip 
radius  ye.  Substituting  the  three  known  values  (xa,ya),  (xa,0a)»  and 
Ue»ye)  into  equation  (13)  yields 


tan  Ga  -  2cxa 

(IS) 

*  bxa  -  “a2 

(19) 

3.  Specified  throat  attachment  angle  0a,  nozzle  length  xfi,  and  nozzle 
area  ratio  e.  In  this  case,  the  nozzle  exit  lip  radius  ye  is  given  by 
equation  (9).  The  coefficients  a,  b,  and  c  are  given  by  equations  (17) 
to  (19). 

Tabular  nozzle.  If  the  supersonic  expansion  contour  is  tabular,  a 
set  of  (x,y)  pairs  spanning  the  region  from  point  A  to  point  E  must  be 
specified.  These  points  must  be  chosen  so  that  a  smooth  transition  to 
the  circular  arc  throat  contour  is  obtained.  The  first  point  in  the 
table  is  a  point  just  downstream  of  the  throat  attachment  point,  point  A. 
The  last  point  in  the  table  specifies  the  nozzle  exit  lip  point,  point  E. 
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3.  NOZZLE  EXTENSION  GEOMETRY 

The  nozzle  extension  consists  of  a  conical  section  beginning  at 
the  end  of  the  basic  nozzle,  point  E  In  Figure  3,  and  ending  at  point 
F.  The  scarfing  angle  6  Is  the  angle.  In  the  meridional  plane, 
between  the  nozzle  axis  and  the  missile  axis.  For 

a  conventional  propulsive  nozzle,  the  scarfing  angle  0  Is  zero.  In 
general,  the  Intersection  of  a  body  of  revolution  such  as  the  nozzle 
contour  with  a  second  body  of  revolution  such  as  the  missile  outer 
envelope  will  not  yield  a  curve  of  Intersection  that  lies  within  a 
plane.  However,  If  the  diameter  of  the  missile  Is  much  larger  than  the 
diameter  of  the  nozzle,  the  curve  of  Intersection  lies  In  a  surface 
that  approaches  a  plane  as  the  diameter  of  the  missile  Increases.  In 
the  present  analysis,  the  curve  of  Intersection  Is  assumed  to  lie  In  a 
plane. 

The  geometry  of  the  nozzle  extension  Is  completely  specified  by  the 
attachment  point  (xe,y0),  the  angle  of  the  conical  extension  6f,  and 
the  length  of  the  conical  extension  Xf . 

If  the  nozzle  extension  attaches  smoothly  to  the  basic  nozzle 
(l.e.,  8e  *  6f),  then  a  conttnous  flow  occurs  across  the  transition 
region.  If  6f  is  less  than  6e,  then  an  oblique  shock  wave  emanates  from 
point  S  and  propagates  downstream  In  the  nozzle  extension*  Both  possi¬ 
bilities  are  accounted  for  in  the  present  analysis.  If  0f  is  greater 
than  0e,  then  a  centered  expansion  wave  emanates  from  point  E  and  prop¬ 
agates  downstream  in  the  nozzle  extension.  That  possibility  is  not 
considered  in  the  present  analysis. 
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SECTION  III 
FLOWFIELD  MODEL 


1.  INTRODUCTION 

The  flowfield  model  consists  of  four  parts. 

(a)  The  transonic  flow  analysis  in  the  throat  region  of  the  basic 
nozzle. 

(b)  The  supersonic  flow  analysis  in  the  basic  nozzle. 

(c)  The  determination  of  the  shock  wave  (or  Mach  line  if  0e  =  0f) 
which  emanates  from  the  point  of  attachment  of  the  nozzle 
extension  to  the  basic  nozzle  (point  E  in  Figure  3)  and  propagates 
downstream  into  the  nozzle  extension. 

(d)  The  supersonic  flow  analysis  in  the  nozzle  extension  downstream  of 
the  oblique  shock  wave. 

The  specification  of  these  four  flowfield  models  is  presented  in  the 
following  paragraphs. 

The  flow  is  assumed  to  originate  upstream  of  the  nozzle  in  a 
uniform  flow  region  having  constant  stagnation  pressure  and  temperature 
and  Tt,  respectively.  The  flowing  fluid  is  assumed  to  be  a  therm¬ 
ally  and  calorically  perfect  gas  (i.e.t  constant  molecular  weight  and 
specific  heats).  The  presence  of  condensed  phases  and  chemical  reac¬ 
tions  are  neglected,  and  boundary  layer  effects  are  considered 
negligible.  Consequently,  the  flowfield  is  isentropic  everywhere 
except  within  the  oblique  shock  wave,  which  Is  described  by  the 
standard  oblique  shock  wave  relationships. 
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The  flowfield  in  the  transonic  region  and  in  the  basic  nozzle  is 
irrotational  (i.e. ,  constant  entropy  and  stagnation  enthalpy  throughout) 
since  the  flow  originates  in  a  uniform  flow  region  and  is  isentropic. 

The  flowfield  is  rotational  downstream  of  the  oblique  shock  wave  due  to 
the  entropy  gradient  produced  by  the  curved  oblique  shock  wave. 

The  overall  numerical  algorithm  is  discussed  in  detail  in  Section  V. 
A  brief  introduction  to  the  overall  numerical  algorithm  is  presented  here 
to  identify  the  various  flowfield  models  that  are  required  to  determine 
the  performance  of  a  scarfed  propulsive  nozzle. 

The  flow  in  the  throat  region  of  the  nozzle  is  assumed  to  be  com¬ 
pletely  specified  by  a  perturbation  analysis  that  depends  o..’y  on  the 
geometry  of  the  nozzle  throat.  From  this  analysis,  the  flow  properties 
along  an  initial-value  line  spanning  the  nozzle  throat  are  determined,  as 
illustrated  in  Figure  4. 

Right-running  Mach  lines  are  then  emanated  from  each  point  along  the 
initial-value  line,  starting  with  the  points  adjacent  to  the  axis,  and 
continued  until  they  intersect  the  nozzle  axis.  The  last  such  right-run¬ 
ning  Mach  line  emanating  from  the  nozzle  throat  wall  point  defines  the 
extent  of  the  initial -value  problem.  This  region  is  illustrated  in 
Figure  5. 

Point  locations  are  prespecified  along  the  circular  arc  initial 
expansion  contour  downstream  of  the  nozzle  throat  wall  point.  A  right¬ 
running  Mach  line  is  emanated  from  each  prespecified  wall  point  and  con¬ 
tinued  until  it  intersects  the  nozzle  axis.  The  last  such  right-running 
Mach  line  defines  the  extent  of  the  initial  expansion  flowfield.  This 
region  is  illustrated  in  Figure  6. 
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Left-running  Mach  lines  are  then  emanated  from  each  point  along  the 
last  right-running  Mach  line  in  the  initial  expansion  flowfield  and  con¬ 
tinued  until  the,y  intersect  the  nozzle  wall.  The  first  such  left-running 
Mach  line  begins  at  the  first  interior  point  on  the  aforementioned  right- 
running  Mach  line,  the  second  left-running  Mach  line  begins  from  the 
second  interior  point,  and  so  on.  This  process  is  continued  until 
either  the  end  of  the  basic  nozzle  is  reached,  or  the  last  point  on  the 
right-running  Mach  line,  which  lies  on  the  nozzle  axis,  has  been  reached. 
In  this  latter  case,  left-running  Mach  lines  are  emanated  from  the 
nozzle  axis  and  continued  until  they  intersect  the  nozzle  wall.  In 
either  case,  a  left-running  Mach  line  will  eventually  reach  the  end  of 
the  basic  nozzle  contour,  point  E.  This  flowfield  is  illustrated  in 
Figure  7. 

Due  to  the  finite  change  in  wall  slope  at  point  E,  an  oblique  shock 
wave  forms  at  that  point.  That  right-running  shock  wave  propagates  into 
the  flowfield  and  intersects  the  left-running  Mach  lines  as  they  propa¬ 
gate  toward  the  wall.  The  flowfield  up  to  the  shock  wave  is  irrotation- 
al,  ar.d  the  flowfield  downstream  of  the  shock  wave  is  rotational.  The 
flow  property  changes  across  the  shock  wave  depend'  on  the  strength  of 
the  shock  wave,  which  is  determined  by  both  the  upstream  flowfield  and 
the  downstream  flowfield.  After  passing  through  the  oblique  shock  wave, 
the  left-running  Mach  lines  continue  until  they  intersect  the  scarfed 
nozzle  extension.  This  process  is  continued  until  the  end  of  the  scarfed 
nozzle  extension,  point  F,  is  reached.  This  flowfield  is  illustrated  in 
Figure  8, 


/////// 


Once  the  entire  scarfed  nozzle  flowfield  has  been  calculated  as 
described  above,  the  performance  of  the  scarfed  nozzle  can  be  computed. 
That  computation  is  described  in  Section  IV. 
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2.  TRANSONIC  FLOW  MODEL 


A.  FLQWFIELD 

The  flowfield  in  the  transonic  throat  region  of  the  basic  nozzle 
is  assumed  to  be  completely  determined  by  the  geometry  of  the  nozzle 
throat.  The  influence  of  the  upstream  subsonic  geometry  is  assumed  to 
be  negligible.  This  is  a  good  approximation  when  the  nozzle  inlet  is 
"well  behaved."  Figure  9  illustrates  the  throat  geometry. 

As  illustrated  in  Figure  9.  the  sonic  line  originates  on  the 
nozzle  wall  upstream  of  the  throat  (point  T),  spans  the  throat  region, 
and  intersects  the  nozzle  axis  downstream  of  the  throat.  When  the 
throat  upstream  radius  of  curvature  ptu  is  large  compared  to  the  throat 
radius  y^,  linearized  flow  analyses  give  reasonable  predictions  for  the 
transonic  flowfield.  Zucrow  and  Hoffman  (1)  present  a  discussion  of 
such  analyses. 

In  the  present  case,  the  analysis  developed  by  Kliegel  and  Levine 
(2)  is  employed.  The  velocity  components  u  and  v  in  the  x  and  y 
directions,  respectively,  are  given  by 
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Figure  9.  Nozzle  throat  geometry 
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a*  =  (yRT*)’/2  =  (•^r)1/2  (22) 

R  =  ftu  ;-f2R)V2»  -  J_ 

yt  ’  h+1  /  yt’ y  yt  (23) 

Tt  Is  the  stagnation  temperature,  y  and  R  are  the  gas  specific  heat 
ratio  and  gas  constant,  respectively,  and  u-j(x,y),  etc.,  are  polynomial 
expressions  In  x  and  y  given  by  equations  (15.100)  to  (15.105)  In 
Reference  (1). 

In  the  present  analysis,  the  initial-value  line  is  a  parabola 
extending  from  the  nozzle  throat  wall  point,  point  T  in  Figure  9,  to 
the  point  on  the  nozzle  axis  where  M  a  1.0,  as  Illustrated  in  Figure  9. 
The  particular  parabola  employed  is  the  curve  along  which  the  flow 
angle  9  [9  *  tan'^v/u)]  is  equal  to  zero  in  Sauer's  linearized  flow 
analysis  [Reference  (3),  sea  also  Reference  (1)].  That  parabola  is 
specified  by 


x  a  e  - 


213+6] 


where  6*0  denotes  planar  flow,  6  =  1  denotes  axi symmetric  flow,  and 
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The  Initial-value  line  is  determined  by  specifying  the  y 
coordinates  of  NI  points  across  the  throat,  and  calculating  the 
corresponding  values  of  x  from  equation  (24).  Figure  10111  us trates  a 
typical  result.  These  sets  of  (x,y)  pairs  are  substituted  into  equa¬ 
tions  (20)  and  (21)  to  determine  the  velocity  components  u(x,y)  and 
v(x,y).  Then  the  velocity  magnitude  V  is  determined  from 

V  «  (u2  +  v2)1/2  (27) 


From  the  energy  equation  for  the  adiabatic  flow  of  a  thermally  and 
calorically  perfect  gas  [see  Reference  (4)], 
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where  cp  is  the  constant  pressure  specific  heat.  The  speed  of  sound  a 
is  given  by 


a  »  (yRT),/2 


(29) 


and  the  Mach  number  H  is  defined  as 

N  »  i  (30) 

a 

From  the  definition  of  the  stagnation  pressure  P^.  for  a  thermally  and 
calorically  perfect  gas  [see  Reference  (4)3 » 
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Figure  10.  Supersonic  initial -value  line. 
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From  the  thermal  equation  of  state,  the  density  p  is  given  by 


Thus,  all  the  flow  properties  along  the  initial -value  line  have  been 


determined. 


B.  PERFORMANCE  PARAMETERS 


The  nozzle  mass  flow  rate  ifi  is  obtained  by  numerically  integrating 


the  differential  mass  flow  rate  dift  across  the  initial -value  line  (IVL). 


Thus, 


diti  =  pu2irydy  -  pv2uydx 


iii  =  I  drti 
^  IVL 


The  nozzle  discharge  coefficient  is  defined  as 


r  = 

0  ifVo 


where  the  reference  ideal  one-dimensional  isentropic  choked  mass  flow 


rate  m^0  is  given  by  [see  Reference  (4)] 
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where  the  parameter  r  is  defined  as 


r-,feA 


The  Initial -value  line  thrust  FjVL  is  the  sum  of  the  axial 
components  of  the  pressure  forces  and  the  momentum  flux  across  the 
initial-value  line.  Thus, 


dF  =  (P  -  P J  2iry  dy  +  u  dift 

a 


■L 


F=  I  dF 
yIVL 


where  is  vhe  atmospheric  pressure.  The  thrust  efficiency  Tip  is 


defined  as 


F  F, 


where  the  reference  ideal  one-dimensional  isentropic  choked  thrust  is 
given  by 


F,-o  ■  <p*  -  V  At  +  a*  *1-0 


where  the  throat  sonic  velocity  a*  is  given  by  equation  (22)  and  the 


throat  choking  pressure  P*  Is  given  by  [see  Reference  (4)] 


r  ■  pt 


The  initial-value  line  specific  impulse  I  is  defined  as 


i  .s. 

sp  m 


The  specific  impulse  efficiency  ri|  is  defined  as 


sp,l-D 


where  the  reference  ideal  one-dimensional  specific  impulse  is  defined  as 


lsp,l-D 


By  combining  equations  (35),  (40),  and  (44),  it  is  obvious  that 


r'F  °  C0  nI 


The  above  analysis  completely  determines  all  of  the  properties  of 
interest  along  the  supersonic  initial -value  line.  Some  numerical 
examples  are  presented  in  Reference  (1). 
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3..  IRROTATIONAL  FLOW  MODEL 


The  portion  of  the  supersonic  fiowfield  between  the  supersonic 


initial-value  line  and  the  oblique  shock  wave  emanating  from  the  wall 


discontinuity  where  the  basic  nozzle  and  the  nozzle  extension  meet 
(i.e,,  point  E)  is  an  irrotational  fiowfield.  The  basic  equations 
governing  an  adiabatic  inviscid  flow  are  the  continuity  equation,  the 
Euler  momentum  equations,  and  the  energy  equation.  These  equations 
are  derived  in  Reference  (5).  For  steady  two-dimensional  flow,  those 


equations  are 


bux  +  pvy  +  UPX  +  Vpy  +  “  »  0 
Puux  +  Pvuy  +  Px  s  0 

PUVX  +  PWy  +  Py  “  0 

2  2 

utlx  +  vfly  +  u(r>x  +  v^y  “  0 


(47) 

(48) 

(49) 

(50) 


where  6  =  0  for  planar  flow  and  5  =  1  for  axisymaetrlc  flow,  and  b  is 
the  enthalpy  of  the  fluid.  For  isentropic  (i.e.,  adiabatic  and 
frictionless)  flows,  it  is  convenient  to  replace  the  energy  equation, 
equation  <50),  by  the  speed  of  sound  equation  [see  Reference  (5)]  to 
eliminate  the  enthalpy  h  from  the  system  of  equations.  That  equation 
is 
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uPx  +  vPy  ‘  a2(upx  +  vpy>  =  0  <51 ) 

Combining  equations  (47)  and  (51)  to  eliminate  the  derivatives  of 
density  yields 

w2(ux  +  vy  +  ~)  +  uPx  +  vPy  =  0  (52) 

Equations  (47),  (48),  (49),  and  (52)  comprise  a  set  of  four  coup¬ 
led  partial  differential  equations  for  determining  the  four  flow 
properties  u,  v,  P,  and  p.  Those  equations  are  the  basis  of  the 
rotational  flow  model  employed  downstream  of  the  oblique  shock  wave. 

That  flow  model  is  developed  further  in  Section  I I I. 5. 

If  the  flowfield  originates  in  a  region  of  parallel  uniform  flow 
and  remains  isentropic,  then  the  vorticity  of  the  flowfield  is 
everywhere  zero  [see  Reference  (5)].  Thus, 

t  =  ?xV  =  0  (53) 

which,  for  two-dimensional  flow,  gives 

uy  -  v,  =  G  (54) 

A  single  equation  containing  derivatives  only  of  u  and  v  can  be 
obtained  by  multiplying  equations  (48)  and  (49)  by  u  and  v,  respect¬ 
ively,  adding  them  together,  and  subtracting  equation  (52).  Th;  result 
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p(u  -a  )UX  +  p(vz  -  a2)vy  +  puv(uy  4.*x)  -  Ml  =  o  (55) 

Dividing  by  the  density  p  and  introducing  equation  (54)  gives 

2 

(u2  -  a2)ux  +  (v2  -  a2)vy  +  2uvuy  -  M-V-  =  o  (56) 

Equation  (56)  is  sometimes  called  the  gasdynamic  equation.  When 
considered  in  conjunction  with  equation  (54),  the  irrotationallty 
condition,  a  set  of  two  partial  differential  equations  is  obtained  for 
determining  the  two  velocity  components  u  and  v.  The  remaining  flow 
properties  (i.e.,  V,  T,  P,  p,  and  M)  can  be  determined  as  discussed  in 
Section  III. 2  for  transonic  flow. 

The  flow  model  employed  for  the  irrotational  flowfield  between  the 
nozzle  throat  and  the  oblique  shock  wave  consists  of  equations  (j6)  and 
(54),  which  are  repeated  and  renumbered  below. 

(u2  -  +  {v2  -  a2)vy  +  2uvuy  -  *  0  (57) 

uv  "  vx  =  0  (58) 


B.  METHOD  OF  CHARACTERISTICS 
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equations  combine  to  form  ordinary  differential  equations.  That  Is, 
the  partial  differential  equations  become  directional  derivatives  along 
the  characteristic  curves.  The  basic  principles  of  the  method  of 
characteristics  are  discussed  in  Reference  (6).  Application  of  the 
method  of  characteristics  to  steady  two-dimensional  irrotational  super¬ 
sonic  flow  is  discussed  in  Reference  (7).  A  brief  summary  of  those 
results  is  presented  below. 

C.  CHARACTERISTIC  EQUATIONS 

For  steady  two-dimensional  irrotational  supersonic  flow,  two 
families  of  characteristic  curves  exist;  the  left-running  and  right- 
rurmirig  Mach  lines.  The  slope  X  of  a  characteristic  curve  is  defined 
as 
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A  quadratic  equation  can  be  found  for  determining  X.  As  shown  in 
Reference  (6;,  that  equation  is 


(59) 
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(a 2  -  a2)X2  -  2uvX  +  (v2  -  a2)  «  0 


(60) 
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Solving  equation  (60)  for  X  yields 


2  jf} 

^  „  /dyN  _  uv  ±  avM  -1 
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u  -  a 


Equation  (61)  can  be  simplified  using  the  definitions 


(61) 
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u  =  V  cose,  v  =  V  sine,  6  =  tan  (v/u),  a  =  sin”  (1/M)  (62) 


where  6  is  the  flow  angle  and  a  is  the  Mach  angle.  After  considerable 
manipulation,  the  following  result  is  obtained. 

(^■)±  =  X+  =  tan  (6  ±  a)  (63) 

Equation  (63),  which  defines  the  characteristic  curves,  is  called  the 
characteristic  equation.  The  subscript  +  denotes  the  C+  characteristic, 
or  i eft-running  Mach  line,  so  called  because  it  runs  off  to  the  left  of 
the  streamline  when  looking  in  the  downstream  direction.  Similarly, 
the  subscript  -  denotes  the  C_  characteristic,  or  right-running  Mach 
line.  Figure  11  illustrates  these  characteristic  curves,  or  Mach  lines, 
at  a  point  in  a  flowfield.  Two  Mach  lines  pass  through  every  point  in 
a  flowfield,  resulting  in  two  infinite  families  of  left-running  and 
right-running  Mach  lines. 

D.  COMPATIBILITY  RELATIONS 

Equations  (57)  and  (58),  when  combined  and  written  as  directional 
derivatives  along  the  Mach  lines,  yield  a  single  ordinary  differential 
equation,  called  the  compatibility  relation,  which  is  valid  along  each 
family  of  Mach  lines.  Consider  any  continuous  function  f(x,y).  The 
total  derivative  of  f  is  given  by 

df  8  f  dx  +  f  dy  (64) 


which  may  be  written  as 


lrf*  +  fy&=f*  +  xf, 


The  partial  derivatives  In  equations  (57)  and  (58)  can  be  put  In 
the  form  of  equation  (65)  by  multiplying  equation  (58)  by  the  term 


(u2  -  a2)  X  -  2uv 


and  adding  the  result  to  equation  (57)  to  obtain 


(u2  -  a2)ux  +  (v2  -  a2)vy  +  2uvuy  -  + 


C(u2  -a2)x  -  2uv](u  -v  )  =  0 

y  * 


Rearranging  equation  (67)  yields 


(u2  -  a2)[ux  +  Xuy]  +  [2uv  -  (u2  -  a2)X]  vx  + 


(v2  -  a2)v  -  ~ r 1  *  0 


y  y 


From  equation  (60), 


(v2-a2)  -  [2iiv  -  (u2-a2)X]  X 


( 


Substituting  equation  (69)  into  equation  (68)  gives 


(70) 
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(u2  -  a2)  Lux  +  Xuy]  +  [2uv  -  (u2  -  a2)X][vx  +  X  y]  - 


The  partial  derivatives  in  equation  (70)  are  in  the  form  of  equation 
(65).  Consequently,  equation  (70)  may  be  written  as 
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(u2-a2)du  +  [2uv  -  (u2-a2)X]dv  -  ^-^dx  =  0 


(71) 


Equation  (71)  is  a  total  differential  equation  which  is  valid 
along  the  Mach  lines.  It  is  called  the  compatibility  relation.  When 
applied  along  both  left-running  and  right-running  Mach  lines,  equation 
(71)  yields  two  equations  for  determining  changes  in  u  and  v  along  the 
two  families  of  Mach  lines. 

E.  SUMMARY  OF  RESULTS 

In  summary,  for  steady  two-dimensional  irrotational  supersonic 
flow,  the  characteristic  equations  are 
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)+  =  X+  =  tan  (6  ±  o) 


and  the  compatibility  relations  are 


(u2-a2)  dut  +  [2uv  -  (u‘-a2)  X±3dv±  -  dx±  =  0 


(72) 


(73) 
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F.  NUMERICAL  METHOD  OF  CHARACTERISTICS 
The  characteristic  equations  and  compatibility  relations  for  steady 
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two-dimensional  irrotational  supersonic  flow  are  presented  in  the  pre¬ 
ceding  section.  Those  equations  are  ordinary  differential  equations. 
They  can  be  integrated  numerically  by  the  second-order  accurate 
modified-Euler  predictor-corrector  method.  Such  a  procedure  is  called 
the  numerical  method  of  characteristics.  The  details  of  its  implemen¬ 
tation  are  presented  in  References  (6)  and  (7).  That  is  the  procedure 
employed  to  calculate  the  flowfield  between  the  throat  of  the  basic 
nozzle  and  the  oblique  shock  wave  at  the  entrance  to  the  scarfed 
extension. 

The  basic  procedure  is  to  start  from  a  supersonic  initial-value 
line  across  the  nozzle  throat,  such  as  illustrated  in  Figure  10,  The 
left-running  and  right-running  Mach  lines  emanating  from  the  Initial - 
value  line  points  are  constructed  by  solving  equation  (72).  The 
velocity  components  u  and  v  are  determined  at  the  intersections  of  the 
Mach  lines  by  solving  equation  (73).  Since  equations  (72)  and  (73)  are 
coupled  and  nonlinear,  they  must  be  solved  numerically  and  simultaneous 
ly.  This  procedure  defines  *he  interior  point  unit  process.  When  a 
left-running  Mach  line  impinges  on  .he  nozzle  wall,  the  boundary 
conditions  of  known  wall  location  and  slope  must  be  applied.  A  right- 
running  Mach  line  is  reflee tec  into  the  flowfield  from  that  point  of 
impingement.  That  procedure  defines  the  -vail  point  unit  process. 
Similarly,  when  a  right-running  Mach  line  reaches  the  nozzle  axis,  the 
boundary  conditions  of  y  =  v  =  0  must  be  applied.  A  left-running  Mach 
line  is  reflected  into  the  flowfield  from  that  point  of  impingement. 
That  procedure  defines  the  axis  point  unit  process. 

This  procedure  is  continued  until  the  entire  flowfield  between  the 
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nozzle  throat  and  the  oblique  shock  wave  has  been  crisscrossed  with 
left-running  and  right-running  Mach  lines  and  the  flow  properties  have 
been  calculated  at  all  the  intersections  of  the  two  families  of  Mach 
lines.  A  coarse  schematic  of  the  resulting  Mach  line  network  Is 
illustrated  in  Figure  12.  Figure  12(a)  illustrates  the  Mach  line 
network  when  right-running  Mach  lines  are  used  to  construct  the  flow- 
field.  Figure  12(b)  illustrates  the  Mach  line  network  when  left-running 
Mach  lines  are  used  to  construct  the  flowfield. 


4.  SHOCK  WAVE  MODEL 


A.  PHYSICAL  MODEL 

The  nozzle  geometric  model  considered  in  the  present  investigation 
consists  of  a  basic  nozzle  followed  by  a  nozzle  extension,  as  illus¬ 
trated  in  Figure  13.  As  discussed  in  Section  II. 1,  the  basic  nozzle  may 
be  conical,  quadratic,  or  tabular.  The  nozzle  extension  is  conical. 

At  the  junction  of  the  basic  nozzle  and  the  nozzle  extension,  the 
slope  of  the  wall  of  the  basic  nozzle  may  be  larger  than,  equal  to,  or 
less  than  the  slope  of  the  nozzle  extension,  resulting  in  the  gener¬ 
ation  of  an  oblique  shock  wave,  a  continuous  flow,  or  a  centered 
expansion  fan,  respectively,  as  illustrated  in  Figure  13.  The  third 
case,  where  the  slope  of  the  basic  nozzle  is  less  than  the  slope  of  the 
nozzle  extension,  is  not  considered  in  the  present  study.  The  second 
case,  where  the  slope  of  the  basic  nozzle  is  equal  to  the  slope  of  the 
nozzle  extension  and  the  flow  is  continuous,  is  treated  as  the  special 
case  of  an  oulique  shock  wave  of  infinitesimal  stregth.  Consequently, 
in  the  present  analysis,  it  is  always  assumed  that  an  oblique  shock 
wave  emanates  from  the  junction  of  the  basic  nozzle  and  the  nozzle 
extension.  The  flow  model  for  determining  the  location  of  and  the  flow 
properties  across  the  oblique  shock  wave  are  discussed  in  this  section. 

The  strength  of  an  oblique  shock  wave  is  always  greater  than  a 
continuous  flow.  Consequently,  the  oblique  shock  emanating  from  the 
junction  point  will  propagate  upstream  into  the  flowfield  determined  by 
the  basic  nozzle,  as  Illustrated  in  Figure  14.  The  strength  of  the 
oblique  shock  wave  is  influenced  by  both  the  upstream  flowfield 
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determined  by  the  basic  nozzle  and  the  downstream  flowfield  determined 
by  the  nozzle  extension. 

B.  GOVERNING  EQUATIONS 

The  basic  geometry  of  an  oblique  shock  wave  Is  illustrated  In 
Figure  15.  The  upstream  velocity  is  denoted  by  Vj,  the  flow  turning 
angle  between  the  upstream  and  downstream  flow  directions  is  denoted  by 
6,  and  the  downstream  velocity  is  denoted  by  Vg.  The  shock  wave  angle 
is  denoted  by  s.  All  of  the  flow  properties  (i.e.,  V,  H,  P,  T,  p,  P^., 
and  Tt)  are  known  upstream  of  tue  oblique  shock  wave.  The  problem  is 
to  determine  the  corresponding  flow  properties  downstream  of  the  shock 
wave  for  a  specified  value  of  the  flow  turning  angle  6.  A  detailed 
analysis  of  the  oblique  shock  wave  is  presented  in  Reference  (8). 

A  coordinate  system  normal  and  tangential  to  the  oblique  shock 
wave  is  employed.  The  equations  of  continuity,  momentum,  and  energy 
are  applied  in  that  coordinate  system.  A  major  conclusion  of  that 
analysis  is  that  the  tangential  component  of  velocity  does  net  change 
across  the  oblique  shock  wave  (i.e.,  *  VT2  *  VT),  and  that  the 

components  of  velocity  normal  to  the  shock  wave  (1  .e, ,  and  V^}  are 
governed  by  the  equations  for  a  normal  shock  wave. 

If  the  wave  angle  e  is  known,  the  upstream  normal  Kach  number  Is 
given  by 


sir.  e 
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From  Figure  IS,  the  downstream  normal  Nach  number  is  given  by  [see 
Reference  (3)3 


(75) 


m 

Pjfl 


m 

M 

m 


few 

M 

m 


hz 


From  Figure  10,  the  downstream  Mach  number  is 
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For  a  given  wave  angle  e,  the  corresponding  flow  turning  angle  6  is 
given  by  [see  Reference  (8)] 


1  _  fytl  M1 

tan  6  17  MN1  ‘  1 


tan  e 


(77) 


If  the  flow  turning  angle  6  is  specified  instead  of  the  wave  angle  e, 
then  equations  (74)  to  (77)  must  be  solved  by  iteration  (e.g.,  using 
the  secant  method).  In  either  case,  for  a  given  a..J  either  e  or  6, 
the  corresponding  M^,  Mg,  and  either  e  or  6  can  be  determined. 

The  remaining  flow  property  ratios  are  obtained  from  the  corres¬ 
ponding  results  for  a  normal  shock  wave  [see  Reference  (8)]. 
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v2  =  (Vfc  +  Vf) 
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C.  NUMERICAL  SOLUTION  PROCEDURE 
Two  different  oblique  shock  wave  situations  occur  in  a  scarfed 
nozzle.  The  first  situation  is  at  the  junction  of  the  basic  nozzle  and 
the  nozzle  extension,  point  E  in  Figure  14,  where  the  oblique  shock  wave 
emanates.  At  that  point,  the  flow  turning  angle  6  is  specified  as  the 
difference  between  the  flow  angle  at  the  end  of  the  basic  nozzle  and 
the  flow  angle  at  the  beginning  of  the  nozzle  extension.  The  procedure 
described  in  the  preceding  paragraphs  can  be  apolied  directly  to  deter¬ 
mine  the  properties  of  the  oblique  shock  wave  at  that  point. 

The  second  situation  is  the  general  point  on  the  oblique  shock 
wave  where  both  the  upstream  and  downstream  flowfield  affect  the  shock 
wave.  Since  the  strength  of  an  oblique  shock  wave  is  greater  than  the 
strength  of  a  Mach  line,  the  oblique  shock  wave  emanating  from  point  E 
propagates  into  the  upstream  flowfield  at  a  steeper  angle  than  the  angle 
of  the  right-running  Mach  lines  emanating  from  the  basic  nozzle  contour. 
Since  a  right-running  oblique  shock  wave  turns  the  flow  toward  the  axial 
direction  { i  .e . 5  toward  the  shock  wave)  and  the  Mach  number  and  Mach 
angle  are  smaller  downstream  of  an  oblique  shock  wave,  the  slopes  of  the 
right-running  Mach  lines  on  the  downstream  side  of  the  shock  wave  are 
greater  that  the  slope  of  the  shock  wave.  Consequently,  the  location 
and  stength  of  the  oblique  shock  wave  are  also  influenced  by  the  down¬ 
stream  flowfield.  Thus,  the  portion  of  the  downstream  flowfield  that 
interacts  with  the  oblique  shock  wave  must  be  determined  along  with 
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the  determination  of  the  shock  wave  itself.  Figure  16  illustrates  the 
propagation  of  the  oblique  shock  wave  into  the  upstream  flowfield  and 
the  overtaking  of  the  shock  wave  on  the  downstream  side  by  a  Mach  line 
from  the  nozzle  extension. 

The  region  of  the  flowfield  affecting  the  solution  at  a  point  on 
the  oblique  shock  wave,  denoted  as  point  4,  is  outlined  in  Figure  16. 

The  upstream  flowfield  is  assumed  to  be  specified,  and  the  propagation 
of  the  oblique  shock  wave  from  point  E  to  point  5  has  been  calculated, 
as  has  the  right-running  Mach  line  on  the  downstream  side  of  the  shock 
wave  from  the  nozzle  extension  to  the  left-running  Mach  line  from  point  5. 
An  enlarged  illustration  of  this  region  is  presented  in  Figure  17. 

The  determination  of  the  location  of  point  4  and  the  flow  properties 
on  either  side  of  the  oblique  shock  wave  at  that  point  requires  an  iter¬ 
ative  procedure  having  five  major  steps.  The  first  step  is  to  find  the 
location  of  point  4  as  the  intersection  of  the  oblique  shock  propagated 
from  point  5  at  the  average  wave  angle  e  with  the  left-running  Mach  line 
propagated  from  point  2  in  the  known  upstream  flowfield.  This  is  accom¬ 
plished  in  an  iterative  manner  by  first  assuming  a  value  for  c^  (note 
that  is  known),  performing  the  remaining  four  major  steps  outlined 
below,  and  then  varying  e ^  iteratively  until  the  overall  procedure  con¬ 
verges. 

The  second  step  is  the  application  of  the  method  of  characteristics 
for  irrotational  flow  (see  Section  III. 3)  to  determine  the  flow  proper¬ 
ties  at  point  4  on  the  upstream  side  of  the  oblique  shock  wave.  The 
third  step  is  the  determination  of  the  flow  properties  at  point  4  on  the 
downstream  side  of  the  oblique  shock  wave  by  applying  the  oblique  shock 
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wave,  analysis  presented  at  the  beginning  of  this  section.  The  fourth 
step  is  the  application  of  the  method  of  characteristics  for  rotational 
flow  (see  Section  III. 5)  to  determine  the  flow  properties  at  point  4', 
which  is  the  intersection  of  the  left-running  Mach  line  from  point  4  on 
the  downstream  side  of  the  oblique  shock  wave  with  the  right-running 
Mach  line  from  point  V. 

All  of  the  above  four  steps  are  performed  for  a  specified  value  of 
e4»  the  oblique  shock  wave  angle  at  point  4,  which  together  with  the 
known  wave  angle  e5  at  point  5,  defines  the  average  wave  angle  e  and 
thus  the  location  of  point  4.  A  check  on  the  correctness  of  is  made 
as  the  fifth  step  in  the  overall  procedure.  The  pressure  P4d  on  the 
downstream  side  of  the  oblique  shock  wave  Is  computed  in  step  3  above. 

A  right-running  Mach  line  is  projected  rearward  from  point  4  to  intersect 
a  previous  Mach  line  (or  the  wall  of  the  nozzle  extension  for  the  first 
point  on  the  oblique  shock  wave  downstream  of  point  E)  at  point  3'.  The 
compatibility  relation  valid  along  that  Mach  line  is  solved  to  obtain  a 
second  value  for  P4fJ.  The  discrepancy  detween  these  two  values  of  P4d 
is  a  measure  of  the  error  in  e^.  The  secant  method  can  be  used  to  vary 
e4  to  drive  the  difference  between  the  two  values  of  to  within  any 
desired  convergence  limit. 

The  procedure  described  above  is  repeated  at  successive  points 
along  the  oblique  shock  wave  until  the  s**ock  wave  approaches  the  nozzle 


5.  ROTATIONAL  FLOW  MODEL 

A.  GOVERNING  EQUATIONS 

The  portion  of  the  flowfield  downstream  of  the  oblique  shock  wave 
emanating  from  the  junction  of  the  basic  noz2le  and  the  nozzle  exten¬ 
sion  is  a  rotational  flowfield  due  to  the  transverse  entropy  gradient 
caused  by  the  curved  oblique  shock  wave.  The  basic  equations  governing 
a  rotational  isentropic  flowfield  are  presented  in  Section  III. 3  [equa¬ 
tions  (47),  (48),  (49),  and  (51)].  Those  equations  are  repeated  and 
renumbered  below. 


,  ,  .  6Pv  _  n 

pux  pvy  upx  vpy  +  ~  "  0 


(83) 


puux  +  PVUy  *■  Px  =  0 


(84) 


puvx  +  pvVy  +  Py  =  0 


(85) 


uPx  +  vPy  "  a  +  vpy^  "  0 


(86) 


Due  to  the  presence  of  an  entropy  gradient  Vs,  the  vorticity  C  is 
not  zero,  and  equation  (53)  is  not  valid.  Consequently,  the  simplifi¬ 
cations  obtained  in  equations  (57)  and  (58)  are  not  applicable,  and  the 
full  set  of  equations  presented  above,  equations  (83)  to  (86),  must  be 
solved.  Those  equations  are  solved  by  the  numerical  method  of 
characteristics  (see  Section  III.3.B). 

G.  CHARACTERISTIC  EQUATIONS 

For  steady  two-dimensional  rotational  supersonic  flow,  three 
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families  of  characteristic  curves  exist;  the  left-running  and  right- 
running  Mach  lines  and  the  streamlines.  The  slope  X  of  a  character¬ 
istic  curve  is  defined  as 


d i  -  x 
dx  X 
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A  fourth-order  equation  can  be  found  for  determining  X.  As  shown  in 
Reference  (9),  that  equation  is 


(uX  -  v)2  [(uX  -  v)2  -  a2(l  +  X2)]  =  0 


The  fourth-order  equation  is  the  product  of  two  second-order  equations. 
Solving  either  of  those  two  second-order  equations  gives  two  of  the 
characteristic  slopes,  and  solving  the  other  second-order  equation 
gives  the  other  two  characteristic  slopes. 

Solving  the  first  second-order  term  yields 


(uX  -  vr  =  0 
o 


X  =  (&)  «  1 
o  Wo  u 


which  defines  the  streamline.  Since  Xq  appears  two  times,  the  stream¬ 
line  is  a  repeated  characteristic.  The  subscript  o  will  be  used 
henceforth  to  denote  the  streamline. 

Rearranging  the  second  second-order  term  yields 
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Equation  (91)  Is  Identical  to  equation  (60) ,  which  defines  the  Mach 
lines.  As  shown  In  Section  III.3.C,  the  Mach  lines  are  given  by 


(jj£)±  s  X±  «  tan  (0  ±  a)  (92) 

The  subscripts  ±  are  used  to  denote  the  left-running  and  right-running 
Mach  lines,  respectively. 

Consequently,  for  steady  two-dimensional  supersonic  rotational 
flow,  four  characteristic  curves  are  obtained:  the  streamline  repeated 
twice  and  the  left-running  and  right-running  Mach  lines.  Figure  18 
illustrates  these  characteristic  curves  at  a  point  in  a  flowfleld.  A 
streamline  and  two  Mach  lines  pass  through  every  point  in  a  flowfleld, 
resulting  in  three  infinite  families  of  characteristic  curves. 

C.  COMPATIBILITY  RELATIONS 

Equations  (83)  to  (86),  when  combined  and  written  as  directional 
derivatives  along  the  streamlines  and  Mach  lines,  yield  two  ordinary 
differential  equations  valid  along  the  streamlines  and  one  ordinary 
differential  equation  valid  along  each  Mach  line.  Thus,  four  compatf 
bility  relations  are  determined  for  steady  two-dimensional  supersonic 
rotational  flow. 

As  shown  In  Section  III.3.D,  the  total  derivative  of  any  contin¬ 
uous  function  f(x,y)  must  satisfy  the  relation  [see  equation  (65)] 


Figure  18.  Characteristic  curves  for  rotational  flow. 
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The  partial  derivatives  In  equations  (83)  to  (86)  can  be  put  .In  the 
form  of  equation  (93)  by  forming  appropriate  linear  combinations  of 
those  equations. 

First,  consider  the  streamline  where  Xq  satisfies  equations  (89) 
and  (90).  Equation  (86)  may  be  rearranged  as 

U(PX  +  |  Py)  -  Ufi2(PX  +  l  Py)  8  0  (94) 


The  derivatives  In  equation  (94)  are  in  the  form  of  equation  (93), 
where  X  -  XQ  =  v/u.  Consequently,  along  streamlines,  equation  (94)  may 
be  written  as 

dp  -  a2dp  e  0  (on  streamlines)  (95) 

Equation  (95)  is  the  speed  of  sound  equation. 

The  second  compatibility  relation  valid  along  streamlines  can  be 
determined  by  multiplying  equation  (85)  by  X  and  adding  the  result  to 
equation  (84) » 


pu(ux  +  uy)  +  Xpu(vx  +  ~  vy)  ♦  (Px  ♦  XPy)  *  0  (S6) 

The  derivatives  in  equation  (96)  are  in  the  form  of  equation  (93), 
where  X  =  =  v/u*  Consequently,  noting  that  along  streamlines  Xu  =  v, 

equation  (S6)  may  be  written  as 


pu  du  *  pv  dv  +  dP  =  0 


(97) 


Equation  (97)  Is  Bernoulli's  equation  for  steady  two-dimensional 
Isentropic  flow,  which  may  be  written  as 

pV  dV  +  dP  =  0  (on  streamlines)  (98) 

Second,  consider  the  Mach  lines  where  X±  satisfies  equations  (91) 
and  (92).  The  compatibility  relation  valid  along  Mach  lines  is  obtained 
by  forming  the  following  sum: 

(uX±  -  v)  [Eq.  (82)]  +  (-X±)  [Eq.  (80)]  +  [Eq.  (81)]  ♦ 

CEq.  (82)]  »  0  ;») 

a4 

After  considerable  manipulation,  equation  (99)  reduces  to 

(pv)du^  -  (pu)dvt  +  [X±  -  u(uX±  -  v)/a2]dP±  - 
d[pv(uXi  -  v)/y]dxt  a  0  (100) 


Equation  (100)  can  be  simplified  by  using  the  definitions 

u  -  V  cose,  v  =  V  sine,  0  «  tan’^v/u),  c  =  sin"1  (1/H)  (101) 

where  8  is  the  flow  angle  and  a  is  the  Kach  angle.  After  considerable 
manipulation,  the  following  result  is  obtained. 
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dx±  *  0  002) 


where  the  upper  subscripts  on  dP,  de,  and  dx  and  the  tipper  signs  In  ± 
de  and  cos  (6±a)  correspond  to  the  left-running  Kach  line*  and  vice 
versa. 

D.  SUfflflRV  OF  RESULTS 

In  stannary,  for  steady  two-dimensional  rotational  supersonic  flow* 
three  Infinite  families  of  characteristic  curves  exist;  the  streamlines 
and  the  left-running  and  right-running  Kach  lines.  Along  the  stream¬ 
lines. 


003) 


dP0  -  a2  do0  =  0 


OM) 


pV  dV0  +  <iP0  *  0 


{!S5) 


(106) 

(10?) 


E.  KilMERICPl  HETHOD  OF  CHARACTERISTICS 
The  characteristic  equations  and  compatibility  relations  for 
steady  two-dimensional  rotational  supersonic  f!cw  are  presented  in  the 


preceding  section.  Those  equations  are  ordinary  differential  equations. 
They  can  be  integrated  numerically  by  the  second-order  accurate 
modified-Euler  predictor-corrector  method.  Such  a  procedure  is  called 
the  numerical  method  of  characteristics.  The  details  of  Its  implemen¬ 
tation  are  presented  in  Reference  (9).  That  Is  the  procedure  employed 
to  calculate  the  flowfield  in  the  nozzle  extension  downstream  of  the 
oblique  shock  wave  emanating  from  the  junction  of  the  basic  nozzle  and 
the  nozzle  extension. 

The  basic  procedure  Is  to  start  from  the  downstream  side  of  the 
oblique  shock  wave.  The  left-running  and  right-running  Kach  »i  r.es  are 
constructed  by  solving  equation  (108).  The  pressure  P  and  flow  angle  6 
are  determined  at  the  intersections  of  the  Kach  lines  by  solving 
equation  (107).  The  streamlines  passing  through  these  intersection 
points  are  projected  upstream  to  intersect  diagonal  lines  joining 
upstream  intersection  points  by  solving  equation  (103).  The  density  p 
and  velocity  V  are  determined  at  the  intersection  points  by  solving 
equations  (104)  and  (105).  Since  equations  (103)  to  (107)  are  coupled 
and  nonlinear,  they  must  be  solved  numerically  and  simultaneously. 

This  procedure  defines  the  interior  point  unit  process. 

yfrsn  a  left-running  Kach  line  impinges  on  the  wall  of  the  nozzle 
extension,  the  boundary  conditions  of  known  wall  location  and  slope  must 
be  applied.  A  right-running  Kach  line  is  reflected  into  the  flowfield 
from  that  point  of  impingement.  That  procedure  derives  the  wall  point 
unit  process. 

This  procedure  is  continued  until  the  entire  flowfield  in  the 
nozzle  extension  has  been  crisscrossed  with  lefc-rumving  and  right- 


SECTION  IV 


SCARFED  NOZZLE  PERFORMANCE  MODEL 

1 .  INTRODUCTION 

The  performance  of  a  propulsive  nozzle  is  specified  by  its  mass 
flow  rate,  thrust  vector,  and  moment  vector.  For  conventional  nozzles 
the  thrust  vector  lies  along  the  nozzle  axis,  and  the  only  thrust  com¬ 
ponent  is  the  axial  component.  In  that  case,  all  of  the  moments  are 
zero.  For  a  scarfed  nozzle,  side  forces  and  moments  are  present.  The 
present  section  presents  the  performance  model  for  the  scarfed  nozzle 
configurations  considered  in  this  study  (see  Section  III. 2). 


2.  NOZZLE  GEOMETRIC  MODEL 


The  nozzle  geometric  model  considered  in  the  present  investigation 
is  illustrated  schematically  in  Figure  20.  The  basic  nozzle  consists  of 
a  double  circular  arc  threat  following  by  a  supersonic  expansion  con¬ 
tour.  The  supersonic  expansion  contour  may  be  conical,  quadratic,  or 
specified  by  tabular  data.  The  basic  nozzle  geometry  and  flowfield  are 
axi symmetric. 

The  nozzle  extension  is  a  conical  contour  attached  to  the  basic 
nozzle.  The  flow  in  the  nozzle  extension  is  assumed  to  be  axi symmetric. 
That  is,  the  wave  (i.e.,  Ma^h  line  or  shock  wave)  emanating  from  point 
E  Is  assumed  to  fall  downstream  of  the  exit  of  the  scarfed  conical 
extension,  as  illustrated  in  Figure  20.  In  that  case,  the  flowfield  in 
the  scarfed  nozzle  extension  may  be  computed  as  though  the  nozzle 
extension  were  not  scarfed  (i.e.,  as  if  the  dashed  protion  of  the  con¬ 
tour  were  present).  The  flowfield  downstream  of  line  EG,  which  is 
outside  of  the  solid  nozzle  boundary,  is  obviously  three-dimensional. 
However,  the  flowfield  upstream  of  line  EF  is  axisyrraietric. 


Scarfed  nozzle  extension 


Figure  20.  Scarfed  nozzle  geometric  model 


3.  BASIC  NOZZLE  PERFORMANCE  EVALUATION 

The  basic  nozzle  Is  a  conventional  axisymmetrlc  nozzle.  The 
flowfleld  In  the  throat  region  Is  predicted  by  the  linearized  flow 
analysis  presented  In  Section  111.2.  From  that  analysis*  an  initial- 
value  line  Is  obtained  from  which  the  supersonic  flowfield  can  be 
calculated.  In  addition,  (see  Figure  10),  the  nozzle  mass  flow  rate  ift, 
discharge  coefficient  Cp,  and  initial-value  line  thrust  Fjy^  are  obtain¬ 
ed. 


The  flowfleld  in  the  supersonic  region  is  calculated  by  the  method 
of  characteristics  for  steady  two-dimensional  Irrotatlonal  supersonic 
flow,  as  discussed  in  Section  III.3.  A  schematic  of  the  resulting  Mach 
line  network  is  presented  In  Figure  12.  From  that  analysis,  the  pressure 
acting  along  the  wall  Is  known  at  each  point  where  a  Mach  line  inter¬ 
sects  the  wall.  The  thrust  developed  along  the  supersonic  expansion 
contour  is  obtained  by  integrating  (numerically)  the  axial  component  of 
the  force  developed  by  the  pressure  acting  on  the  wall.  Thus, 


(P  -  Pj  2wy  dy 

a 


(108) 


The  total  thrust  developed  by  the  basic  nozzle,  FN,  is  the  sum  of  the 
thrust  developed  across  the  Initial-value  line,  F^,  and  the  thrust 
developed  along  the  supersonic  contour,  F^.  Thus, 


(109) 
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4.  SCARFED  NOZZLE  EXTENSION  PERFORMANCE  EVALUATION 

The  flowfield  in  the  nozzle  extension  is  calculated  by  the  method 
of  characteristics  for  steady  two-dimensional  rotational  supersonic 
flow,  as  discussed  in  Section  III. 4.  A  schematic  of  the  resulting  Mach 
line  network  is  presented  in  Figure  19.  From  that  analysis,  the 
pressure  acting  along  the  wall  of  the  nozzle  extension  is  known  at  each 
point  where  a  Mach  line  intersects  the  wall.  The  thrust  and  moments 
developed  by  the  scarfed  nozzle  extension  are  obtained  by  integrating 
(numerically)  the  differential  force  and  moment  components  developed  by 
the  pressure  acting  on  the  wall. 

In  the  present  analysis,  it  is  assumed  that  the  nozzle  extension  is 
axisymnetric  (in  fact,  conical)  and  that  the  exit  of  the  nozzle  exten¬ 
sion  is  scarfed  by  a  plane  perpendicular  to  the  xy  plane  which  passes 
through  point  E  on  the  top  of  the  nozzle  where  scarfing  begins  and 
through  point  F  on  the  bottom  of  the  nozzle  where  scarfing  ends.  Conse¬ 
quently,  the  flowfield  is  symmetrical  about  the  xy  plane.  Thus,  as 

illustrated  in  Figure  21,  only  two  force  components  exist,  F  and  F  , 

x  y 

and  only  one  moment  component  exists,  The  force  components  and 
moment  are  assumed  to  act  at  the  center  of  the  throat  of  the  nozzle. 

The  geometric  model  employed  to  determine  the  force  components  and 
moment  developed  by  the  scarfed  nozzle  extension  is  illustrated  in  Fig¬ 
ure  22.  The  geometry  of  the  scarfed  nozzle  extension  is  specified  by 
the  location  of  point  E  (xp,  ye),  the  angle  of  the  scarfed  conical 
extension  8.,  and  either  the  scarfing  angle  8  or  the  location  of  point 
F  (x^,  yf). 
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The  locations  of  points  F  and  F  are  determined  as  follows.  The 
equation  of  the  scarfing  line,  line  EF*  is  given  by 

y  Z 1*.  -tan  e  (HO) 

X  -  Xe  p 

where  the  angle  B  is  in  the  range  0  <  &  <  90  deg.  Solving  equation 
(110)  for  y  gives 

y  B  ♦  x*  tan  6  -  x  tan  B  »  a  -  x  tan  B  Oil) 

6  C 

where 


a  a  ye  +  xfitan  6 

012) 

The  equation  of  line  E'F*  is  given  by 

y  -  ye» 

— ■  -tan  8* 

X  -  *e,  f 

013} 

y  *  ye,  +  xe,tan  8f  -  xtan  8f  «  b  -  xtan  flf 

OH) 

where 

b  ■  ye,  +  xe.  tan  8f 

015) 

Substituting  (x^,  >y)  into  equations  (111)  and  (113)  gives 
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yf,  =  a  -  x^ tang 


(116) 


y.,  =  b  -  x.  tan  9. 


(117) 


Solving  equations  (116)  and  (117)  simultaneously  yields 


s  a  -  b 
xf  **  tan  6  -  tan  6, 


(118) 


Equation  (117)  may  then  be  solved  for  y^.  and  y^  =  -y^,. 

Point  G  illustrated  in  Figure  22  is  the  intersection  point  of  a 

plane  perpendicular  to  the  nozzle  axis  with  the  edge  of  the  scarfed 

extension,  line  EF.  The  location  of  point  G  corresponding  to  a  specified 

value  of  x  is  obtained  by  substituting  x  into  equation  (111).  Thus, 

9  9 


yg  =  a  -  xg  tan  S 


(119) 


The  angle  cU)  illustrated  in  Figure22  (b)  is  determined  as  follows. 


C(x)  *  sin''(yg/ywal,) 


(120) 


where  y  ^  is  the  radius  of  the  scarfed  conical  extension  at  point  G. 
The  angle  £(x)  illustrated  in  Figure  22(b)  is  detemined  from 


i'M 


C(x) 


021) 


All  of  the  geometric  properties  of  the  scaned  conical  extension 

are  now  determined.  Those  properties  are  the  coordinates  (x  ,  y  )  and 

0  e 

(x^»  yf).  the  equation  of  the  scarfing  line  EF'  [equation  (114)],  and 
the  angle  v(x). 
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The  performance  of  the  scarfed  conical  extension  Is  determined  by 
Integrating  the  pressure  forces  acting  along  the  wall  from  the  point 
where  the  scaflng  begins,  point  E,  to  the  end  of  the  scarfed  conical 
extension,  point  F.  The  geometric  model  for  Integrating  the  pressure 
forces  Is  presented  In  Figure  23.  The  scarfed  conical  extension  Is 
assumed  to  be  symmetrical  about  the  xy  plane.  Consequently,  the 
pressure  forces  acting  on  only  one-half  of  the  nozzle  must  be  deter¬ 
mined  by  integration.  The  total  pressure  force  Is  then  found  by 
symmetry.  For  symmetry  about  the  xy  plane,  no  net  component  of  the 
pressure  force  acts  in  the  z-di recti on. 

The  components  of  the  pressure  force  acting  in  the  x  and  y  direct¬ 
ions  are  determined  by  defining  a  differential  element  of  area  d3T 
having  a  magnitude  dA  and  acting  In  the  direction  of  the  outward  unit 
normal  rT.  The  unit  vector  system  T,  J,  F  in  the  Cartesian  coordinate 
system  xyz  is  illustrated  in  Figure  23.  The  radial  unit  vector  T  , 
which  is  directed  radially  outward  from  the  x-axis  and  lies  in  a  plane 
perpendicular  to  the  x-axis,  is  also  illustrated  in  Figure  23.  The 
element  of  area  dft  is  given  by 

dA  -  Tr  dAr  -  T  dAx  (122) 


where 


dA„  -  r  d9  dx 

r 


(123) 


dA  =  r  d8  dr 
x 


(124) 
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In  the  Cartesian  coordinate  system,  T  is  given  by 

Tj.  =  J  cos  0  +  k  sin  6  (125) 

where  8  is  measured  counterclockwise  from  the  positive  y-axis.  Combin¬ 
ing  equations  (122)  to  (125)  gives 

d£  »  -T  r  dO  dr  +  J  r  cos  0  d0  dx  +  k  r  sin  0  de  dx  (126) 

The  pressure  force  is  given  by 

d?p=(P-Pa)dA  (127) 

Due  to  the  symmetry  about  the  xy  plane,  the  net  component  of  the  pres¬ 
sure  force  perpendicular  to  that  plane  is  zero.  Consequently,  the 
third  term  on  the  right-hand  side  of  equation  (126)  may  be  dropped  from 
consideration,  and  equations  (126)  and  (127)  may  be  combined  to  yield 

dfp  3  -1(P  -  P^JrdQ  dr  *  j(P  -  cos  6  d?  dx  (128) 

The  total  pressure  force  acting  on  the  scarfed  conical  extension 
is  obtained  by  integrating  equation  (128)  between  plane  EE'  and  plane 
FF*.  Thus, 

VTf**lFy  029) 
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where 


F 


x 


fyf  ref 

J  ye J  ee  (p  -  Pa)r  de  d r 


{P  “  P,)r  ccs  8  d8  dx 

G 


(130) 

031) 


The  pressure  P  acting  on  the  nozzle  wall  is  a  function  of  the 
local  value  of  x,  but  it  is  not  a  function  o,  the  angle  6  at  a  given 
axial  location.  Consequently,  at  a  given  axial  location,  equations 
(130)  and  (131)  may  be  integrated  over  the  range  <  6  <  it,  where  <jj  is 
a  function  of  the  axial  loca4’’**n.  Thus, 


Fx  =  -  2J  OH»}(p  -  Pa)r  dr  (132) 

x 

x  (P  -  Pjr  sin  ip  dx  (133) 

p  a 


The  factor  2  appearing  in  equations  (132)  and  (133)  accounts  for  the  two 
halves  of  the  nozzle,  since  6  is  integrated  from  \p  to  it. 

In  equations  (132)  end  (133),  the  pressure  P  is  a  function  of  the 
axial  location.  Consequently,  those  two  equations  must  be  integrated 
numerically.  At  any  given  axial  location,  the  integral  of  equation 

(132)  over  the  radial  increment  y^  ^  to  y.  and  the  integral  of  equation 

(133)  over  the  axial  increment  x^^  to  are  given  by 


(AFx)i  =  -2P  y  (tt  -  i>)  (y1  -  yi-1) 
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(134) 


5.  SCARFED  NOZZLE  PERFORMANCE  EVALUATION 


The  total  axial  force  acting  on  the  entire  nozzle  is  obtained  by 
adding  the  force  given  by  equation  (139)  to  the  axial  force  F^  [see 
equation  (109)]  developed  by  the  portion  of  the  nozzle  upstream  of  the 
conical  scarfed  extension.  Thus, 


Fx  =  FIVL  +  FS$  +  Fx,  SCE 


(141) 


The  total  side  force  acting  on  the  entire  nozzle  is  simply  the  side 
force  acting  on  the  scarfed  conical  extension,  given  by  equation  (140). 
Thus, 


Fy  =  Fy.  sce 


(142) 


The  scarfed  nozzle  axial  and  side  specific  impulses  are  given  by 

(143) 

(144) 

The  geometric  relationship  between  the  nozzle  axis  and  the  missile 
axis  is  illustrated  in  Figure  24.  The  missile  coordinate  system  is 
denoted  by  X,Y.  The  nozzle  coordinate  system  x,y  is  assumed  to  lie  in 
the  meridional  plane  through  the  missileaxls,  that  is,  in  the  XY  plane. 
The  forces  acting  on  the  missile  are  thus 


^sp^x 


IF. 


11 

ih 


vIsp^y 


F 

LxL 

m 
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Fx  =  Fy  cos  0  -  F  sin  0 


(H5 


In  the  present  analysis,  the  major  item  of  interest  is  the  axial 
thrust  delivered  by  the  rocket  motor  to  the  missile,  F^,  and  how  that 
thrust  is  related  to  the  axial  thrust  developed  by  t.he  nozzle,  F^,  It 
is  assumed  that  several  rocket  motors  (at  least  two)  are  arranged 
symmetrically  around  the  circumference  of  the  missile,  so  that  all  side 
forces  and  turning  moments  exactly  cancel.  If  that  is  net  the  case, 
the  analysis  may  he  easily  extended  to  determine  the  turning  moments 
associated  with  the  scarfed  nozzle.  That  is  not  done,  however,  in  the 
present  analysis. 

The  final  performance  parameter  of  interest  is  the  ratio  r\  of  the 
axial  thrust  delivered  to  the  missile,  F^,  to  the  axial  thrust  that 
would  be  generated  by  the  nozzle  if  the  scarfed  conical  extension  were 
not  scarfed  and  if  it  were  aligned  along  the  axis  of  the  missile.  That 
value  of  thrust  is  determined  by  calculating  the  axial  force  developed 


by  the  conical  extension  without  scan'  ,ng,  F^ 


,  and  adding  that  value 


to  the  axial  force  acting  on  the  basic  nozzle.  The  axial  force  FrF  is 


determined  by  applying  equation  (108)  from  point  E  to  point  F.  Thus, 


fce-Wp 


(149) 


Note  that  if  the  angle  6^  of  the  scarfed  conical  extension  is  zero,  then 


Fce  is  zero.  The  total  reference  axial  thrust  of  the  unscarfed  scarfed 


nozzle  F^  ^  is  given  by 


F  =  F  +  F  +  F 
rRef  IVL  SS  '  CE 


(150) 


Thus,  the  ratio  n  is  given  by 


(151) 


The  factor  n  may  be  regarded  as  the  efficiency  of  the  scarfed  nozzle. 


From  Figure  24,  it  is  obvious  that  thf*  scarfed  nozzle  exit  flow 


area  lies  on  the  missile  skin.  T*\us,  the  pressun  e  forces  acting  on  that 


area  are  normal  to  the  missile  axis  and  do  not  contribute  to  the  missile 


axial  thrust.  The  missile  axial  thrust  depends  on  t*ie  momentum  flux 


crossing  the  nozzle  exit  area,  which  depends  on  the  nvuzzie  mass  flow 


rate  and  nozzle  exit  velocity.  The  exit  velocity  is  mdependem;  of  the 


pressure  level;  it  depends  on  the  nozzle  geometry,  the  gas  thermodynamic 


model  (i.e.,  y  and  R),  and'  the  stagnation  temperature  T^.  The  (dess  flow 
rate  depends  linearly  on  the  pressure  level,  that  is„  the  -stagnation 


pressure.  Consequently,  the  missile  axial  thrust  awl  nozzle  mass  How 


rate  both  depend  linearly  on  the  stagnation  pressure.  However,  the 
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ratio  of  those  two  quantities,  the  missile  axial  specific  impulse,  is 
independent  of  both  the  stagnation  pressure  and  the  atmospheric  pressure. 
Consequently,  missile  axial  specific  impulse  is  the  most  meaningful  per¬ 
formance  parameter  for  a  scarfed  propulsive  nozzle. 
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SECTION  V 

OVERALL  NUMERICAL  ALGORITHM 


1 .  INTRODUCTION 

The  objective  of  the  present  investigation  was  to  develop  a  pro¬ 
cedure  for  predicting  the  performance  of  a  particular  class  of  scarfed 
propulsive  nozzles.  That  class  of  nozzles  consists  of  a  conventional 
basic  nozzle  which  has  a  double  circular  arc  throat  contour  followed  by 
a  conical,  quadratic,  or  tabular  supersonic  contour.  At  the  end  of  the 
basic  nozzle,  a  scarfed  conical  extension  is  attached.  Figure  25  ill¬ 
ustrates  this  type  of  scarfed  propulsive  nozzle.  A  detailed  discussion 
of  the  nozzle  geometric  model  is  presented  in  Section  II. 

The  flowfield  in  the  transonic  throat  region  of  the  nozzle  is  pre¬ 
dicted  by  a  linearized  flow  analysis.  The  flowfield  in  the  supersonic 
portion  of  the  basic  nozzle  is  predicted  by  the  method  of  characteri¬ 
stics  for  steady  two-dimensional  irrotational  supersonic  flow.  The 
transition  between  the  basic  nozzle  flowfield  and  the  nozzle  extension 
flowfield  occurs  across  an  oblique  shock  wave  that  emanates  from  the 
junction  of  the  basic  nozzle  wall  and  the  nozzle  extension.  The  flow- 
field  in  the  nozzle  extension  is  uredicted  by  the  method  of  character¬ 
istics  for  steady  two-dimensional  rotational  supersonic  flow.  A 
detailed  discussion  of  the  aforementioned  flow  models  is  presented  in 
Section  III. 

The  performance  of  the  scarfed  propulsive  nozzle  is  obtained  by 
integrating  the  flowfield  properties  to  determine  the  overall 
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performance  parameters  (i.e.,  mass  flow  rate,  force  components,  and 
turning  moments).  The  mass  flow  rate  is  obtained  by  integrating  across 
the  transonic  initial-value  line.  The  thrust  of  the  basic  nozzle  is 
determined  by  integrating  the  momentum  flux  and  pressure  forces  across 
the  initial-value  line,  and  adding  to  that  value  the  integral  of  the 
pressure  forces  along  the  nozzle  wall.  The  thrust  of  the  nozzle  exten¬ 
sion  is  obtained  by  integrating  the  pressure  forces  along  the  wall  of 
the  nozzle  extension.  The  total  nozzle  thrust  components  are  determined 
by  adding  the  thrust  components  of  the  basic  nozzle  and  the  nozzle 
extension.  The  thrust  delivered  to  the  missile  is  found  by  calculating 
the  component  of  the  nozzle  thrust  components  in  the  direction  of  the 
missile  axis.  Details  of  the  performance  model  are  discussed  in 
Section  IV. 

The  objective  of  the  present  section  is  to  present  the  logic 
employed  to  combine  all  of  the  individual  analyses  presented  in  Sections 
II  to  IV  into  an  overall  numerical  algorithm  for  predicting  the  perfor¬ 
mance  of  scarfed  propulsive  nozzles. 

A  computer  program  has  been  written  for  implementing  the  analysis 
developed  in  this  investigation.  A  discussion  of  that  program  is  pre¬ 
sented  in  Section  VI.  In  the  present  section  describing  the  overall 
numerical  algorithm,  references  are  made  to  the  particular  program 
routines  that  implement  the  various  aspects  of  the  overall  numerical 
algorithm.  Thus,  program  MAINQO  controls  the  overall  logic  flow  by 
calling,  in  sequence,  program  MAINio  to  read  in  the  input  datat 
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program  MAIN20  to  construct  the  flowfleld  in  the  basic  nozzle*  program 
NAIN30  to  construct  the  flowfte^  in  toe  nozzle  extension,  and  program 
HAIN40  to  calculate  the  performance  of  the  scarfed  nozzle. 
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and  passing  through  the  throat  attachment  point,  point  A,  defines  the 
downstream  extent  of  the  flowfield  along  the  initial  expansion  contour. 

Five  different  modes  of  flowfield  construction  are  considered  In 
the  program.  Those  modes  of  operation  are  described  in  the  following 
paragraphs. 

Mode  1 .  The  first  mode  analyzes  the  flowfield  in  a  nozzle  without 
an  extension  by  constructing  a  network  of  right-running  Hach  lines,  as 
illustrated  in  figure  28.  The  Initial  expansion  flowfield  illustrated 
in  Figure  26  is  first  constructed.  Ripi-r inning  Kacb  lines  are  then 
initiated  from  points  along  the  nozzle  supersonic  contour,  contour  AE. 
These  right-running  Kach  lines  are  continued  until  they  Intersect  the 
nozzle  axis.  Right-running  Nach  line  EE*,  emanating  fi*em  the  nozzle 
exit  lip  point,  point  E,  defines  the  downstream  extent  of  the  nozzle 
flowfield. 

Mode  2.  This  mode  analyzes  the  flowfield  in  a  nozzle  without  an 
extension  by  constructing  a  network  of  left-running  Hach  Vines,  as  illu¬ 
strated  in  Figure  2§.  The  initial  expansion  flowfield  illustrated  in 
Figure  26  is  first  constructed.  Ieft*imn1n§  Mach  lines  are  then  orig¬ 
inated  from  line  AA'.  These  left-running  Mach- lines  are  continued  until 
they  intersect  the  nozzle  wall,  contour  AE.  Left-running  fiach  line  EE* « 
which  passes  through  the  nozzle  exit  lip  point,  point  E,  defines  the 
downstream  extent  of  the  nozzle  flowfield. 

Mode  3,  This  mode  analyzes  the  flowfield  in  a  sozzle  without  an 
extension  in  which  an  ea&edded  right-running  oblique  shock  wave  is 
detected  and  tracked.  The  flowfield  is  constructed  along  a  network  of 
left-running  Hach  lines,  as  illustrated  fe  Figure  30.  Tne  initial 


expansion  flowfield  illustrated  in  Figure  26  is  first  constructed. 
Left-running  Mach  lines  are  then  originated  from  line  AA'.  These  left¬ 
running  Mach  lines  are  continued  until  they  intersect  the  nozzle  wall, 
contour  AE.  The  point  of  initiation,  point  S,  of  the  right-running 
oblique  shock  wave  SS '  must  be  specified  in  terms  of  its  (1,0)  charac¬ 
teristic  coordinates.  Subsequent  left-running  Mach  lines  pass  through 
the  oblique  shock  wave,  which  is  illustrated  as  a  double  line  in  Figure 
30.  The  location  and  properties  of  the  oblique  shock  wave  are  deter¬ 
mined  simultaneously  with  the  construction  of  the  left-running  Mach 
line  network.  Left-running  Mach  line  E'E,  which  passes  through  the 
nozzle  exit  lip  point,  point  E,  defines  the  downstream  extent  of  the 
nozzle  flowfield. 

Mode  4.  This  mode  analyzes  the  flowfield  in  a  nozzle  with  a  scarfed 
nozzle  extension  in  which  an  attached  right-running  oblique  shock  wave 
emanates  from  the  junction  between  the  basic  nozzle  and  the  nozzle 
extension,  point  E.  The  flowfield  is  constructed  along  a  network  of 
left-running  Mach  lines,  as  illustrated  in  Figure  31.  The  initial 
expansion  flowfield  illustrated  in  Figure  26  is  first  constructed.  Left¬ 
running  Mach  lines  are  then  originated  from  line  AA'.  These  left-running 
Mach  lines  are  continued  until  they  intersect  the  nozzle  wall,  contour  AE. 
Left-running  Mach  line  E'E,  which  passes  through  the  nozzle  exit  lip 
point,  point  E,  defines  the  flow  properties  at  point  E.  An  attached 
right-running  oblique  shock  wave  is  originated  from  point  E.  Subsequent 
left-running  Mach  lines  pass  through  the  oblique  shoek  wave.  The  loca¬ 
tion  and  properties  of  the  oblique  shock  wave  are  determined 
simultaneously  with  the  construction  of  the  left-running  Mach  line 


92 


network.  Left-running  Mach  line  F'F,  which  passes  through  point  F,  the 
exit  lip  point  of  the  nozzle  extension,  defines  the  downstream  extent 
of  the  flowfield  in  the  scarfed  nozzle  extension. 

Mode  5.  This  mode  is  identical  to  the  fourth  mode  described  in  the 
previous  paragraph,  except  that  the  initial  expansion  flowfield  illu¬ 
strated  in  Figure  27  is  employed  instead  of  the  initial  expansion 
flowfield  illustrated  in  Figure  26.  The  remainder  of  the  flowfield  is 
constructed  exactly  as  described  in  the  previous  paragraph.  The  result¬ 
ing  flowfield  is  illustrated  in  Figure  32. 


2.  PROBLEM  SPECIFICATIONS 


The  first  step  in  the  procedure  is  the  specification  of  the  prob¬ 
lem  to  be  analyzed.  The  nozzle  geometric  model  must  be  specified  as 
described  in  Section  II  and  illustrated  in  Figure  25.  The  nozzle 
throat  is  specified  by  the  throat  radius  y^,  the  throat  upstream  radius 
of  curvature  p  which  is  employed  in  the  transonic  flow  analysis,  and 

the  throat  downstream  radius  of  curvature  p,  .  which  controls  the 
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initial  supersonic  expansion.  The  x  location  of  the  nozzle  throat  is 
chosen  to  be  0.0. 

The  basic  nozzle  supersonic  contour,  y(x),  attaches  smoothly  to 
the  throat  contour  at  the  attachment  angle  0  .  The  length  of  the  basic 
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nozzle  is  x  .  If  the  supersonic  expansion  contour  is  conical,  the  cone 
angle  must  be  equal  to  0  .  A  conical  nozzle  of  length  x  is  completely 

3  0 

specified  by  any  one  of  the  variables  9a,  y&,  or  e.  If  the  supersonic 

contour  is  quadratic,  0  and  x  must  be  specified.  Then,  one  of  the 
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variables  0  ,  y  ,  or  e  completes  the  specification  or  the  quadratic 
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nozzle.  For  the  tabular  nozzle,  a  set  of  (x,y)  pairs  must  span  the 
region  from  point  A  to  point  E.  The  first  point  in  the  table  must  be 
the  first  point  downstream  of  point  A.  The  last  point  in  the  table 
becomes  point  E.  The  basic  nczzle  geometry  is  described  in  the  com¬ 
puter  program  in  subroutine  BflUNOYE. 

The  conical  extension  is  specified  completely  by  the  values  of  x^ 
and  8^.  The  geometry  of  the  conical  extension  is  described  in  the 
computer  program  in  subroutine  B0UNDYU'. 

The  gas  thermodynamic  model  is  specified  by  the  gas  specific  heat 
ratio  y  and  the  gas  constant  R.  Subroutines  THERM0I  and  TKF.RHOR 


evaluate  the  equations  of  state  for  a  thermally  and  calorically  perfect 
gas  flowing  in  an  irrotational  and  rotational  flowfield,  respectively. 

The  noz2le  operating  conditions  are  specified  by  the  stagnation 
pressure  P^,  the  stagnation  temperature  and  the  ambient  pressure  P&. 

The  parameters  described  in  the  above  paragraphs  completely  specify 
a  specific  nozzle  and  its  operating  conditions.  Numerous  other  para¬ 
meters  must  be  specified  to  enable  a  numerical  solution  of  the  flow- 
field.  Those  parameters  are  described  in  the  following  paragraphs  as 
they  are  encountered. 

The  input  data  are  read  into  the  computer  program  by  subroutine 


INPUT  in  overlay  LINK10. 


3.  INITIAL-VALUE  LINE 

The  first  step  in  the  numerical  solution  is  the  determination  of  a 
supersonic  initial -value  Tine  spanning  the  nozzle  threat  region  fros 
which  the  method  of  characteristics  solution  for  the  supersonic  flow- 
field  can  be  initiated.  Two  options  for  obtaining  a  supersonic  initial- 
value  line  are  contained  in  the  computer  program.  The  first  option  is 
an  internally  generated  initial-value  line.  The  second  option  Is  to 


input  a  tabular  initial-value  line  obtained  from  any  other  source. 


A.  Internally  Generated  Initial-Value  tine 
The  internally  generated  initial-value  line  is  based  os  illegal *$. 
analysis  [Reference  (2)].  That  transonic  flowfield  model  is  described 
in  Sectioa  HI. 2.  The  only  additional  information  required  to  specify 
the  location  of  the  initial -value  line  are  the  number  of  points  desired 
and  the  spacing  of  those  points.  Tht  first  point  is  the  point  on  the 
nozzle  axis  (e,  0.0} »  and  the  last  point  is  the  nozzle  throat  point 


(0.0.  yt). 

The  total  number  of  points  or,  the  initial -value  line  H  HI,  Those 
NI  points  cay  be  spaced  uniformly  along  the  initial-value  line*  or  they 


may  be  spaced  according  to  a  geometric  progression  in  Ay  so  that  the 
ratio  of  the  final  Ay  adjacent  to  the  nozzle  wall  to  the  initial  Ay 
adjacent  to  the  nozzle  asis  is  given  by  the  parameter  DVRATJO,  This 
latter  option,  spacing  according  to  a  geometric  progression,  is  useful 


when  either  throat  radius  of  curvature,  or  p^.  is  small  compared 
to  the  thmt  radius  >v,  so  that  large  gradients  in  flow  properties 


occur  adjacent  to  the  nozzle  wall  in  the  threat  region. 
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Figure  33  illustrates  a  typical  initial-value  line  determined  by 
the  present  algorithm 

B*  Tabular  Initial-Value  tine 

The  superssnk  Initial -value  line  may  be  determined  fros?  ai^y 
external  source  and  read  into  the  computer  prefim  In  tabular  form. 

The  order  of  -the  tabular  data  must  be  consistent  wt  tit  the  order  in 
which  the  data  are  stored  and  used  by  the  remainder  of  the  program, 
the  first  tabular  point  oust  be  the  nozzle  ssis  .points  point  TS 
and  the  last  tabular  point  must  be  the  nozzle  .wall  faints  point  T» 

The  mfeif  of  points  and  their  spacing  is  'arbitrary. 

The  tabular  initial-value  line  H  specified  by  tfet’ Aseter  of  points 
0*1  the  line*  Hi*  the  location  of  the  points  {i.u.®  x  and  y),  and  the 
Mach  §  at  each  point*  Fro©  this  information, 

the  rsp&inder  of  -  f he -f  lps?  at  ..each  point  calculated 

internally  (i.e,F  P*  p?  !>.¥.*  u*  and  v). 
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4.  FLOWFIELD  FROM  THE  INITIAL-VALUE  LINE 


The  flowfield  model  in  the  basic  nozzle  is  described  in  Section 
III. 3.  The  numerical  solution  for  this  flowfield  is  obtained  by  apply¬ 
ing  the  method  of  characteristics  for  steady  two-dimensional  irrotation- 
al  supersonic  flow.  That  procedure  is  implemented  by  constructing 
left-running  and  right-running  Mach  lines,  starting  from  the  initial- 
value  line,  until  they  crisscross  the  entire  supersonic  flowfield.  The 
solution  of  the  flowfield  from  the  initial-value  line  is  described  in 
this  section. 

Figure  34  illustrates  schematically  the  application  of  the  unit 
process  for  an  interior  point  at  the  first  two  points  located  on  the 
•  O'  .'’'  portion  of  the  initial -value  line,  for  determining  the  location 
of  and  the  flow  properties  at  the  downstream  point  of  intersection  of 
the  two  Mach  lines  emanating  from  the  two  initial-value  line  points. 
Figure  35  illustrates  schematically  the  application  of  the  unit  process 
for  an  axis  point,  for  determining  the  location  of  and  the  flow  prop¬ 
erties  at  the  downstream  point  of  intersection  of  the  right-running 
Mach  line  from  the  interior  point  and  the  nozzle  axis. 

The  foregoing  procedure  is  repeated  from  the  next  two  points  on 
the  initial -value  line,  thus  extending  the  corresponding  right-running 
Mach  line  to  the  nozzle  axis.  The  procedure  is  repeated  until  the 
complete  region  determined  by  the  initial -value  line  has  been  determined. 
Figure  36  illustrates  the  resulting  network  of  Mach  lines  emanating  from 
the  Initial-value  line.  The  flowfield  from  the  initial-value  line  is 
now  complete.  To  continue  the  solution  further,  the  wall  boundary 
conditions  must  be  employed. 
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Initial -value  line 


Figure  34.  Application  of  the  unit  process  for  an  interior  point  from 
the  first  two  points  on  the  initial -value  line. 


Figure  35.  Application  of  the  unit  process  for  an  axis  point 
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In  the  computer  program,  the  logic  described  In  this  section  for 
determining  the  flowfield  from  the  initial-value  line  is  implemented  by 
subroutine  MOCIVL  in  overlay  LINK20. 
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5.  FLOWFIELD  FROM  THE  CIRCULAR  ARC  THROAT 

The  flowfield  from  the  circular  arc  throat  is  determined  by  pre- 
specifyingthe  number  and  location  of  points  along  the  circular  arc,  and 
then  applying  the  method  of  characteristics  to  determine  the  flow  prop¬ 
erties  at  those  points.  Two  different  Mach  line  network  construction 
procedures  are  available  for  determining  the  flow  properties  along  the 
initial  expansion  contour. 

The  first  Mach  line  network  construction  procedure  constructs  right-* 
running  Mach  lines  from  the  prespecified  points  along  the  circular  arc 
throat  contour.  These  right-running  Mach  lines  are  propagated  across  the 
flowfield  to  the  nozzle  axis  just  as  the  right-running  Mach  lines  from  the 
initial -value  line  points  were  propagated  to  the  nozzle  axis  (see  Section 
V.4).  This  procedure  is  continued  until  the  right-running  Mach  line  from 
point  A,  the  attachment  point  between  the  circular  arc  throat  and  the 
supersonic  expansion  contour,  has  been  constructed.  This  Mach  line  net¬ 
work,  illustrated  in  Figure  26,  is  used  for  MODE  =  1  to  4,  as  discussed 
in  Section  V.l  and  illustrated  in  Figures  28  to  31. 

The  second  Mach  line  network  construction  procedure  constructs  left¬ 
running  Mach  lines  from  the  right-running  Mach  line  at  the  downstream 
extent  of  the  flowfield  from  the  initial-value  line,  illustrated  in 
Figure  36.  These  left-running  Mach  lines  are  propagated  across  the  flow- 
field  to  the  circular  arc  throat  contour.  This  procedure  is  continued 
until  a  left-running  Mach  passes  downstream  of  the  last  prespecified  wall 
point  on  the  circular  arc  throat  contour,  point  A.  This  Mach  line  network, 
illustrated  in  Figure  27,  is  used  for  MODE  =  5,  as  discussed  in  Section 
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V.1  and  illustrated  In  Figure  32. 

The  total  turning  angle  along  the  circular  arc  throat  Is  9a.  A 
number  of  points  are  prespecified  along  the  circular  arc  contour.  The 
first  point  is  a  point  just  downstream  of  the  throat  point  (i.e.,  point 
T),  and  the  last  point  is  the  attachment  point  between  the  circular  arc 
throat  and  the  supersonic  contour  (I.e. »  point  A).  The  only  additional 
information  required  to  specify  the  location  of  these  points  are  the 
number  of  points  desired  and  the  spacing  of  those  points. 

The  total  number  of  points  on  the  circular  arc  throat  (not  including 
the  throat  point,  point  T)  Is  NT.  These  NT  points  may  be  spaced  along 
the  circular  arc  in  equal  angular  increments #  or  they  may  be  spaced 
according  to  a  geometric  progression  in  A 9  so  that  the  ratio  of  the  final 
AD  adjacent  to  point  A  to  the  initial  AD  adjacent  to  point  T  is  given  by 
the  parameters  DARATIO.  This  latter  option,  spacing  according  to  a  geo¬ 
metric  progression,  is  useful  when  the  throat  downstream  radius  of 
curvature  pt(j  is  small  compared  to  y^,  so  that  large  gradients  in  flow 
properties  occur  adjacent  to  the  nozzle  wall  along  the  circular  arc  con¬ 
tour.  Figure  37  illustrates  a  typical  distribution  of  points  along  the 
throat  downstream  circular  arc  contour. 

Figure  38  illustrates  schematically  the  application  of  the  inverse 
wall  point  unit  process  to  determine  the  first  prespecified  point  on  the 
throat  downstream  circular  arc  contour.  A  left-running  Nach  line  is 
projected  rearward  from  the  prespecified  wall  point  to  intersect  the  pre¬ 
vious  right-running  Hach  line,  as  illustrated  in  Figure  38.  If  the  point 
of  intersection  falls  below  the  second  point,  then  the  point  of  intersec¬ 
tion  is  assumed  to  fall  between  the  second  and  third  points,  as 
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Illustrated  In  Figure  39.  This  procedure  Is  continued  until  the  two 
points  bracketing  the  rearward  projected  left-running  Mach  line  are 
located.  In  most  cases,  these  two  points  are  the  first  two  points  on 
the  previous  right-runni  tig  Mach  line.  The  flow  properties  at  the  point 
of  Intersection  are  then  calculated  by  linear  interpolation.  The  flow 
properties  at  the  prespecified  wall  point  ere  then  found  by  applying  the 
compatibility  relation  along  the  left -running  Mach  line  and  the  solid 
wall  boundary  conditions. 

The  above  procedure  for  determining  the  first  inverse  wall  point  is 
the  same  for  both  of  the  Mach  line  network  construction  procedures  descri¬ 
bed  at  the  beginning  of  this  section.  From  here  on,  however,  the  two 
procedures  are  different.  They  are  described  in  the  next  two  sections. 


Notification  of  the  inverse  wall  point  unit  prx>cess 


6.  RIGHT-RUNNING  MACH  LINE  NETWORK  FROM  THE  CIRCULAR  ARC  THROAT 

Once  the  solution  for  the  Inverse  wall  point  has  been  obtained,  as 
illustrated  in  Figure  38,  a  right-running  Mach  line  is  originated  front 
the  inverse  wal!  point  and  continued  until  it  intersects  the  nozzle  axis, 
as  described  in  Section  V.4  for  the  right-running  Mach  lines  emanating 
from  points  on  the  initial -value  line.  The  above  procedure  is  repeated 
at  each  successive  prespecified  wall  point  nn  the  nozzle  circular  arc 
throat  contour  until  the  region  of  the  flowfield  determined  by  that  por¬ 
tion  of  the  nozzle  contour  has  been  determined.  Figure  40  illustrates 
the  resulting  network  of  Kach  lines  emanating  from  the  nozzle  circular 
arc  throat  contour.  In  the  computer  program,  the  logic  described  in  this 
section  is  implemented  in  subroutine  MOCARCR  in  overlay  LINK20. 


Downstream  extent  of  the  initial  exoansion 
flcwfieid 
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7.  CHARACTERISTIC  COORDINATE  SYSTEM 

By  this  point,  it  should  be  obvious  that  the  numerical  solution  is 
being  obtained  along  the  network  of  left-running  and  right-running  Mach 
lines  that  crisscross  the  flcwfield.  One  of  the  major  advantages  of  the 
numerical  method  of  characteristics  is  that  numerical  information  propa¬ 
gates  along  the  Mach  lines*  which  are  the  actual  paths  of  propagation  of 
physical  information.  Consequently*  the  Kach  line  network  provides  an 
instant  view  of  the  domain  of  dependence  and  range  of  influence  of  each 
and  every  point  in  the  flowfield.  Hence,  it  is  very  useful  to  be  able 
to  identify  specific  individual  Kach  lines  in  a  flcwfield. 

The  identification  of  individual  Mach  lines  is  accomplished  in  the 
numerical  solution  by  eoploying  a  characteristic  coordinate  system,  as 
il lustre,  ted  in  Figure  41.  The  characteristic  coordinate  system  consists 
of  the  right-running  Nach  lines,  called  I  characteristics,  and  the  left¬ 
running  Hach  lines,  called  0  characteristics.  The  characteristic  coor¬ 
dinate  system  is  obviously  not  an  orthogonal  coordinate  system,  but  a 
curvilinear  coordinate  system  composed  of  the  right-running  end  left¬ 
running  Hach  lines. 
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In  the  characteristic  coordinate  system*  the  right-running  Mach 
lines  are  numbered  starting  at  one  for  the  point  on  the  initial-value 
line  at  the  nozzle  axis,  and  increasing  by  one  for  each  succeeding 
right-running  Mach  line.  The  left-running  Mach  lines  are  numbered 
starting  atone  for  the  point  on  the  initial-value  line  at  the  nozzle 
throat,  and  Increasing  by  one  for  each  succeeding  left-running  Mach 
line.  By  following  a  constant  value  of  the  I  or  J  characteristic 
coordinate  through  the  flowfield,  a  unique  right-running  or  left¬ 
running  Mach  line  can  be  traced  from  its  point  of  origin  to  its  final 
point  in  the  flowfield.  The  influence  of  initial -value  line  points  and 
wall  points  on  the  flowfield  can  be  easily  tracked  in  this  manner. 

To  illustrate  the  characteristic  coordinate  system  further,  several 
special  points  are  considered.  The  coordinates  of  the  point  on  the 
initial-value  line  at  the  nozzle  axis  are  (1,  NI),  and  the  coordinates 
of  the  point  at  the  nozzle  throat  are  (MI,  1).  The  coordinates  of  all 
of  the  initial-value  line  points  are  specified  on  Figure  41,  where 
HI  =  5.  Note  that  points  (1,1)  to  (1.4),  (2,1)  to  (2,3),  (3,1),  (3,2), 
and  (4,1)  are  not  defined. 

The  first  prespecified  wall  point  has  an  I  value  of  NI  +  1.  The 
corresponding  0  value  is  that  of  the  left-running  Hach  line  just  above 
the  rearward  projected  left-running  Mach  line  through  the  wall  point. 
Consequently,  the  J  value  of  a  prespecified  wall  point  is  the  same  as 
the  J  value  of  the  previous  wall  point  when  the  rearward  projected 
left-running  Hach  line  intersects  the  previous  right-running  Hach  line 
between  the  first  two  points  on  that  Mach  line.  When  that  intersection 
point  falls  between  the  second  and  third  point,  the  0  value  of  the 
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prespecified  wall  point  Is  one  greater  than  the  J  value  of  the  previous 
prespecified  wall  point.  The  coordinates  of  the  prespecified  wall 
points  are  presented  on  Figure  41.  Clearly,  as  left-running  Mach  lines 
intersect  the  wall  and  dissappear  from  the  Mach  line  network,  the 
characteristic  coordinate  system  reflects  this  situation  by  the  increa¬ 
sing  J  values  of  the  wall  points. 

On  the  nozzle  axis,  right-running  Mach  lines  from  the  initial- 
value  line  or  the  nozzle  wall  intersect  the  axis  and  left-running  Mach 
lines  are  initiated  from  the  intersection  point  and  propagated  into  the 
flowfield.  Thus,  new  left-running  characteristics  having  increasing  J 
values  are  generated  from  the  nozzle  axis.  The  coordinates  of  the  axis 
points  are  specified  on  Figure  41. 

The  characteristic  coordinates  of  all  of  the  interior  points  are 
simply  the  I  and  J  values  of  the  right-running  and  left-running  Mach 
lines  that  intersect  at  each  point.  The  coordinates  of  the  interior 
points  along  the  right-running  Mach  line  from  the  last  prespecified 
wall  point  illustrated  in  Figure  41  are  specified  to  illustrate  the 
coordinates  for  interior  points.  The  coordinates  of  several  other 
Interior  points  are  specified  on  Figure  41  for  illustration. 
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8.  LEFT-RUNNING  MACH  LINE  NETWORK  ADJACENT  TO  THE  CIRCULAR  ARC  THROAT 


The  first  step  in  determining  a  left-running  Mach  line  network  along 
the  initial  expansion  contour  is  the  same  as  the  procedure  illustrated  in 
Figure  38.  After  determining  the  first  inverse  wall  point,  the  rights 
running  Mach  line  emanating  from  that  point  is  extended  into  the 
flowfield  to  find  its  point  of  intersection  with  the  left-running  Mach 
line  just  below  the  interpolated  point  on  the  previous  right-running  Mach 
line.  That  left-running  Mach  line  becomes  the  control  left-running  Mach 
line  in  the  following  procedure. 

So  far,  this  procedure  is  identical  to  the  procedure  presented  in 
Section  V.6.  However,  instead  of  continuing  the  right-running  Mach  line 
to  the  nozzle  axis  as  desribed  in  Section  V.6,  the  next  inverse  wall 
point  is  determined,  as  illustrated  in  Figure  42.  The  right-running  Mach 
line  emanating  from  that  inverse  wall  point  is  then  extended  into  the 
flowfield  to  intersect  the  control  left-running  Mach. 

This  procedure  is  continued  from  succeeding  prespecified  inverse 
wall  points  until  the  point  of  intersection  of  the  rearward  projected 
left-running  Mach  line  from  the  next  inverse  wall  point  and  the  previous 
right-running  Mach  line  lies  below  the  control  left-running  Mach  line,  as 
illustrated  ir  Figure  42.  At  that  point,  that  particular  control  left¬ 
running  Mach  line  is  complete,  and  the  entire  procedure  is  repeated, 
starting  from  the  next  point  on  the  last  right-running  Mach  line  at  the 
downstream  extent  of  the  initial-value  flowfield. 

The  above  procedure  is  repeated  from  succeeding  points  on  the  right¬ 
running  Mach  line  at  the  downstream  extent  of  the  initial-value  flowfield 
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unit  process. 


until  the  last  prespecified  inverse  wall  point  has  been  determined.  The 
final  control  left-running  Mach  line  is  then  extended  forward,  by  the 
direct  wall  point  unit  process,  to  intersect  the  supersonic  turning  con¬ 
tour.  Figure  43  illustrates  the  resulting  network  of  Mach  lines.  This 
Mach  line  network  is  employed  when  analyzing  a  scarfed  nozzle  by  MODE  = 
5.  In  the  computer  program,  the  logic  described  in  this  section  is 
implemented  in  subroutine  MOCARCL  in  overlay  LINK20. 


9.  RIGHT-RUNNING  MACH  LINE  NETWORK  IN  THE  BASIC  NOZZLE 

The  flowfield  from  the  supersonic  contour  (i.e.,  conical,  quadra¬ 
tic,  or  tabular  contour)  is  determined  by  extending  a  left-running  Mach 
line  from  the  second  point  on  the  right-running  Mach  line  emanating 
from  the  previous  wall  point  until  it  intersects  the  supersonic  contour. 
The  flow  properties  at  the  point  of  intersection  are  determined  by  the 
method  of  characteristics.  A  right-running  Mach  line  is  then  emanated 
from  the  wall  points.  The  right-running  fiach  line  is  propagated  across 
the  flowfield  to  the  nozzle  axis  just  as  the  right-running  Mach  lines 
from  the  initial-value  line  (see  Section  V.4)  and  right-running  Mach 
lines  from  the  throat  circular  arc  contour  (see  Section  V.6)  were 
propagated  to  the  nozzle  axis.  This  procedure  is  continued  until  the 
intersection  of  the  left-running  Mach  line  from  the  previous  right¬ 
running  Mach  line  and  the  supersonic  contour  falls  on  the  projection  of 
the  supersonic  contour  beyond  the  end  point  of  the  basic  nozzle,  point 
E.  The  flow  properties  at  point  E  are  determined  by  the  inverse  wall 
point  method  described  in  Section  V.4.  The  final  right-running  Mach 
line  from  point  E  is  propagated  to  the  nozzle  axis,  thus  completing  the 
flowfield  determined  by  the  basic  nozzle  supersonic  contour. 

Figure  44  illustrates  schematically  the  application  of  the  direct 
wall  point  unit  process  to  determine  the  first  point  on  the  supersonic 
contour  just  downstream  of  the  attachment  point,  point  A,  where  the 
throat  circular  arc  contour  joins  the  supersonic  contour.  The  flow 
properties  at  the  direct  wall  point  are  then  found  by  applying  the 
compatibility  relation  along  the  left-running  Mach  line  and  the  solid 
wall  boundary  conditions.  This  procedure  is  implemented  in  the  computer 
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Figure  44.  Application  of  the  direct  wall  point  unit  process. 
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program  by  subroutine  DRWALLI  and  BOUNDYE  in  overlay  LINK20. 

The  complete  flowfield  along  right-running  Mach  lines  emanating  from 
the  supersonic  contour  of  the  basic  nozzle  is  illustrated  schematically 
in  Figure  45.  This  Mach  line  network  is  employed  when  analyzing  a  basic 
nozzle  (MODE  =  1)  or  a  basic  nozzle  with  an  embedded  right-running 
oblique  shock  wave  (MODE  =3).  In  the  computer  program,  the  logic  de¬ 
scribed  in  this  section  is  implemented  in  subroutine  MOCRRC  in  overlay 
LINK20. 

The  procedure  described  above  is  very  straightforward  in  a  contir.* 
uous  flowfield.  However,  when  compression  waves  are  present  they  tend  to 
coalesce,  and  if  enough  coalescence  occurs,  oblique  shock  waves  are 
formed.  The  coalescence  of  Mach  lines  and  the  appearance  of  oblique 
shock  waves  complicates  both  the  flowfield  itself  and  the  logic  procedure 
described  above  for  calculating  the  flowfield.  In  the  present  investi¬ 
gation,  the  assumption  is  made  that  no  strong  oblique  shock  waves  occur 
in  the  basic  nozzle. 

When  right-running  Mach  lines  coalesce  (i.e.,  cross),  the  original 
Mach  line  is  retained  and  the  new  Mach  line  is  terminated  downstream  of 
the  point  of  crossing.  This  procedure  is  illustrated  in  Figure  46. 

Since  the  computational  logic  is  based  on  following  right-running  Mach 
lines  across  the  flowfield  from  the  nozzle  wall  to  the  nozzle  axis,  each 
Mach  line  must  reach  the  nozzle  axis.  The  above  manner  of  handling  cross¬ 
ing  Mach  lines  achieves  that  goal.  The  new  Mach  line  thus  consists  of 
the  new  points  calculated  before  the  point  of  crossing  and  the  old  points 
downstream  of  the  point  of  crossing.  Thus,  the  new  Mach  line  completely 
spans  the  nozzle  flowfield  from  the  nozzle  wall  to  the  nozzle  axis.  In 


the  computer  program  output,  the  points  downstream  of  the  point  of 
crossing  are  written  out  as  the  solution  for  both  right-running  Mach 
lines. 

When  left-running  Mach  lines  cross,  the  new  Mach  line  is  terminated. 
Since  the  computational  logic  is  based  on  following  right-running  Mach 
lines  across  the  flowfield  instead  of  following  left-running  Mach  lines, 
no  further  action  is  necessary  when  left-running  Mach  lines  cross.  This 
procedure  is  illustrated  in  Figure  47.  Essentially,  the  left-running 
Mach  line  that  crossed  its  upstream  neighbor  has  been  assumed  to  cease 
to  exist.  This  would  create  errors  in  the  (I,J)  characteristic  coordi¬ 
nate  system  if  it  were  not  accounted  for.  In  the  present  computer 
program,  when  a  left-running  Mach  line  crosses  its  upstream  neighbor, 
the  solution  is  aborted  and  the  x  location  of  the  solution  point  is 
specified  by  -1.0.  The  aborted  point  is  given  a  location  in  the  solution 
storage  arrays  and  the  (I,J)  characteristic  coordinate  system.  That 
left-running  Hach  line  then  acts  like  a  phantom  Mach  line  until  it  passes 
out  of  the  flowfield  at  the  nozzle  wall,  as  illustrated  in  Figure  4S. 

In  any  particular  flowfield,  several  phantom  left-running  Mach  lines  may 
exist  simultaneously,  adjacent  to  each  other  or  separated  by  real  left¬ 
running  Kach  lines.  In  the  computer  program  output,  these  phantom 
points  are  written  out  with  their  I,J  characteristic  coordinates  and  a 
blank  line  of  data. 

The  occurence  of  crossed  right-running  Mach  lines  and  phantom 
left-running  Mach  lines  considerably  complicates  the  program  logic. 
However,  the  procedure  described  above  works  very  well  and  gives  very 
accurate  results  as  long  as  the  Mach  line  crossings  are  not  too  numerous. 
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10.  LEFT-RUNNING  MACH  LINE  NETWORK  IN  THE  BASIC  NOZZlE 


A  left-running  Mach  line  network  can  be  constructed  in  the  basic 
nozzle  starting  from  either  line  AA‘,  illustrated  in  Figure  40,  or  line 
TT" ,  illustrated  in  Figure  43.  Modes  2,  3,  and  4,  described  in  Section 
V.2,  construct  left-running  Mach  line  networks  starting  from  line  AA\ 
Mode  5  constructs  a  left- running  Mach  line  network  starting  from  line 
TT".  Those  two  network  construction  procedures  are  described  in  this 
section. 

When  the  left-running  Mach  line  network  construction  procedure 
illustrated  in  Figure  40  is  employed,  the  flowfield  is  constructed  along 
left-running  Mach  lines  emanating  from  the  last  right-running  Mach  line 
at  the  downstream  extent  of  the  initial  expansion  flowfield.  That  right¬ 
running  Mach  line  is  called  the  control  right-running  Mach  line.  Start¬ 
ing  from  the  flowfield  illustrated  in  Figure  40,  a  left-running  Mach  line 
is  initiated  from  the  second  point  on  the  control  right-running  Mach  line 
That  left- running  Mach  line  is  continued  upward  until  it  intersects  the 
nozzle  wall,  as  illustrated  in  Figure  49(a). 

This  procedure  is  repeated  from  successive  points  along  the  control 
right-running  Mach  line  until  one  of  two  events  occurs:  (1)  a  left¬ 
running  Mach  line  from  the  last  point  {i.e.,  the  point  on  the  axis)  on 
the  control  right-running  Mach  line  has  been  generated,  or  (2)  the 
intersection  of  the  left-running  Mach  line  and  the  supersonic  contour 
falls  on  the  projection  of  the  supersonic  contour  beyond  the  end  of  the 
basic  nozzle,  point  E.  In  the  first  case0  an  axis  point  is  determined  as 
illustrated  in  Figure  35,  and  the  next  left-running  Mach  line  emanates 
from  that  axis  point.  In  the  second  case,  the  flow  properties  at  point  E 
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are  determined  by  the  inverse  wall  point  method  described  in  Section  V.4. 
This  completes  the  flowfield  determined  by  the  basic  nozzle  contour. 

When  the  left-running  Mach  line  network  construction  procedure 
illustrated  in  Figure  43  is  employed,  the  flowfield  is  constructed  along 
left-running  Mach  lines  emanating  from  the  last  right-running  Mach  line 
at  the  downstream  extent  of  the  flowfield  from  the  initial-value  line, 
line  TT'.  That  right-running  Mach  line  is  called  the  control  right¬ 
running  Mach  line.  Enough  left-running  Mach  lines  have  already  been 
initiated  from  the  control  right-running  Mach  line  to  determine  the  flow 
properties  along  the  initial  expansion  contour,  contour  TA.  That  proce¬ 
dure  is  discussed  in  Section  V.8  and  illustrated  in  Figure  42. 

The  procedure  described  in  Section  V.8  is  continued  from  successive 
points  along  the  control  right-running  Mach  line.  Each  successive  left¬ 
running  Mach  line  is  continued  upward  until  it  intersects  the  nozzle 
wall,  as  illustrated  in  Figure  49(b).  The  procedure  from  here  on  is 
identical  to  that  described  above  where  line  AA'  is  the  control  right¬ 
running  Mach  line. 

The  flowfield  determined  by  either  of  the  two  aforementioned  Mach 
line  network  construction  procedures  is  illustrated  schematically  in 
Figure  50.  This  Mach  line  network  construction  procedure  is  employed 
when  analyzing  a  basic  nozzle  (MODE  =  2),  a  basic  nozzle  with  an  embedded 
right-running  oblique  shock  wave  (MODE  =  3),  or  a  scarfed  nozzle  with  an 
attached  right-running  oblique  shock  wave  (MODE  =  4  or  5),  In  the  com¬ 
puter  program,  the  logic  described  in  this  section  is  implemented  in 
subroutine  NOCLRCI  in  overlay  UNK20. 

The  procedure  described  above  is  very  straightforward  in  a 
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continuous  flowfield.  However,  when  compression  waves  are  present  they 
tend  to  coalesce,  and  if  enough  coalescence  occurs,  oblique  shock  waves 
are  formed.  The  coalescence  of  Mach  lines  and  the  appearance  of  oblique 
shock  waves  complicates  both  the  flowfield  itself  and  the  logic  procedure 
described  above  for  calculating  the  flowfield.  In  the  present  investi¬ 
gation,  the  assumption  is  made  that  no  strong  oblique  shock  waves  occur 
in  the  basic  nozzle. 

When  left-running  Mach  lines  coalesce  (i.e.,  cross),  the  original 
Mach  line  is  retained  and  the  new  Mach  line  is  terminated  downstream  of 
the  point  of  crossing.  This  procedure  is  illustrated  in  Figure  51. 

Since  the  computational  logic  is  based  on  following  left-running  Mach 
lines  across  the  flowfield  from  the  downstream  extent  of  the  initial 
value  flowfield  or  the  nozzle  axis  to  the  nozzle  wall,  each  left-running 
Mach  line  must  reach  the  nozzle  wall.  The  above  manner  of  handling 
crossing  Mach  lines  achieves  that  goal.  The  new  Mach  line  thus  consists 
of  the  new  points  calculated  before  the  point  of  crossing  and  the  old 
points  after  the  point  of  crossing.  Thus,  the  new  Mach  line  completely 
spans  the  nozzle  flowfield  from  the  downstream  extent  of  the  initial 
value  flowfield  or  the  nozzle  axis  to  the  nozzle  wall.  In  the  computer 
program  output,  the  points  after  the  point  of  crossing  are  written  out  as 
the  solution  for  both  left-running  Mach  lines. 

When  right-running  Mach  lines  cross,  the  new  Mach  line  is  terminated. 
Since  the  computational  logic  is  based  on  following  left-running  Mach 
lines  across  the  flowfield  instead  of  following  right-running  Mach  lines, 
no  further  action  is  necessary  when  right-running  Mach  lines  cross.  This 
procedure  is  illustrated  in  Figure  52.  Essentially,  the  right-running 
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Mach  line  that  crossed  its  upstream  neighbor  has  been  assumed  to  cease 
to  exist.  This  would  create  errors  in  the  (1,0)  characteristic  coordi¬ 
nate  system  if  it  were  not  accounted  for.  In  the  present  computer  pro¬ 
gram,  when  a  right-running  Mach  line  crosses  its  upstream  neighbor,  the 
solution  is  aborted  and  the  x  location  of  the  solution  point  is  specified 
by  -1.0.  The  aborted  point  is  given  a  location  in  the  solution  storage 
arrays  and  the  (I , J)  characteristic  coordinate  system.  That  right-running 
Mach  line  then  acts  like  a  phantom  Mach  line  until  it  passes  out  of  the 
flowfield  at  the  nozzle  axis,  as  illustrated  in  Figure  53.  In  any  part¬ 
icular  flowfield,  several  phantom  left-running  Mach  lines  may  exist 
simultaneously,  adjacent  to  each  other  or  separated  by  real  left-running 
Mach  lines.  In  the  computer  program  output,  these  phantom  points  are 
written  out  with  their  I , J  characteristic  coordinates  and  a  blank  line 
of  data. 

The  occurence  of  crossed  left-running  Mach  lines  and  phantom  right¬ 
running  Mach  lines  considerably  complicates  the  program  logic.  However, 
the  procedure  described  above  works  very  well  and  gives  very  accurate 
results  as  long  as  the  Mach  line  crossings  are  not  too  numerous. 


11.  OBLIQUE  SHOCK  HAVE 

An  oblique  shock  emanates  from  the  junction  point  of  the  basic 
nozzle  and  the  nozzle  extension  point  E  in  Figure  25.  This  oblique 
shock  wave  propagates  upstream  into  the  flowfield  emanating  from  the 
supersonic  contour  of  the  basic  nozzle.  The  flowfield  model  for  this 
oblique  shock  wave  is  presented  in  Section  III. 3.  The  implementation  of 
that  flowfield  model  into  the  overall  numerical  algorithm  is  discussed 
in  the  present  section. 

The  first  step  in  the  determination  of  the  oblique  shock  wave  is 
the  determination  of  the  wave  angle  e  and  the  flow  properties  on  the 
downstream  side  of  the  oblique  shock  wave  at  point  E  itself.  This  is 
accomplished  by  assuming  a  value  for  the  wave  angle  e,  solving  for  the 
flow  turning  angle  6  from  equation  (77) ,  and  comparing  the  resulting 
flow  angle  (i.e.,  8Q  -  d)  with  the  angle  of  the  conical  extension  (I.e,, 
Gf).  The  wave  angle  e  is  varied  iteratively  until  (Qe  -  6)  converges  to 
within  a  prespecified  convergence  tolerance.  The  downstream  Mach 
number  Mg  is  then  determined  from  equations  (74)  to  (76),  and  the 
remaining  flow  properties  are  determined  from  equations  (78)  to  (82). 
This  completes  the  solution  for  the  oblique  shock  wave  at  point  E. 

One  small  modification  in  the  characteristic  network  is  made 
before  the  oblique  shock  wave  solution  is  initiated.  The  interpolated 
point  at  the  intersection  of  the  reward  projection  of  the  left-running 
Kach  line  through  point  E  is  substituted  into  the  solution  arrays  in 
place  of  the  point  just  below  it,  as  illustrated  in  Figure  $4, 
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The  next  step  In  determining  the  flowfield  downstream  of  the  oblique 
shock  wave  Is  to  initiate  a  left-running  Hach  line  from  the  next  point  on 
the  right-running  Mach  line  at  the  downstream  extent  of  the  upstream 
flowfleld  or  from  the  nozzle  axis  If  the  upstream  flowfield  has  reached 
the  axis.  That  left-running  Mach  line  Is  then  propagated  toward  the 
scarfed  nozzle  extension  until  It  reaches  the  oblique  shock  wave.  It 
penetrates  the  oblique  shock  wave  and  continues  until  it  intersects  the 
scarfed  nozzle  exte^ion,  as  Illustrated  in  Figure  55. 

The  numerical  procedure  for  obtaining  the  solution  at  the  intersec¬ 
tion  of  the  left-running  Hach  line  and  the  oblique  shock  wave  Is 
presented  In  Section  III.4.C.  The  application  of  that  model  to  the  first 
point  on  the  oblique  shock  wave  downstream  of  point  E  is  illustrated  in 
Figure  56.  The  circled  points  are  temporary  points  Inserted  on  rearward 
projected  right-running  Mach  lines.  The  shock  wave  point  is  located  at 
the  downstream  intersection  of  the  oblique  shock  wave  from  point  E  and 
the  left-running  Mach  line  from  the  upstream  flowfield.  The  oblique 
shock  wave  Is  Illustrated  schematically  as  a  double  line  to  Indicate  that 
there  are  two  sets  of  flow  properties  at  each  point  on  the  shock  wave; 
the  upstream  properties  determined  by  application  of  the  upstream  inter- 
io  point  unit  process  and  the  downstream  properties  determined  by  the 
oblique  shock  wave  equations.  The  first  wall  point  downstream  of  point 
E  is  also  determined  during  the  iteriative  procedure  for  the  shock  wave 
angle  at  the  shock  wave  point. 

After  the  first  point  on  the  oblique  shock  wave  downstream  of  point 
E  has  been  determined,  the  next  point  on  the  shock  wave  is  determined  by 
the  procedure  discussed  in  Section  III.4.C.  That  procedure  is  illustrated 
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Figure  56.  First  point  on  the  oblique  shock  wave  downstream  of  point  E. 
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In  Figure  57.  This  procedure  Is  repeated  from  subsequent  points  at  the 
downstream  extent  of  the  upstream  flowfield  until  the  entire  flowfield 
In  the  scarfed  nozzle  extension  has  been  determined. 

Since  the  strength  of  the  oblique  shock  wave  Is  greater  than  the 
strength  of  a  Mach  line,  the  oblique  shock  wave  propagates  at  a  steeper 
angle  than  the  upstream  right-running  Mach  line.  Eventually  the  shock 
wave  will  overtake  the  upstream  Mach  line. 

When  the  oblique  shock  wave  overtakes  the  upstream  right-running 
Mach  line,  that  crossed  right-running  Mach  line  is  simply  terminated. 

The  solution  for  the  oblique  shock  wave  is  then  obtained  from  the  right¬ 
running  Mach  line  preceeding  the  crossed  right-running  Mach  line.  The 
crossed  right-running  Mach  line  is  then  treated  as  a  phantom  right¬ 
running  Mach  line  In  the  characteristic  coordinate  system,  as  discussed 
in  Section  V.10  and  illustrated  in  Figure  53.  This  procedure  of  termi¬ 
nating  the  crossed  right-running  Mach  line  and  dropping  back  to  the 
previous  right-running  Mach  line  is  repeated  until  the  shock  wave  point 
is  successfully  located.  This  procedure  is  illustrated  in  Figure  58. 

The  right-running  Mach  line  downstream  of  the  oblique  shock  wave  is 
steeper  than  the  shock  wave  because  of  the  compression  and  deceleration 
of  the  flow  caused  by  the  shock  wave.  Eventually  the  downstream  Mach 
line  will  overtake  t m  shock  wave.  When  the  downstream  right-running 
Mach  line  overtakes  the  oblique  shock  wave,  that  downstream  right-running 
Mach  is  terminated.  It  is  then  treated  as  a  phantom  right-running  Mach 
line  in  the  characteristic  coordinate  system,  as  discussed  in  Section  V.10 
and  illustrated  in  Figure  53.  The  next  solution  point  Is  than  determined 
at  the  intersection  of  the  control  left-running  Hacti  line  down  stream  of 
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Figure  57,  Application  of  the  oblique  shock  wave  point  unit  process. 


12.  FLOWFIELD  IN  THE  SCARFED  NOZZLE  EXTENSION 


The  flowfleld  in  the  scarfed  nozzle  extension  consists  of  the  flow- 
field  downstream  of  the  oblique  shock  wave  emanating  from  point  E,  the 
junction  point  of  the  basic  nozzle  and  the  nozzle  extension.  Due  to  the 
presence  of  an  entropy  gradient  caused  by  the  curved  oblique  shock  wave, 
that  flowfleld  Is  rotational.  The  flowfleld  model  for  a  rotational  flow 
is  presented  in  Section  III. 5.  The  implementation  of  that  procedure  to 
determine  the  flowfleld  in  the  nozzle  extension  is  described  in  this 
section. 

The  flowfield  in  the  scarfed  nozzle  extension  is  determined  by 
initiating  left-running  Mach  lines  from  the  upstream  flowfield,  as 
illustrated  in  Figure  55,  and  propagating  those  Mach  lines  to  the  scarfed 
nozzle  extension.  This  procedure  is  repeated  for  succeeding  left-running 
Mach  lines  until  one  of  the  three  following  situations  arises: 

1.  The  oblique  shock  wave  overtakes  the  right-running  Mach  line  at  the 
downstream  extent  of  the  upstream  flowfield  before  the  end  of  the 
scarfed  nozzle  extension  has  been  reached, 

2.  The  oblique  shock  wave  intersects  the  nozzle  axis  before  the  end  of 
the  scarfed  nozzle  extension  has  been  reached. 

3.  The  left-running  Mach  line  intersects  the  extension  of  the  scarfed 
nozzle  extension  beyond  point  F.  In  this  case,  the  indirect  wall 
point  unit  process  is  applied  at  point  F,  as  illustrated  in  Figure 
60,  and  the  scarfed  nozzle  flowfield  is  complete.  This  is  the  normal 
mode  of  operation. 

Situations  (1)  and  (2)  are  illustrated  in  Figures  61  and  62 
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respectively.  When  the  oblique  shock  wave  reaches  the  nozzle  axis,  the 
calculations  are  terminated.  In  an  actual  flowfield,  the  oblique  shock 
wave  generally  transitions  to  a  Mach  disc  and  an  embedded  subsonic  zone 
In  the  neighborhood  of  the  nozzle  axis.  That  phenomena  cannot  be  pre¬ 
dicted  by  the  method  of  characteristics.  In  fact,  the  prediction  of  the 
Mach  disc  and  the  subsonic  flow  region  Is  extremely  difficult,  and 
considerably  beyond  the  scope  of  the  present  investigation.  Consequently, 
in  the  present  Investigation,  the  oblique  shock  wave  is  simply  terminated 
at  the  last  interior  point  above  the  nozzle  axis. 

Practically  speaking,  this  situation  does  not  arise  frequently. 

When  It  does  occur,  the  last  point  actually  calculated  is  generally  quite 
close  to  point  F.  In  this  situation,  the  approximation  is  made  that 
a  uniform  flow  exists  from  the  last  predicted  wall  point  to  point  F. 

Since  point  F  is  at  the  last  remnant  of  the  scarfed  nozzle  extension  and 
the  pressure  level  in  this  portion  of  the  flowfield  is  very  nearly  atmos¬ 
pheric  in  many  cases,  the  aforementioned  assumption  is  generally  quite 
good.  When  the  last  calculated  point  is  far  upstream  of  point  F  and/or 
the  internal  pressure  level  is  greatly  different  from  atmospheric  pres¬ 
sure,  the  aforementioned  assumption  may  Introduce  a  small  but  significant 
error  into  the  nozzle  side  force,  and  an  even  smaller  and  probably  insig¬ 
nificant  error  into  the  nozzle  axial  force.  In  that  case,  some  judgement 
as  to  the  accuracy  of  the  solution  must  be  made  based  on  experience. 

The  flowfield  described  in  Situation  (1)  above  can  sometimes  be 
continued  further  by  using  the  MODE  «  5  option.  That  option  esiploys  a 
left-running  Mach  line  network  construction  procedure  for  the  flowfield 
along  the  initial  expansion  contour.  Consequently,  the  oblique  shock 


wave  can  propagate  all  the  way  to  the  nozzle  axis  before  the  calculations 
must  be  terminated. 

The  ijsplssentation  of  the  logic  described  in  this  section  for 
determining  the  flowfield  in  the  scarfed  nozzle  extension  is  contained  in 
subroutine  HQCLRCR  in  overlay  LINK30. 


13.  SCARFED  NOZZLE  PERFORMANCE 

The  scarfed  nozzle  performance  model  is  presented  in  Section  IV. 

The  application  of  that  model  is  discussed  in  this  section.  The  imple¬ 
mentation  of  this  model  in  the  computer  program  is  contained  in  overlay 
LINK40. 

The  entire  flowfield  in  the  scarfed  nozzle  can  be  calculated  by  the 
procedures  presented  in  the  preceding  sections.  The  performance  of  the 
basic  nozzle  is  specified  in  terms  of  the  nozzle  mass  flow  rate  A  and  the 
axial  thrust  F^. 

Two  steps  are  required  to  evaluate  the  performance  of  the  scarfed 
nozzle  extension.  The  first  step  is  the  calculation  of  the  angle  tp, 
given  by  equation  (121).  That  calculation  is  performed  in  subroutine 
GECH  in  overlay  LINK30.  The  second  step  is  the  calculation  of  the 
scarfed  nozzle  thrust  components  SCE  and  Fy#  SCE>  given  by  equations 
(139)  and  (140).  That  calculation  is  performed  in  subroutine  FORNOH  in 
overlay  UNK30, 

The  final  step  in  the  overall  numerical  algorithm  is  the  calculation 
of  the  performance  of  the  entire  nozzle.  That  performance  Is  specified 
by  Fx,  Fy,  (Isp)^r  tISp)y*  ancJ  tha  effective  scarfing  angle  Those 

parameters  are  defined  in  Section  IV. 4.  Figure  63  illustrates  the 
relationship  between  the  thrust  components  of  the  scarfed  nozzle  and  the 
thrust  components  transmitted  to  the  missile.  The  calculation  of  the 
overall  performance  parameters  is  accomplished  in  subroutine  F0RHQM  in 
overlay  UNK30. 
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SECTION  VI 


COMPUTER  PROGRAM 


1.  INTRODUCTION 

A  computer  program  has  been  developed  to  calculate  the  performance 
of  a  scarfed  propulsive  nozzle.  A  brief  description  of  that  program  is 
presented  in  this  section.  The  overlay  structure  of  the  program  is 
described,  and  a  brief  description  of  the  function  of  each  program  and 
subroutine  is  presented. 

2.  PROGRAM  OVERLAY  STRUCTURE 

The  overall  program  consists  of  52  routines  (5  program  routines  and 
47  subroutines).  The  program  is  organized  into  an  overlay  scheme  to 
minimize  the  amount  of  computer  memory  required.  The  overlay  structure 
is  presented  in  Figure  64.  Two  overlay  levels  are  used:  the  resident 
overlay  level  and  the  primary  overlay  level. 

OVERLAY  (0,0)  is  the  resident  overlay  which  controls  the  overall 
execution  of  the  program.  Several  commonly  used  subroutines  are  included 
in  this  overlay.  OVERLAY  (1,0)  perofrms  data  input,  parameter  initiali¬ 
zation,  and  generation  of  the  initial-vaiue  line  across  the  nozzle 
throat.  OVERLAY  (2,0)  performs  the  evaluation  of  the  flowfield  in  the 
basic  nozzle.  OVERLAY  (3,0)  constructs  the  oblique  shock  wave  emanating 
from  the  junction  between  the  basic  nozzle  and  the  nozzle  extension  and 
evaluates  the  flowfield  in  the  nozzle  extension.  OVERLAY  (4,0)  calcu¬ 
lates  the  tnrust  components  in  the  scarfed  extension,  and  determines  the 
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performance  of  the  overall  scarfed  propulsive  nozzle. 

In  the  following  sections,  a  brief  description  is  given  of  the 
function  of  each  subroutine  in  the  computer  program.  This  information 
supplements  the  information  available  within  the  program  in  the  form  of 
comment  statements. 


3.  OVERLAY  (0,0) 


MAINOO.  This  program  routine  Is  the  main  control  routine  In 
OVERLAY  (0,0),  the  resident  overlay.  MAINOO  first  calls  OVERLAY  (1,0) 
for  data  input,  parameter  initialization,  and  generation  of  the  initial- 
value  line.  MAINOO  then  calls  OVERLAY  (2,0)  to  calculate  the  flowfield 
in  the  basic  nozzle,  OVERLAY  (3,0)  to  calculate  the  oblique  shock  wave 
and  the  flowfield  in  the  nozzle  extension,  and  OVERLAY  (4,0)  to  calcu¬ 
late  the  nozzle  performance. 

THERMQI.  This  subroutine  calculates  the  temperature  T,  pressure  P, 
density  p,  speed  of  sound  a,  and  Mach  number  M,  corresponding  to  a 
specified  value  of  the  velocity  magnitude  V,  for  the  irrotational  flow 
of  a  thermally  and  calorically  perfect  gas. 

THERMOR.  This  subroutine  calculates  the  temperature  T,  speed  of 
sound  a,  and  Mach  number  M,  corresponding  to  specified  values  of  the 
pressure  P,  density  p,  and  velocity  magnitude  V,  for  the  rotational  flow 
of  a  thermally  and  calorics! ly  perfect  gas. 

OUTPUT.  This  subroutine  writes  out  the  flow  properties  I,  0,  x,  y, 
u,  v,  M,  V,  0,  P,  p,  T,  H,  m,  and  F  at  a  point  in  the  flowfield. 

THRUST.  This  subroutine  calculates  the  increment  in  thrust  along 
the  wall  due  to  the  pressure  forees  acting  on  the  wall,  the  total  nozzle 
thrust,  and  the  nozzle  specific  impulse,  as  described  in  Section  IV. 2. 

The  isentropic  one-dimensional  values  are  calculated  for  comparison,  and 
the  nozzle  performance  parameters  F,  F^_q»  hp  !$p»  and  are 

written  out. 

DW00TR.  This  subroutine  calculates  the  Increments  in  mass  flow  rate 
and  thrust  between  two  adjacent  points  on  a  right-running  Mach  line 
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using  the  trapezoidal  rule. 

DWDOTL.  This  subroutine  calculates  the  increment  in  mass  flow  rate 
and  thrust  between  two  adjacent  points  on  a  left-running  Mach  line  using 
the  trapezoidal  rule. 

BQUNOVW.  This  subroutine  contains  the  description  of  the  super¬ 
sonic  contour  of  the  basic  nozzle  (see  Section  II. 2).  The  point  of 
intersection  of  a  left-running  Mach  line  and  the  supersonic  contour  is 
iocated  in  this  subroutine. 

BOUNDYE.  This  subroutine  contains  the  description  of  the  nozzle 
extension  contour  (see  Section  II. 3).  The  point  of  intersection  of  a 
left-running  Mach  line  and  the  ncz2le  extension  is  located  in  this  sub¬ 
routine. 

INTERI .  This  subroutine  implements  the  method  of  characteristics 
solution  for  an  interior  point  in  an  irrotational  flowrield,  as  illu¬ 
strated  in  Figure  27. 

AXIS!.  This  subroutine  implements  the  method  of  characteristics 
solution  for  an  axis  point  in  an  irrotational  flowrield,  as  illustrated 
in  Figure  ZB. 

PDATA.  This  subroutine  saves  the  coordinates  of  the  points  along 
right-running  ?4ach  lines  for  subsequent  plotting  of  the  right-running 
Mach  line  network. 

PRINT.  This  subroutine  prints  out  the  data  saved  by  subroutine 
PDATA. 

PLOTWAl ,  This  subroutine  initializes  the  plotting  subroutines  and 
plots  the  wall  contour. 

ROTRRC.  This  subroutine  plots  the  right-running  Mach  line  network, 

using  the  data  saved  by  subroutine  PDATA. 
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4.  OVERLAY  (1,0) 


HAIN10.  This  program  is  the  main  control  routine  in  OVERLAY  (1,0). 
Subroutine  INPUT  is  called  to  read  in  and  write  out  the  input  data.  Sub¬ 
routine  IVLINE  is  called  to  generate  the  supersonic  initial-value  line. 

INPUT.  This  subroutine  reads  in  the  input  data  and  problem  specifi¬ 
cations  through  namelists  DATA  and  WALL.  All  of  the  input  variables  are 
defined  by  comment  cards.  Default  values  are  assigned  to  most  of  the 
input  variables  before  the  namelists  are  read  in.  Several  parameters 
are  initialized,  and  the  units  conversion  factors  are  specified.  The 
program  description,  the  job  title,  and  the  problem  specifications  are 
written  out.  If  requested,  subroutine  COMPRES  is  called  to  compress  the 
tabular  wall  contour. 

IVLINE.  This  subroutine  selects  the  type  of  supersonic  initial-value 
line  to  be  employed. 

KLIEG.  This  subroutine  internally  generates  a  supersonic  initial-value 
line  by  the  perturbation  analysis  developed  by  Kliegel  and  Levine  (2).  The 
location  (i.e.,  x  and  y  coordinates)  of  the  points  on  the  initial -value 
line  are  specified,  and  the  velocity  components  (i.e.,  u  and  v)  are  calcu¬ 
lated. 

IVLTAB.  This  subroutine  reads  in  the  coordinates  (i.e.,  x  and  y),  the 
Mach  number  M,  and  the  flow  angle  e  along  a  tabular  initial-value  line.  The 
velocity  components  (i.e.,  u  and  v)  are  internally  calculated  at  each  point. 

IVLPERF.  This  subroutine  calculates  all  of  the  performance  param¬ 
eters  for  the  initial-value  line  from  the  specified  stagnation  pressure 
Pt  and  temperature  Tt,  the  velocity  components  u  and  v,  and  the  coordinates 
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x  and  y  of  the  points  on  the  Initial-value  line.  The  calculated  parameters 
include  the  Mach  number  M,  the  velocity  magnitude  V»  the  flow  angle  e(  the 
static  pressure  P,  density  pt  and  temperature  T,  the  stagnation  pressure 
Pj.  and  temperature  Tt»  the  mass  flow  rate  m  flowing  between  each  point  and 
the  nozzle  axis,  and  the  thrust  associated  with  that  mass  flow  rate.  The 
nozzle  overall  performance  parameters  are  also  written  out.  These  include 
the  mass  flow  rate  ra,  the  discharge  coefficient  CD»  the  thrust  F,  the 
thrust  efficiency  np  the  specific  impulse  Isp,  and  the  specific  impulse 
efficiency  np 

COMPRES.  This  subroutine  compresses  the  axial  coordinates  of  the 
tabular  wall  data. 

ATTACH.  This  subroutine  attaches  the  compressed  tabular  contour 
smoothly  to  the  circular  arc  initial -expansion  contour. 
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5.  OVERLAY  (2,0) 


MAIN20.  This  program  is  the  main  control  routine  in  OVERLAY  (2,0). 
Subroutines  MOCIVL,  MOCARCR ,  MOCARCL ,  MOCRRC,  and  MOCLRCI  are  called  to 
calculate  the  flowfield  emanating  from  the  initial -value  line,  the 
nozzle  throat  downstream  circular  arc  contour,  and  the  supersonic  expan¬ 
sion  contour. 

MOCIVL.  This  subroutine  is  the  control  subroutine  for  calculating 
the  flowfield  emanating  from  the  supersonic  initial-value  line.  Sub¬ 
routines  INTERI  and  AXISI  are  called,  as  described  in  Section  V.4,  to 
calculate  the  flowfield.  Subroutine  MOVE!  is  called  to  transfer  data 
between  computational  points  and  storage  arrays,  as  required.  Subroutine 
DWDOTR  is  called  to  integrate  the  mass  flow  rate  and  thrust  along  the 
right-running  Mach  lines  as  they  are  constructed.  Subroutine  OUTPUT  is 
called  to  write  out  the  flow  properties  at  each  solution  point. 

MOCARCR.  This  subroutine  is  the  control  subroutine  for  calculating 
the  flowfield  along  right-running  Mach  lines  emanating  from  the  throat 
downstream  circular  arc  contour  of  the  basic  nozzle  as  described  in 
Section  V.6.  Subroutines  INUALll,  INTERI,  and  AXISI  are  called  to  cal¬ 
culate  the  flowfield.  Subroutine  M0VE1  is  called  to  transfer  data 
between  computational  points  and  solution  arrays,  as  required.  Subrou¬ 
tine  THRUST  is  called  to  integrate  the  thrust  oenerated  along  the  nozzle 
wall,  as  described  in  Section  IV. 2.  Subroutine  DWDOTR  is  called  to 
integrate  the  mass  flow  rate  and  thrust  along  the  right-running  Mach 
lines  as  they  are  constructed.  Subroutine  OUTPUT  is  called  to  write  out 
the  flow  properties  at  each  solution  point. 
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MOCARCL.  This  subroutine  is  the  control  subroutine  for  calculating 
the  flowfield  along  left-running  Mach  lines  along  the  initial  expansion 
contour  of  the  basic  nozzle,  described  in  Section  V.8.  Subroutines 
INWALLI,  IHTERI,  and  AXISI  are  called  to  calculate  the  flowfield.  Sub¬ 
routine  H0VE1  is  called  to  transfer  data  between  computational  points 
and  solution  arrays,  as  required.  Subroutine  THRUST  is  called  to  inte¬ 
grate  the  thrust  generated  along  the  nozzle  wall,  as  described  in 
Section  IV. 2.  Subroutine  DWDOTL  is  called  to  integrate  the  mass  flow 
rate  and  thrust  along  the  left-running  Mach  lines  as  they  are  construct¬ 
ed.  Subroutine  OUTPUT  is  called  to  write  out  the  flow  properties  at 
each  solution  point. 

MOCRRC.  This  subroutine  is  the  control  subroutine  for  calculating 
the  flowfield  along  right-running  Mach  lines  emanating  from  the  super¬ 
sonic  contour  of  the  basic  nozzle  as  described  in  Section  V.9. 

Subroutines  DRMALU,  INTERI,  and  AXISI  are  called  to  calculate  the  flow- 
field.  Subroutine  MGVE1  is  called  to  transfer  data  between  computational 
points  and  solution  arrays,  as  required.  Subroutine  THRUST  is  called  to 
integrate  the  thrust  generated  along  the  nozzle  wall,  as  described  in 
Section  IV. 2.  Subroutine  OWDOTR  is  called  to  integrate  the  mass  flow 
rate  and  thrust  along  the  right-running  Hach  lines  as  they  are  constructed. 
Subroutine  OUTPUT  is  called  to  write  out  the  flow  properties  at  each 
solution  point. 

NQCtRCI,  This  subroutine  is  the  control  subroutine  for  calculating 
the  flowfield  along  left-running  Mach  lines  that  intersect  the  supersonic 
contour  of  the  basic  nozzle,  as  described  in  Section  V.10.  Subroutines 
DRiiAltl,  INTERI*  and  AXISI  are  called  to  calculate  the  flowfield. 
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Subroutine  M0VE1  Is  called  to  transfer  data  between  computational  points 
and  solution  arrays,  as  required.  Subroutine  THRUST  Is  called  to  Inte¬ 
grate  the  thrust  generated  along  the  nozzle  wall,  as  described  In 
Section  IV. 2.  Subroutine  DWDOTL  Is  called  to  integrate  the  mass  flow 
rate  and  thrust  along  the  left-running  Mach  lines  as  they  are  constructed. 
Subroutine  OUTPUT  is  called  to  write  out  the  flow  properties  at  each 
solution  point. 

RESET.  This  subroutine  resets  points  along  a  right-running  Mach 
line  after  it  crosses  the  previous  right-running  Mach  line,  so  that  the 
remaining  points  to  be  calculated  are  defined  to  be  the  remaining  points 
on  the  previous  right-running  Mach. line,  as  illustrated  in  Figure  39. 

SHIFT.  This  subroutine  resets  points  along  a  left-running  Mach 
line  after  it  crosses  the  previous  left-running  Mach  line,  so  that  the 
remaining  points  to  be  calculated  are  the  remaining  points  on  the  pre¬ 
vious  left-running  Mach  line,  as  illustrated  in  Figure  44. 

HOVE! .  This  subroutine  transfers  data  between  grid  points  (points 
1  to  4)  and  the  storage  arrays  [X(15Q),  Y(15Q),  etc.]  as  requested  by 
the  logic  control  subroutines  (HGCIVL,  N0CARCR,  MOCARCL,  KOCRRC,  and 
HDCIRCI). 

ORtfAULl.  This  subroutine  implements  the  method  of  characteristics 
solution  for  a  direct  wall  point  in  an  irrotational  flowfield,  as 
illustrated  in  Figure  37. 

I  KHALI! .  This  subroutine  implements  the  method  of  characteristics 
solution  for  an  inverse  wall  point  in  an  irrotational  flowfield,  as 
illustated  in  Figures  31  and  32. 
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POINTEE.  This  subroutine  defines  the  properties  at  the  beginning 
of  an  embedded  shock  wave  at  point  (IE.  JE)  in  the  characteristic  coor¬ 
dinate  system.  This  feature  is  used  in  MODE  a  3  and  only  in  MODE  =  3. 


6.  OVERLAY  (3,0) 


MAIN30.  This  program  is  the  main  control  program  for  OVERLAY  (3,0), 
which  calculates  the  rotational  flowfield  downstream  of  the  oblique  shock 
wave.  Subroutine  MOCLRCL  is  called  to  control  the  logic  for  calculating 
the  flowfield.  If  an  attached  oblique  shock  wave  exists  at  point  E 
(MODE  =  4  or  5),  subroutine  POINTE  is  called  first  to  determine  the  prop¬ 
erties  at  the  shock  wave  at  point  E. 

MOCLRCR.  This  subroutine  contains  the  logic  for  constructing  the 
oblique  shock  wave  emanating  from  the  junction  of  the  basic  nozzle  and 
the  nozzle  extension  (see  Section  V.ll). 

Subroutine  POINTE  is  called  to  determine  the  properties  of  the  oblique 
shock  wave  at  point  E  itself.  Then  subroutines  INTERI,  SHOCK,  INTERR, 
and  DRWALLR  are  called  in  sequence  to  determine  the  solution  for  an 
oblique  shock  wave  point,  as  illustrated  in  Figure  48.  Subroutine  THRUST 
is  called  to  calculate  the  thrust  at  the  wall  points.  Subroutine  MOVE2 
is  called  to  transfer  data  between  computational  points  and  storage 
arrays,  as  required.  Subroutine  DWDOTL  is  called  to  integrate  the  mass 
flow  rate  and  thrust  along  the  left-running  Mach  lines  as  they  are  con¬ 
structed.  Subroutine  OUTPUT  is  called  to  write  out  the  flow  properties 
at  each  solution  point.  When  the  shock  wave  reaches  the  downstream 
extent  of  the  upstream  flowfield,  as  illustrated  in  Figure  53  or  54,  the 
solution  is  terminated  and  the  properties  at  the  exit  of  the  scarfed 
nozzle  extension  are  extrapolated. 

SHIFT2.  This  subroutine  resets  points  along  a  left-running  Mach 
line  after  it  crosses  the  previous  left-running  Mach  line,  so  that  the 
remaining  points  to  be  calculated  are  the  remaining  points  on  the  previous 
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left-running  Mach  line,  as  illustrated  in  Figure  44. 

M0VE2.  This  subroutine  transfers  data  between  grid  points  (points 
1  to  5)  and  the  storage  arrays  [X(150),  Y(15Q),  etc.]  as  required  by  the 
logic  control  subroutine  MOCLRCR  and  the  other  computational  subroutines. 

POINTE.  This  subroutine  calculates  the  properties  of  the  downstream 
side  of  the  oblique  shock  wave  at  point  E,  as  described  in  Section  V . 1 1 . 

SHOCK.  This  subroutine  controls  the  logic  for  calculating  the  sol- 
ution  at  a  shock  wave  point,  as  discussed  in  Section  V.ll. 

SHOCKU.  This  subroutine  calculates  the  properties  on  the  upstream 
side  of  the  oblique  shock  wave  as  discussed  in  Section  V.ll  and  illu¬ 
strated  in  Figures  49  and  50. 

SHOCKO.  This  subroutine  calculates  the  properties  on  the  downstream- 
side  of  the  oblique  shock  wave  as  discussed  in  Section  V.ll  and  illu¬ 
strated  in  Figures  49  and  50. 

RRCHAR.  This  subroutine  applies  the  compatibility  relation  along  a 
rearward  projected  right-running  Mach  line  from  the  downstream  side  of 
the  oblique  shock  wave  to  check  for  overall  convergence  of  the  oblique 
shock  wave  calculation  procedure.  This  procedure  is  described  in  Section 
V.ll  and  illustrated  schematically  in  Figures  49  and  50. 

DRWALLR.  This  subroutine  implements  the  method  o?  characteristics 
solution  for  a  direct  wall  point  in  a  rotational  flow-field,  as  illustrated 
schematically  in  Figure  48. 

INWALL?,.  This  subroutine  implements  the  method  of  characteristics 
solution  for  an  inverse  wall  poiRt  in  a  rotational  flowfield,  as  discussed 
in  Section  V. 12  and  illustrated  schematically  in  Figure 


168 


INTERR.  This  subroutine  implements  the  raet^^f  characteristics 
solution  for  an  interior  point  in  a  rotational  tsm^isid,  as  discussed  in 
Section  V.12  and  illustrated  schematically  in  Figure  57. 

AX I SR.  This  subroutine  implements  the  method  of  characteristics 
solution  for  an  axis  point  in  a  rotational  flowfield.  In  the  present 
program,  right-running  Mach  lines  from  the  nozzle  extension  never  reach 
the  nozzle  axis,  as  discussed  in  Section  V.ll  and  illustrated  schemati¬ 
cally  in  Figures  53  and  54.  Consequently r.  subroutine  AXISR  is  not  called 
in  the  present  program. 
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7.  OVERLAY  (4,0) 

MAIN40.  This  program  is  the  main  control  program  in  OVERLAY  (4,0). 
The  purpose  of  OVERLAY  (4,u)  is  to  implement  the  calculation  of  the 
scarfed  nozzle  performance,  as  discussed  in  Section  V.13.  MAIN40  calls 
subroutine  GEOM  to  calculate  the  angle  i p,  given  by  equation  (121).  Sub¬ 
routine  F0RM0M  is  called  to  calculate  the  performance  of  the  scarfed 
nozzle  and  che  missile  performance. 

GEOM.  This  subroutine  calculates  the  angle  \pt  given  by  equation 

021). 

F0RM0M.  This  subroutine  calculates  the  performance  of  the  scarfed 
nozzle  and  the  overall  missile  performance,  as  discussed  in  Section  V.13. 


SECTION  VII 


INPUT  PARAMETERS 


1 .  INTRODUCTION 

The  problem  specifications  are  defined  by  input  parameters  which  are 
read  in  by  namelists  DATA  and  WALL  in  subroutine  VNPUT  and  by  namelist 
IVSL  in  subroutine  IVLTAB,  all  in  OVERLAY  (1,0).  Only  those  parameters 
that  are  pertinent  for  a  particular  problem  need  be  read  in.  Many  para¬ 
meters  have  default  values,  and  need  not  be  specified  unless  values 
different  from  the  default  values  are  to  be  considered.  AH  of  the 
input  parameters  are  discussed  In  this  section.  Where  appropriate, 
default  values  are  specified.  The  default  values  specify  Sample  Case 
No.  1. 

2.  TITLE  CARD 

The  first  card  of  each  data  deck  is  a  title  card  consisting  of  80 
alphanumeric  characters  of  identifying  information.  This  card  may  be 
blank,  or  contain  any  combination  of  allowable  FORTRAN  characters.  This 
card  must  be  the  first  card  of  every  data  deck,  even  if  the  card  is 
blank.  The  format  of  the  card  is  (8A10). 


3.  NAMELIST  DATA 

The  parameters  specified  by  namelist  DATA  are  described  in  this 
section.  All  of  these  parameters  have  default  values,  which  are  reset 
to  their  original  values  before  each  data  case.  The  first  14  parameters 
(IUNITS  to  NWRITE)  are  logic  control  parameters.  The  next  two  para¬ 
meters  (G  and  RG)  specify  the  gas  thermodynamic  model.  The  following 
three  parameters  (PS,  TS,  and  PA}  define  the  nozzle  operating  conditions. 
The  next  tan  parameters  (DELTA  to  XMAX)  specify  the  geometry  of  the 
basic  nozzle.  The  following  three  parameters  (AF,  XF,  and  BETA)  specify 
the  geometry  of  the  scarfed  conical  extension.  The  next  two  parameters 
(DVRATIO  and  DARATIO)  are  used  to  refine  the  spacing  of  the  character¬ 
istic  network.  The  final  five  parameters  (ICORI  to  IDUMP}  are 
convergence  control  parameters  and  dump  flags. 

There  are  actually  eight  more  parameters  included  in  namelist  DATA, 
These  are  ICMP,  IPLOT,  ICPLOT,  YMAX ,  XDES,  YDES,  IE,  and  JE,  These 
paramccr-^  are  not  used  in  the  analysis  of  a  conventional  nozzle  (MODE  - 
1  or  2)  or  in  the  analysis  of  a  scarfed  nozzle  (MODE  =  4  or  5).  They 
are  used  in  the  analysis  of  a  compressed  nozzle  with  an  embedded  right¬ 
running  oblique  shock  wave.  Those  eight  parameters  are  not  described  in 
this  report. 

IUNITS  A  positive  integer  variable  denoting  the  unit  system  employed  in 
the  analysis. 

IUNITS  Unit  system 

1  English  engineering  (EE)  units 


System  International  (SI)  units 


The  default  value  of  IUNITS  is  1. 

A  positive  integer  variable  denoting  the  type  of  analysis  to  be 


performed. 


M-OPE  Type  of  analysis  to  be  performed 

0  Transonic  initial-value  line  only 

1  Conventional  nozzle  following  right-running  Mach 

lines 

2  Conventional  nozzle. following  left-running  Mach 

1  i  nes 

3  Conventional  nozzle  with  an  embedded  RR  shock  wave 

4  Scarfed  nozzle  following  RRCS  and  l.RCS 

5  Scarfed  nozzle  following  LRCS 
The  default  value  of  Hd>DE  is  4. 

IWALL  A  positive  integer  variable  denoting  the  type  of  basic  nozzle 
contour  to  be  analyzed. 


IWALL 


Type  of  basic  nozzle  contour 
Conical  contour 
Quadratic  contour 
Tabular  contour  read  in  from  TAPE5 
Tabular  contour  read  in  from  TAPE8 


The  default  value  of  IWALL  is  1. 

JWALL  A  positive  integer  variable  which  controls  punching  out  the 
basic  nozzle  contour. 


JWALL 


Basic  nozzle  contour  punch  option 
Don't  punch 
Do  punch 
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The  default  value  of  JUALL  is  0. 

NI  A  positive  Integer  variable  denoting  the  nusber  af  points  on 
the  initial-value  line  across  the  nozzle  throat  The  default 
va'iye  of  NI  is  11. 

KT  A  positive  integer  variable  denoting  the  r.umi>£r  of  indirect  wall 
points  on  the  nozzle  throat  downstream  circular  arc  contour. 

The  default  yalue  of  NT  is  15. 

IVS  A  positive  Integer  variable  dennttny  the  type  of  initial-value 
line  to  be  employed. 

IVS  Type  of  initial -value  line 

1  KliegeVs  analysis  (internally  generated) 

b  Tabular  data  read  in  fmi  TAPE5 

7  Tabular  d*ta  read  in  fr«n  TAPE? 

The  default  value  of  IVS  is  1. 

IOUT  A  positive  integer  variable  which  controls*  the  width  of  the 
output. 

TOUT  Width  of  output 

0  80  column  output 

1  132  column  output 

IWRITE  A  positive  integer  variable  specifying  the  right-running 

characteristic  on  which  KWRITE  is  set  to  2,  The  default  value 
of  IWRITE  is  0. 

JWRITE  A  positive  integer  variable  specifying  the  loft-running  charac¬ 
teristic  on  which  KWRITE  is  set  equal  to  2.  The  default  value 
of  JWRITE  is  0. 
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KWRITE  A  positive  integer  variable  which  controls  the  amount  of  flow- 
field, output. 

KWRITE  Output  option 

1  Write  out  only  first  and  last  points  on  each  Mach 
■"  line 

2  Write  out  all  points  on  each  Mach  line,  including 
deleted  points. 

3  Write  out  all  points  on  each  Mach  line,  excluding 
deleted  points 

4  Write  out  only  the  first,  second,  and  last  pages 
The  default  value  of  KWRITE  is  1. 

LWRITE  A  positive  integer  variable  specifying  the  right-running 

characteristic  at  v\*hich  all  of  the  10  dump  flags  IDUMP(I)  are 
set  equal  to  1,  Th*  default  value  of  LWRITE  is  0. 

MWRITE  A  positive  integer  variable  specifying  the  left-running 

characteristic  at  which  all  10  dump  flags  IDUMP(I)  are  set 
equal  to  1.  The  default  value  of  MWRITF  is  0. 

G  The  gas  specJfic  heat-  ratio  y  (dimensionless).  The  default 

value  of  C  is  » .2. 

RG  The  gas  constant  R  (fi-sbf/Tbm-R  or  fcO/kg-K).  The  default  value 
of  RG  is  65.0. 

PS  The  nozzle  inlet  stagnation  pressure  P^  (lbf/in.2  or  N/m^l.  The 
default  value  of  Ps  is  1000.0. 

TS  The  nozzle  inlet  stagnation  temperature  i't  (R  or  K) .  The  default 
value  of  TS  is  5000.0. 
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PA  The  ambient  (i.e.,  atmospheric)  pressure  Pa  (lbf/1n.2  or  N/ra2). 
The  default  value  of  PA  Is  14,696. 

The  nozzle  geometric  model  Is  Illustrated  In  Figure  65.  The  nozzle 
throat  consists  of  a  double  circular  arc  contour.  The  supersonic  contour 
of  the  basic  nozzle  attaches  smoothly  to  the  throat  downstream  circular 
arc  contour  (i.e.,  at  point  A,  ya  and  0a  for  the  circular  arc  and  the 
supersonic  contour  are  the  same).  The  length  of  the  basic  nozzle  is  xe. 
The  basic  nozzle  contour  may  be  conical  (IWALL  =  1).  The  conical 
contour  may  be  specified  by  any  of  the  following  five  sets  of  parameters. 

1.  The  throat  attachment  angle  6a  and  the  nozzle  length  xg. 

2.  The  throat  attachment  angle  0a  and  the  exit  lip  radius  ye, 

3.  The  throat  attachment  angle  0&  and  the  nozzle  area  ratio  e, 

4.  The  nozzle  legth  xe  and  exit  lip  radius  ye. 

5.  The  nozzle  length  x.  and  nozzle  area  ratio  e. 

The  basic  nozzle  contour  may  be  a  quadratic  (IWALL  ®  2).  The 
quadratic  contour  may  be  specified  by  any  of  the  following  three  sets  of 
parameters. 

1.  The  throat  attachment  angle  ©a,  the  exit  lip  angle  0e,  and  the  nozzle 
length  xe. 

2.  The  throat  attachment  angle  6a,  the  nozzle  length  xe,  and  the  exit 
lip  radius  ye. 

3.  The  throat  attachment  angle  ea,  the  nozzle  length  xe,  and  the  nozzle 
area  ratio  e. 

The  default  values  of  all  these  parameters  are  0.0.  The  desired 
option  is  selected  simply  by  specifying  nonzero  values  for  the  appro¬ 
priate  variables.  The  remaining  variables  are  internally  calculated 
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(Vyt> 


from  those  specified. 

The  remaining  option  for  specifying  the  supersonic  contour  of  the 
basic  nozzle,  the  tabular  wall  option,  is  discussed  in  Section  VII. 4. 

The  conical  extension  is  completely  specified  by  the  cone  angle  6f 
and  the  x  coordinate  of  the  end  of  the  extension,  Xf.  The  conical 
extension  is  attached  to  the  basic  nozzle  at  point  E. 

DELTA  A  real  variable  denoting  the  type  of  coordinate  system  to  be 
considered. 

DELTA  Coordinate  system 

0.0  Two-dimensional  planar 

1.0  Two-dimensional  ax i symmetric 

The  default  value  of  DELTA  is  1.0. 

YT  The  nozzle  throat  radius  yt  (In,  or  m).  The  default  value  of  YT 
is  1.0. 

RTU  The  nozzle  throat  upstream  circular  arc  radius  of  curvature 
(in.  or  m).  The  default  value  of  RTU  is  1.0. 

RTD  The  nozzle  throat  downstream  circular  arc  radius  of  curvature 
(in.  or  m).  The  default  value  of  RTD  is  0.01. 

AA  The  nozzle  throat  attachment  angle  ea  (deg).  The  default  value 
of  AA  is  0.0. 

AE  The  basic  nozzle  exit  lip  angle  at  point  E,  6e  (deg).  The 
default  value  of  AE  is  0.0. 

EPS  The  basic  nozzle  area  ratio  e  B  (ye/yt)^  (dimensionless).  The 
default  value  of  EPS  is  0.0. 

YE  The  basic  nozzle  exit  lip  radius  ye  (in.  or  m).  The  default 
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value  of  YE  is  0.0. 

XE  The  x  coordinate  of  the  length  of  the  basic  nozzle  xe  (in.  or  m). 
The  default  value  of  XE  is  O.Q. 

XMAX  A  real  variable  denoting  the  x  location  at  which  right-running 
Mach  lines  may  be  terminated  to  save  computational  effort  when 
the  flowfield  beyond  the  end  of  the  no.’zle  is  not  required*  as 
illustrated  in  Figure  66.  XMAX  is  only  used  when  analyzing  a 
basic  nozzle  following  right-running  Mach  lines  (MODE  =1),  In 
that  case,  XMAX  must  be  at  least  as  large  as  XE,  the  length  of 
the  basic  nozzle.  For  the  analysis  of  scarfed  nozzles  (MODE  =  4 
or  5),  XMAX  must  be  a  large  number  so  that  the  entire  flowfield 
from  the  basic  nozzle  is  obtained.  The  default  value  of  XMAX  is 
1000.0. 

The  conical  nozzle  extension  is  also  illustrated  in  Figure  65.  The 
conical  extension  may  be  specified  by  either  of  the  following  options. 

1.  The  cone  angle  9f  and  the  x  coordinate  of  the  end  of  the  conical 
extension  x^. 

2.  The  cone  angle  9f  and  the  nozzle  scarfing  angle  0. 

The  default  values  of  Xf  and  6  are  both  0.0.  The  desired  option  is 
selected  by  simply  specifying  a  nonzero  value  for  either  x^  or  6.  The 
other  variable  is  internally  computed  from  the  one  specified. 

AF  The  conical  extension  wall  angle  0f  (deg).  The  default  value  :v; 
AF  is  0.0. 

XF  The  x  coordinate  of  the  end  of  the  conical  extension  Xf  (/*r 
m).  The  default  value  of  XF  is  0.0. 
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Reduced  flowfield 


BETA  The  nozzle  scarfing  angle.  The  default  value  of  BETA  is  0.0. 

The  next  two  variables,  DYRATlO  and  DARATIO,  are  used  to  refine  the 
spacing  of  the  characteristic  network.  DYRATlO  is  discussed  in  Section 
V.3,  and  DARATIO  is  discussed  in  Section  V.5. 

DYRATlO  A  real  variable  specifying  the  ratio  of  the  final  Ay  on  the 

initial-value  line  adjacent  to  the  nozzle  wall  to  the  initial  Ay 
adjacent  to  the  nozzle  axis  (dimensionless).  The  default  value 
of  DYRATlO  is  1.0. 

DARATIO  A  real  variable  specifying  the  ratio  of  the  final  AO  on  the 
nozzle  throat  downstream  circular  arc  adjacent  to  the  throat 
attachment  point,  point  A  in  Figure  37,  to  the  initial  AO  adja¬ 
cent  to  the  nozzle  throat,  point  T  in  Figure  37  (dimensionless). 
The  default  value  of  DARATIO  is  1.0. 

The  remaining  five  parameters  are  convergence  control  parameters 
and  dump  flags. 

ICOR  A  positive  integer  variable  specifying  the  number  of  applications 
of  the  corrector  in  the  modified-Euler  predictor-corrector 
numerical  solution  of  the  characteristic  equations  and  compati¬ 
bility  relations.  ICOR  s  0  yields  the  Euler  predictor  method, 
ICOR  =  1  yields  the  modified-Euler  predictor-corrector  method, 
and  ICOR  >  2  yields  the  modified-Euler  predictor-corrector  method 
with  (ICOR  -  1)  iterations.  The  default  value  of  ICOR  is  2, 

El  A  positive  real  variable  denoting  the  tolerance  (in.  or  m)  for 
location  convergence  in  the  numerical  method  of  characteristics 
calculations.  The  default  value  of  El  is  0.0,  causing  the  unit 
process  calculations  to  terminate  after  ICOR  applications  of  the 


Integration  algorithm, 

E 2  A  positive  real  variable  denoting  the  fractional  convergence 

tolerance  (dimensionless)  for  flow  property  convergence  in  the 
numerical  method  of  characteristics.  The  default  value  of  E2  is 
0.0,  causing  the  unit  process  calculations  to  terminate  after 
ICOR  applications  of  the  integration  algorithm. 

I0UHP  A  one-dimensional  Integer  variable  array  dimensioned  at  10,  Each 
element  of  IDUMP  activates  output  dumps  of  selected  parameters 
during  the  flowfield  calculations  for  use  in  debugging  problems, 

A  value  of  0  suppresses  the  dumps*  and  a  value  of  1  activates  the 
du&ps.  The  default  values  of  all  10  elements  of  IDOHP  are  0. 
IPUMP  (I)  Subroutines  with  dumps  activated 

1  H0CLRCI,  HAIN30,  H0CLRCR 

2  RESET,  SHIFT,  SHIFTS 

3  POINTEE.  POINTE,  SHOCK,  SK0CKU,  SHOCKS 

4  C0MPRES,  ATTACH 

5  RRCHAR 

6  BOUNOYW,  BGUNDYE,  OR WALL I,  0RWALLR 

7  inter;,  interr 

8  INWALL I,  INWALLR 

9  AXISI,  AXISR 

10  not  used 

This  completes  the  specification  of  all  of  the  input  data  read  in  by 
namelist  DATA. 
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4.  NAMELIST  WALL 

The  parameters  specified  by  namelist  WALL  are  described  in  this 
section.  These  parameters  specify  the  tabular  wall  contour  option 
(IWALL  =  5  or  8  in  namelist  DATA)  for  the  basic  nozzle.  For  the  analy¬ 
tical  wall  options  (IWALL  =  1  or  2),  this  namelist  should  be  omitted 
from  the  data  deck.  The  tabular  nozzle  option  is  discussed  in  Section 
II. 3. d. 

The  tabular  nozzle  geometry  is  illustrated  in  Figure  67,  The  wall 
table,  when  complete,  will  have  NWALL  pairs  of  x  and  y  coordinates  in 
the  arrays  XW  and  VW,  respectively.  The  first  point  in  the  table  will  be 
the  throat  attachment  point,  point  A.  The  coordinates  of  point  A  (i.e., 
xa  and  ya)  are  calculated  internally  and  stored  in  XW(1)  and  YW(1). 
Consequently,  the  values  of  and  8a  (i.e.,  RID  and  AA)  in  namelist 
DATA  must  be  specified.  The  remaining  wall  points  are  read  in  through 
namelist  WALL.  The  first  input  point  is  XW(2)  and  YW(2).  The  succeeding 
points  follow  in  order.  The  last  point  specified  in  the  tabular  data  is 
defined  internally  to  be  point  E  (i.e.,  x0  and  y0).  The  number  of  points 
actually  read  in  through  namelist  WALL  is  denoted  by  NWALL.  The  total 
number  of  points  in  the  wall  table  is  NWALL  +  1,  due  to  the  insertion  of 
point  A  as  the  first  point  in  the  table.  The  input  value  of  NWALL  is 
increased  by  one  internal  to  the  program  to  reflect  the  total  number  of 
points  in  .  r  wall  table  The  computer  program  locally  fits  a  quadratic 
polynomial  to  the  tabular  data  when  calculating  the  intersection  point  of 
a  left-running  Kach  line  ar.d  the  tabular  nozzle  wall. 


NWALL  An  Integer  variable  specifying  the  number  of  points  in  the  tab¬ 
ular  wall  table  (must  be  99  or  less).  No  default  value  is 
specified  for  NWALL. 

XW  A  one-dimensional  real  variable  array  dimensioned  at  99.  XW(I) 

K  (I  8  1,  NWALL)  denotes  the  axial  location  (in.  or  m)  of  the  lift 
tabular  wall  point.  No  default  values  are  specified  for  XW(I). 

YW  A  one-dimensional  real  variable  array  dimensioned  at  99.  YW(I) 
(I  »  1,  NWALL)  denotes  the  radial  location  (in.  or  m)  of  the  1th 
tabular  wall  point.  No  default  values  are  specified  for  Ytt(I), 
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5.  NAMELIST  IVSL 

The  parameters  specified  by  namelist  IVSL  are  described  in  this 
section.  These  parameters  specify  the  tabular  initial -value  line  option 
(IVS  a  5  or  7  in  namelist  DATA).  For  the  internally  generated  initial- 
value  line  option  (IVS  =  1),  this  namelist  should  be  omitted  from  the 
data  deck.  The  tabular  initial -valurltne  option  is  discussed  in  Section 
V.J.B. 

The  tabular  initial-value  line  is  illustrated  in  Figure  68.  The 
f>rst  point  on  the  tabular  initial-value  line  must  be  the  nozzle  axis 
point,  and  the  last  point  must  be  the  nozzle  throat  point,  point  T.  The 
number  of  points  NI  (specified  in  namelist  DATA)  and  their  spacing  is 
arbitrary.  The  maximum  number  of  tabular  initial-value  line  points  is 
30. 

XIV  A  one-dimensional  real  variable  array  dimensioned  at  30.  XIV(I) 

(1  =  1,  NI)  denotes  the  axial  location  (in.  or  m)  of  the  Ith^ 
tabular  initial -value  line  point.  No  default  values  are  speci¬ 
fied  for  XIV(I). 

RIV  A  one-dimensional  real  variable  array  dimensioned  at  30.  RIV(I) 

(I  “  1,  NI)  denotes  the  radial  location  (in.  or  m)  of  the  Itn 
tabular  initial-value  line  point.  No  default  values  are  specified 
for  RIV(I) . 

HIV  A  one-diwensional  real  variable  array  dimensioned  at  30.  HIV(I) 

(1  =  1,  NI)  denotes  the  Nach  number  (dimensionless)  at  the  Ith 
point  on  the  tabular  initial-value  line.  No  default  values  are 
specified  for  HIV(I). 
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TIV  A  one-dimensional  real  variable  array  dimensioned  at  30.  TIV(I) 
(I  «  1,  Nil  denotes  the  flow  angle  (deg)  at  the  Ith  point  on  the 
tabular  initial-value  line.  No  default  values  are  specified  for 
TIV  Cl ) * 


c.  i? 


*1  'r 


r  •  •• 


•  ' .  •  V 


«.**•*■. 


168 


•  * 

*  ►  *.•  ' 
,‘V- 


/  s 


SECTION  VIII 
SAMPLE  CASES 


1.  INTRODUCTION 

Fifteen  sample  cases  are  presented  in  this  section  to  illustrate  the 
application  of  the  computer  program  for  calculating  the  flowfield  in 
scarfed  propulsive  nozzles.  For  each  sample  case,  a  discussion  of  the 
problem  is  presented,  the  input  data  are  discussed,  and  the  data  deck 
required  to  execute  the  program  is  presented.  Selected  portions  of  the 
computer  output  are  presented  for  some  of  the  sample  cases. 

The  first  eight  sample  cases  are  all  concerned  with  the  same  scarfed 
nozzle  configuration.  That  nozzle  consists  of  a  conical  nozzle  and  a 
cylindrical  extension.  The  nominal  case  is  specified  by  a  nozzle  cone 
angle  0a  a  15  deg,  a  nozzle  area  ratio  e  -  10,  a  cylindrical  extension 
with  6f  =  0  deg,  and  a  scarf  angle  0  "  30  deg. 

Sample  Case  No.  1  includes  a  complete  discussion  of  the  problem  and 
the  input  data  for  the  case  in  which  the  conical  contour  is  specified  by 
the  cone  angle  0a  and  the  area  ratio  e.  Sample  Cases  No.  2  to  4  are 
identical  to  Sample  Case  No.  1,  except  that  the  conical  contour  is 
specified  by  the  exit  lip  radius  ye,  the  .sngth  xe,  and  tabular 

data,  respectively.  Sample  Case  No.  5  illustrates  the  tabular  initial- 
value  line  option.  Sample  Case  No.  6  illustrates  an  optional  Mach  Vine 
network  construction  feature.  Sample  Case  No.  7  illustrates  the 
variable  initial-value  line  spacing  option  and  the  variable  inverse  wail 


point  spacing  option  on  the  nozzle  throat  downstream  circular  arc. 

Sample  Case  No.  8  is  the  same  as  Sample  Case  No.  1,  except  that  the 
length  Xf  of  the  conical  extension  is  specified  instead  of  the  scarf 
angle  8.  Sample  Case  No.  9  has  the  same  nozzle  as  Sample  Case  No.  1,  but 
the  scarfed  nozzle  extension  is  a  cone  with  6^  c  5  deg. 

Sample  Uses  No.  10  to  12  illustrate  various  options  for  analyzing 
a  simple  conical  nozzle  without  an  extension.  Sample  Case  No.  10  con¬ 
siders  the  nominal  case  modified  by  specifying  the  angle  Of  of  the 
conical  extension  to  be  the  same  as  the  nozzle  cone  angle  6a.  Sample 
Case  No.  11  analyzes  the  combined  nozzle  and  extension  considered  in 
Sample  Case  No.  10  as  a  single  conical  nozzle  without  an  extension. 

Sample  Case  No.  12  illustrates  an  optional  Mach  line  network  construction 
feature  for  a  single  nozzie.  This  same  feature  is  illustrated  in  Sample 
Case  No.  6  for  the  nominal  case  with  a  cylindrical  extension. 

Sample  Cases  No.  13  to  15  are  all  concerned  with  the  same  scarfed 
nozzle  configuration.  That  nozzle  consists  of  a  quadratic  nozzle  contour 
and  a  cylindrical  nozzle  extension.  Sample  Case  No.  13  is  the  nominal 
case  for  the  quadratic  nozzle  contour.  The  nominal  case  is  specified  by 
a  throat  attachment  angle  9a  s  25  deg,  an  area  ratio  e  °  10,  and  a 
nozzle  length  xe  3  8.07104661.  The  scarfed  nozzle  extension  is  cylindri¬ 
cal  and  the  scarf  angle  6  °  30  deg.  This  nozzle  has  the  identical 
envelope  as  the  conical  nozzle  nominal  case.  Sample  Case  No.  1.  Sample 
Cases  No.  14  and  15  are  identical  to  Sample  Case  No.  13,  except  that  the 
quadratic  nozzle  contour  is  specified  by  the  nozzle  exit  lip  radius  ye  and 
the  nozzle  exit  lip  angle  9e,  respectively. 
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2.  SAMPLE  CASE  NO.  1 

This  sample  case  is  the  basic  problem  considered  in  all  of  the 
sample  cases.  The  nozzle  to  be  analyzed  is  a  conical  nozzle  having  a 
scarfed  cylindrical  extension  as  illustrated  in  Figure  69.  The  English 
Engineering  (EE)  system  of  units  is  employed  in  all  of  the  sample  cases 
(IUNITS  =  1).  The  problem  to  be  considered  is  the  analysis  of  a  scarfed 
propulsive  nozzle  (K(DE  =  4).  The  basic  nozzle  is  a  conical  nozzle 
(IWALL  *  1).  Punched  output  of  the  wall  contour  is  not  desired  (OWALt  » 
0).  Eleven  equally  spaced  points  are  desired  on  the  initial-value  line 
(NI  =  11  and  DYRATIO*  1.0),  and  15  equally  spaced  points  are  desired  on 
the  nozzle  throat  downstream  circular  arc  (NT  =  15  and  DARATIO®  1.0). 

The  initial-value  line  is  to  be  generated  internally  by  Kliegel's  method 
CIVS  =1).  Every  point  on  each  Mach  line  is  to  be  written  out  (KWRITE  ® 
2),  and  KWRITE  is  to  remain  at  the  value  of  2  for  the  entire  calculation 
(IWRITE  «  JWRITE  «  0). 

The  gas  specific  heat  ratio  y  is  1.2  (GAftiA  »  1.2)  and  the  gas  con¬ 
stant  R  is  65.0  (ft-lbf)/(lbra-R). 

The  nozi-lw  inlet  stagnation  pressure  and  temperature  T^  are 
100U.9  lbf/in.2  and  5000.0  R,  respectively  (PS  -  1000.0  and  TS  =  5000.0). 
The  ambient  pressure  Pa  is  14.696  lbf/in.^  (PA  =  14.696). 

The  nozzle  is  axisyuaetric  (DELTA  =  1.0),  The  throat  radius  y^  is 
1.0  in.  (VT  *  1.0),  the  nozzle  throat  upstream  radius  of  curvature  is 
1.0  in,  (RTU  »  1.0),  and  the  nozzle  throat  downstream  radius  of  curvature 
Pkj  is  0.01  in.  (RTO  -  0.01).  The  basic  nozzle  contour  is  conical  (IWALL 
-  1),  and  the  conical  geometry  is  specified  by  the  nozzle  throat 
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attachment  angle  ea  of  15  degrees  (AA  =  15.0)  and  an  area  ratio  e  =  10 
(EPS  -  10.0).  The  parameters  0a  and  e  completely  specify  the  geometry 
of  the  conical  nozzle.  Hence,  the  exit  lip  radius  ya  and  the  nozzle 

c 

length  xe  are  left  at  their  default  values  of  0.0  (YE  =  0.0  and  XE  = 

0.0).  The  complete  flowfield  from  the  basic  nozzle  must  be  calculated 
since  a  scarfed  extension  is  present.  Thus,  XMAX  is  left  at  its  default 
value  of  1000.0. 

The  cone  angle  6f  of  the  scarfed  conical  extension  is  0,0  deg  (AF  = 
0.0).  The  scarf  angle  B  =  30.0  deg  (BETA  =  30.0).  The  parameters  Of  and 
B  completely  specify  the  geometry  of  the  scarfed  conical  extension. 

Hence,  the  length  x*  or  the  scarfed  conical  extension  is  left  at  its 
default  value  of  0.0  (XF  =  0.0). 

Two  applications  of  the  corrector  (ICOR  =  2)  are  to  be  made  in  the 
modified-Euler  predictor-corrector  numerical  integration  procedure  used 
in  the  method  of  characteristics  unit  processes.  All  flowfield  points 
are  to  be  corrected  two  times,  so  the  convergence  tolerances  are  set 
equal  to  0.0  (El  =  0.0  and  E2  =  0.0). 

This  completes  the  specification  of  namelist  DATA.  Since  the 
initial-value  line  and  the  basic  nozzle  contour  are  internally  generated 
(IVS  3  1  and  IWALL  =1),  name’ists  IVSL  and  WALL  are  omitted  from  the 
data  deck.  The  data  deck  for  Sample  Case  No.  1  is  presented  below. 

SAMPLE  CASE  NO.  1.  NOMINAL  CASE  WITH  AA=15,  EPS=10,  AF=0,  AND  BETA=30 
SDATA  MODE =4,  rUALl-1 ,  NI=11,  NT=15,  KWR1TE«2. 

G=1 .2,  RG=65.0,  P$=1000.Q.  TS=5000.0,  PA=14.696, 

YT=1 .0,  RTU=1 .0,  RTB=0.01 ,  AA=15.0,  EPS=10.0,  AF=0.0,  BETA=30.0  $ 
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The  computer  output  for  Sample  Case  No.  1  Is  presented  in  Append.  A. 
The  first  page  presents  the  program  abstract,  the  job  title,  the  problem 
specifications,  and  the  thermodynamic  model.  The  second  page  describes 
the  nozzle  geometry.  The  next  two  pages  present  a  complete  listing  of 
namelist  DATA.  The  fifth  page  presents  the  initial-value  line  proper-  ■ 
ties  and  performance  parameters. 

The  three  following  pages  present  the  flow  properties  at  each  point 
along  the  right-running  Nach  lines  emanating  fro®  the  initial-value  line, 
as  illustrated  in  Figure  36.  The  first  point  in  each  set  of  data  is  the 
initial-value  ’ine  point.  That  point  is  followed  by  the  succeeding 
points  along  the  right-running  Hach  line.  The  last  point  in  each  set  of 
data  is  the  point  where  the  right-running  Hach  line  intersects  the  nozzle 
axis.  Since  11  points  were  specified  on  the  initial-value  line  (NI  •*  11), 
11  right-running  characteristics  are  present  in  this  initial  portion  of 
the  flowfield  (I  «  1  to  11). 

The  next  nine  pages  present  the  flow  properties  at  each  point  along 
the  right-running  Kaeh  lines  emanating  from  the  prespecified  points  on 
the  initial  expansion  contour,  that  Is,  the  circular  arc  throat  down¬ 
stream  radius  of  curvature,  as  illustrated  in  Figure  40,  Since  15  points 
were  specified  along  the  initial -expansion  contour  (NT  «  15),  15  right¬ 
running  characteristics  are  present  in  this  portion  of  the  flowfield 
(I  c  12  to  26).  The  nozzle  performance  parameters  at  each  prespecified 
wall  point  are  listed  ahead  of  the  corresponding  Kach  line  data. 

The  next  10  pages  present  the  flow  properties  at  ti  a  point  along 
the  left-running  characteristics  emanating  from  points  along  the  final 
right-running  characteristic  from  the  initial -expansion  contour,  as 
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illustrated  in  Figure  50.  The  nozzle  performance  parameters  at  each 
direct  wall  point  at  the  end  of  each  left-running  Mach  line  are  listed 
at  the  end  of  each  set  of  Mach  line  points.  These  left-running  Mach 
lines  are  generated  until  the  end  of  the  basic  nozzle  contour  is  reached. 
In  this  sample  case,  30  left-running  Mach  lines  are  required  (J  =  1  to 
30).  The  last  point  (I  =  55,  J  =  30)  is  the  exit  lip  point  of  the  basic 
nozzle,  point  E. 

The  next.  12  pages  present  the  flow  properties  at  each  point  along 
the  left-running  characteristics  that  reach  the  scarfed  nozzle  extension. 
All  of  these  Mach  lines  pass  through  the  right-running  oblique  shock  wa»e 
emanating  from  the  junction  between  the  basic  nozzle  contour  and  the 
scarfed  nczzle  extension,  point  E.  The  first  point  (I  =  55f  J  =  30)  is 
the  shock  wave  point  at  point  E  itself.  Each  shock  wave  point  has  two 
lines  of  data:  the  first  line  contains  the  flow  properties  on  the 
upstream  side  of  the  shock  wave  and  the  second  line  contains  the  flow 
properties  on  the  downstream  side  of  the  shock  wave.  IHe  nozzle  perform¬ 
ance  parameters  at  each  point  on  the  scarfed  nozzle  extension  at  the  end 
of  each  left-running  Mach  line  are  listed  at  the  end  of  each  set  of  Hach 
line  points.  These  left-running  Ha eh  lines  are  generated  until  the  end 
of  the  scarfed  nozzle  extension  is  reached.  In  this  sample  case,  20  left- 
running  Kacb  lines  are  required  (J  =  31  to  SO).  The  last  point  {!  =  71, 

0  *  SO)  is  the  exit  lip  point  of  the  scarfed  nozzle  extension*  point  E. 

The  next  two  pages  present  a  smssaury  of  the  nozzle  performance  para¬ 
meters  at  each  point  along  the  nozzle  wail  contour  calculated  daring  the 
method  of  characteristics  analysis. 

The  next  page  sueoarizes  the  performance  of  the  scarfed  nozzle 
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extension. 

The  final  page  of  output  presents  a  summary  of  the  overal 1  scarfed 
nozzle  performance  parameters  and  a  summary  of  the  overal 1  missile 
performance  parameters.  The  final  result  of  the  analysis  Is  the  axial 
specific  Impulse  delivered  to  the  missile*  denoted  by  ISPXM.  For  the 
present  sample  case,  this  value  is  206.785  (lbf-see)/lbm. 

This  computer  run  required  approximately  25  sec  of  CPU  time  on  a 
CDC  6500  computer  and  generated  approximately  2100  lines  of  printed 
output. 

In  many  cases,  the  amount  of  output  described  above  Is  not  needed. 

A  considerable  reduction  In  output  is  achieved  by  setting  KWRITE-  Iv¬ 
in  that  case,  every  section  of  output  described  about  Is  still  obtained, 
but  only  the  first,  last,  and  shock  wave  points  are  written  out  along 
each  characteristic.  For  Sample  Case  No.  1  with  KWRITE  ~  1,  approximate 
ely  21  sec  of  CPU  time  were  required  and  approximately  500  lines  of 
printed  output  were  generated.  Even  less  printed  output  is  obtained  by 
setting  KWRITE  «  4.  In  that  case,  only  the  first  two  pages  describing 
the  probl^n  specifications  and  the  last  page  summarizing  the  overall 
performance  parameters  arc  written  out.  For  Sample  Case  No.  1  with 
KWRITE  -  4,  approximately  15  sec  of  CPU  time  were  required  and  approxi¬ 
mately  70  lines  of  printed  output  were  generated.  The  portion  of  the 
output  for  Sample  Case  No.  1  that  is  abbreviated  by  setting  KWRITE  ■  1 
is  presented  In  Append,  B. 


3.  SAMPLE  CASE  NO.  2 

Sample  Case  No.  2  is  identical  to  Sample  Case  No.  1,  except  that  the 
conical  portion  of  the  basic  nozzle  is  specified  by  the  exit  radius  ye 
instead  of  the  cone  angle  6a,  and  the  amount  of  output  is  reduced.  The 
identical  basic  nozzle  contour  is  obtained  for  ye  =  3.16227766  in.  Thus, 
the  only  changes  to  the  input  data  are  to  delete  AA  =  15.0,  to  add  YE  = 
3.16227766,  and  to  change  KWRITE  to  1.  The  data  deck  for  Sample  Case  No. 
2  is  presented  below.  The  output  is  identical  to  the  abbreviated  output 
for  Sample  Case  No.  1. 

SAMPLE  CASE  NO.  2.  NOMINAL  CASE  WITH  YE=3. 16227766  SPECIFIED 
$DATA  W 1WALL=1,  NI=11,  NT=15,  KWRITE=1 , 

6*1.2,  RG=65.0,  PS=1000.0,  TS=5000.0,  PA=14.696, 

YT=1 .0,  RTU=1 ,0,  RTD*0.01,  AA=15.0,  YE=3. 16227766,  AF=0.0,  BETA=30.0  $ 
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4.  SAMPLE  CASE  NO.  3 

Sample  Case  No.  3  is  identical  to  Sample  Case  No.  1,  except  that  the 
conical  portion  of  the  basic  nozzle  is  specified  by  the  nozzle  length  xe 
instead  of  the  cone  angle  Ga.  The  identical  basic  nozzle  is  obtained  for 
xe  =  8.07104661  in.  Thus,  the  only  changes  to  the  input  data  are  AA  = 

0.0  and  XE  ®  8.07104661.  The  data  deck  for  Sample  Case  No,  3  is  present¬ 
ed  below.  The  output  is  identical  to  the  abbreviated  output  for  Sample 
Case  No.  1. 


SAMPLE  CASE  NO,  3.  NOMINAL  CASE  WITH  XE=8. 071 04661  SPECIFIED 
$DATA  M0DE=4,  IWALL=1 ,  NI=11,  NT=15,  KWRITE=1, 

G=1 .2,  RG=65.0,  PS=1000.0,  TS=5000.0,  PA-^14.696, 

YT=1 .0,  RTU=1.0,  RTD=0.01 ,  AA=15.0,  XE=8. 071 04661 ,  AF=0.0,  BETA*30.0  $ 
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5.  CAMPLE  CASE  NO.  4 

Sample  Case  No.  4  is  identical  to  Sample  Case  No..  1,  except  that  the 
supersonic  contour  of  the  basic  nozzle  is  read  in  as  tabular  data  from 
TAPF5  (IWALL=5) .  For  the  tabular  wall  option,  the  basic  nozzle  throat 
contour  must  still  be  specified,  since  the  first  point  in  the  wall  contour 
tabie  is  computed  internally  as  point  A.  Thus,  the  throat  attachment 
angle  must  be  explicitly  specified  as  15.0  deg  (AA=15.0).  The  basic 
nozzle  exit  lip  point,  point  E,  is  internally  defined  to  be  the  last  point 
in  the  wall  table.  Thus,  xe  and  ye  or  e  (i.e.,  XE,  YE,  and  EPS)  do  not 
need  to  be  input.  Thus,  the  only  change  in  the  input  parameters  required 
in  name!* .t  DaTA  is  IWALL=5. 

For  the  tabular  wall  option,  namelist  WALL  must  be  specified.  Only 
two  points  are  required  to  specify  a  conical  wall  (NWALL=2).  The  first 
point  can  La  any  point  downstream  of  the  throat  attachment  point,  point  A. 
For  a  contoured  nozzle,  the  first  point  should  be  close  to  point  A  to 
achieve  an  accurate  representation  of  the  nozzle  wall.  For  a  conical 
wall,  the  location  of  the  first  ooint  is  irrelevant,  as  long  as  it  lies 
between  the  points  A  and  E.  That  point  is  specified  as  (4.03681741, 
2.08130520).  which  is  exactly  midway  between  point*  A  and  E.  The  last 
point  is  point  E  (8.07104661,  3. 16227766) , 

The  data  deck  for  Sample  Case  No.  4  is  presented  below. 


SAMPLE  CASE  NO.  4.  NOMINAL  CASE  WITH  TABULAR  NOZZLE  WALL 
$DATA  M0DE=4,  IWALL=5,  HI-11,  NT*15,  KWRITE-1, 

6=1.2,  R6=65.0,  PS=1000.0,  TS=5000.0,  PA=14.696, 

YT=1.0,  RTU=1.0,  RTD=0.01»  AA=15,0,  AF*0.0,  BETA-3Q.O,  $ 

$WALL  NWALL=2,  XW(2)=4. 03681 741,  8.07104661, 

YW(2)«2. 08130920,  3.16227766  $ 

The  output  for  Sample  Case  No.  4  is  the  same  as  the  abbreviated 
output  for  Sample  Case  No.  1,  except  that  the  tabular  wall  data  are 
written  out  on  the  second  page  of  output.  That  page  Is  presented  In 
Append.  C. 


6.  SAMPLE  CASE  NO.  5 

Sample  Case  No.  5  is  identical  to  Sample  Case  No.  1,  except  that  the 
initial-value  line  is  read  in  as  tabular  data  from  TAPES  (IVS=5).  The 
initial-value  data  are  taken  from  the  output  of  Sample  Case  No.  1.  The 
only  change  in  the  input  required  in  namelist  DATA  is  IVS=5. 

For  the  tabular  initial -value  line  option,  namelist  IVSL  must  be 
specified.  The  point  locations  and  flow  properties  are  taken  directly 
from  the  output  of  Sample  Case  No.  1. 

The  data  deck  for  Sample  Case  No.  5  is  presented  below. 

The  output  for  Sample  Case  No.  5  is  the  same  as  the  abbreviated  out¬ 
put  for  Sample  Case  No.  1. 

SAMPLE  CASE  NO.  5.  NOMINAL  CASE  WITH  TABULAR  INITIAL -VALUE  LINE 
$DATA  M0DE=4,  IWALL=1,  NI=11,  NT=15,  IVS=5,  KURITE=1 , 

6=1.2,  RG=65.0,  PS=1000.0,  TS=5000.G,  PA=14.696, 

YT=1 .0,  RTU=1 .0,  RTD=0.01 ,  AA=15.0,  AF=0.0,  BETA=30.0  $ 

$IVSL  XI (1 )=  0.262,  0.260,  0.252,  0.239,  0.220,  0.197,  0.168, 

0.134,  0.094,  0.050,  0.000, 

YI (1 )=  0.000,  0.100,  0.200,  0.300,  0.400,  0.500,  0.600, 

0.700,  0.800,  0.900,  1.000, 

MI (1 )=  1.031,  1.032,  1.035,  1.041,  1.049,  1.060,  1.075, 

1.095,  1.121,  1.156,  1.203, 

TI(1 )=  0.000,  0.110,  0.210,  0.300,  0.370,  0.410,  0.400, 

0.350,  0.260,  0.130,  0.000  $ 
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7.  SAMPLE  CASE  NO.  6 


Sample  Case  No.  6  illustrates  the  option  of  starting  the  left-running 
characteristic  network  construction  from  the  last  right-running  character¬ 
istic  from  the  initial-value  line  instead  of  from  the  last  right-running 
characteristic  from  the  initial  expansion  contour.  This  option  is  useful 
when  the  shock  wave  intersects  the  last  right-running  characteristic  from 
the  initial  expansion  contour.  This  option  is  specified  by  M0DE=5.  All 
other  parameters  are  identical  to  the  data  for  Sample  Case  No.  1.  The 
data  deck  for  Sample  Case  No.  6  is  presented  below.  The  output  for 
Sample  Case  No.  6  is  presented  in  Append.  D. 

SAMPLE  CASE  NO.  6.  NOMINAL  CASE  WITH  M0DE=5 
$DATA  M0DE=5,  IWALL=1 ,  NI=11,  NT=15,  KWRITE=2, 

G=1 .2,  RG=65.0,  PS=10Q0.0,  TS=5000.0,  PA=14.696, 

YT=! .0,  RTIH.O,  RTD=0.01,  AA=15.0,  EPS=10.0,  AF=0.0,  BETA=30.0  $ 

The  first  six  pages  of  output  for  Sample  Case  No.  6  are  identical  to 
the  output  for  Sample  Case  No.  1.  The  next  two  pages  present  the  solution 
along  the  first  left-running  characteristic  (J=2)  emanating  from  the  down¬ 
stream  extent  of  the  initial -value  line.  The  next  ten  pages  present  the 
solution  along  the  28  left-running  characteristics  (J=  3  to  30)  that  reach 
the  wall  of  the  basic  nozzle.  Although  calculated  in  a  different  order 
from  Sample  Case  No.  1,  the  same  data  points  are  calculated  and  the  same 
solution  is  obtained.  The  remainder  of  the  output  is  identical  to  the 
output  for  Sample  Case  No.  1. 


8.  SAMPLE  CASE  NO.  7 

Sample  Case  No.  7  illustrates  the  use  of  the  variable  step  size 
options  on  the  initial-value  line  and  the  throat  downstream  circular  arc 
contour.  In  all  other  respects.  Sample  Case  No.  7  is  identical  to 
Sample  Case  No.  1.  The  value  of  Ay  on  the  initial -value  line  at  the 
nozzle  wall  is  to  be  one  half  as  large  as  the  value  of  Ay  adjacent  to  the 
nozzle  axis  (DYRATI0=0.5).  The  value  of  AS  adjacent  to  the  throat 
attachment  point,  point  A,  is  to  be  twice  as  large  as  the  value  of  AS 
adjacent  to  the  nozzle  throat  (DARATI0=2.0).  The  data  deck  for  Sample 
Case  No.  7  is  presented  below. 


SAMPLE  CASE  NO.  7.  NOMINAL  CASE  WITH  DYRATI0=0.5  AND  DARATI0=2,0 
$DATA  MQDE=4,  IWALL=1 ,  NI=11,  NT=15,  KWRITE=1,  DYRATI0=0.5,  DARATI0=2.0, 
G-l .2,  RG=65.0,  PS*1000.0,  TS=5GOO.O,  PA=14.696, 


YT=1 .0,  RTtM.O,  RTD=0.Q1,  AA=15.0,  AF=0.0,  BETA=30.0  $ 


9.  SAMPLE  CASE  NO.  8. 

Sample  Case  No.  8  is  the  same  as  Sample  Case  No,  1,  except  that  the 
x  coordinate  of  the  end  of  the  scarfed  nozzle  extension,  point  F,  is 
specified  instead  of  the  scarf  angle  0.  The  same  scarfed  nozzle  extension 
is  obtained  for  Xf  ®  19.02549776  in.  The  only  change  in  the  input  data 
deck  is  the  deletion  of  BETA  =  30.0  and  the  addition  of  XF  =  19.02549776. 
The  data  deck  for  Sample  Case  No.  8  is  presented  below.  The  output  for 
Sample  Case  No.  8  is  the  same  as  the  abbreviated  output  for  Sample  Case 
No.  1. 

SAMPLE  CASE  NO,  8.  NOMINAL  CASE  WITH  XF=1 9.02549776  SPECIFIED 
$DATA  N0DE=4,  IWALL=1,  NI=11,  NT=15,  KWRITE=1, 

G=1 .2,  RG=65.0,  PS=1000.0,  TS=5000.0,  PA=14.696, 

YT=1.0,  RTtKI.O,  RTD=0.01 ,  AA=15.0,  EPS=10.Q,  AF»0.0,  XF=1 9. 02549776  $ 


10.  SAMPLE  CASE  NO,  $ 

Sample  Case  No.  9  considers  the  same  basic  nozzle  and  operating 
conditions  as  Sample  Case  No.  1.  In  Sample  Case  No.  9,  however,  the 
scarfed  nozzle  extension  is  a  cone  with  an  angle  8<c  =  5.0  deg  (AF=5.0). 
The  data  deck  for  Sample  Case  No.  9  is  presented  below. 

SAMPLE  CASE  NO.  9.  NOMINAL  CASE  WITH  AF=5.0 
$DATA  M0DE=4,  IWALL=1,  NI=11,  NT=15,  KWRITE=1, 

6=1.2,  RG=65.0,  PS=1000.0,  TS=5000.0,  PA=14.696, 

YT=1 .0,  RTU=1 .0,  RTD=0.01 ,  AA=15.0,  EPS=10.0,  AF=5,0,  BETA=30.0  $ 

The  output  for  Sample  Case  No.  9  is  identical  to  the  output  of 
Sample  Case  No.  1  through  the  basic  nozzle  flowfield.  The  oblique  shock 
wave  and  the  flowfield  in  the  nozzle  extension  are  different.  Append.  E 
presents  that  portion  of  the  output  downstream  of  the  basic  nozzle. 


11.  SAMPLE  CASE  HO,  10 

Sample  Case  No.  10  illustrates  a  scarfed  conical  nozzle  where  the 
cone  angle  of  the  basic  nozzle  ea  and  the  cone  angle  of*  the  conical 
extension  6|  are  identicals  in  this  case  15,0  deg  (AAa15.0  and  hFb15.0). 
Thus,  the  flow  turning  angle  at  point  E,  the  juction  of  the  basic  nozzle 
and  the  nozzle  extension,  is  zero,  in  that  case,  the  oblique  shock  wave 
emanating  from  point  E  becomes  a  Mach  line,  and  the  flow  properties 
across  the  junction  are  continuous.  The  data  deck  for  Sample  Case  No.  10 
is  presented  below. 

SAMPLE  CASE  NO.  10.  NOMINAL  CASE  WITH  AAaAFal5.0 
$0ATA  MQDEa4,  IWALL=1,  NI=11,  NT=15,  KWRITE*!, 

0=1.2,  RG=65.0,  PS=1000.0,  TS«5000.0,  PAa14.696, 

YT-l.O,  RTUa1.0»  RT0a0.Ql,  AAa15.0,  EPS*10.0,  AF=15.Q,  BETAa3G.O  $ 

The  output  for  Sample  Case  No.  10  is  identical  to  the  output  for 
Sample  Case  No.  1  through  the  basic  nozzle  flowfield.  The  flowfield  in 
the  nozzle  extension  is  different.  Append.  F  presents  that  portion  of 
the  output  downstream  of  the  basic  nozzle.  Note  that  the  oblique  shock 
wave  emanating  from  point  E,  the  junction  point  between  the  basic  nozzle 
and  the  conical  extension,  is  basically  a  Mach  line.  The  strength  of  the 
oblique  shock  wave  does  not  remain  infinitesimal,  as  it  theoretically 
should.  This  discrepancy  is  due  to  numerical  errors  introduced  by  using 
the  rotational  flow  method  of  characteristics  downstream  of  the  pseudo 
shock  wave  to  calculate  an  irrotational  flowfield.  However,  the  errors 
are  relatively  minor. 
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12.  SAMPLE  CASE  NO.  11 

Sample  Case  No.  11  is  identical  to  Sample  Case  No.  1,  except  that 


the  basic  nozzle  has  a  length  of  28.51233074  in.  (XE=28. 51 233074)  and 
there  is  no  nozzle  extension  (M0DE=1).  Consequently,  t,:c  problem  being 
analyzed  is  simply  a  15  deg  conical  nozzle  28.51233074  in.  long.  The 
flowfield  in  this  nozzle  is  identical  to  the  flowfield  in  the  basic 
nozzle  and  nozzle  extension  considered  in  Sample  Case  No.  10.  The  output 
is  slightly  different,  however,  since  the  present  nozzle  does  not  have  a 
special  junction  point  at  x  =  8.C7104661  in.  A  clos-  comparison  of  the 
output  from  the  two  sample  cases  clearly  shows  that  the  flowfields  are 
essentially  identical.  The  data  deck  for  Sample  Case  No.  11  is  presented 
below. 


SAMPLE  CASE  NO.  11.  CONICAL  NOZZLE  CORRESPONDING  TO  CASE  10,  M0DE=1 
$DATA  MODE=l ,  IWALL=1 ,  NI=11,  NT=15,  KWRITE*!, 

Gal .2,  RG=65.0,  PS=10QO.O,  TS-5000.0,  PA*14.696, 

YT=1 .0,  RTU=1 .0,  RTD=0.01,  AA=15.0,  XE=28. 51233074  $ 

The  output  for  Sample  Case  No.  11  is  identical  to  the  output  for 
Sample  Case  No.  1  through  left-running  Mach  line  29.  The  flowfield  down¬ 
stream  of  that  Mach  line  is  the  same  as  the  flowfield  determined  in 
Sample  Case  No.  10.  The  right-running  and  left-running  Mach  lines  are 
not  identical,  however,  since  the  special  point  at  x  =  8.07104661  in. 
considered  in  Sample  Case  No.  10  is  not  considered  in  Sample  Case  No.  11. 
The  flowfield  downstream  of  left* running  Mach  line  29  in  Sample  Case  No. 
11  is  presented  in  Append.  6. 
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13.  SAMPLE  CASE  NO.  12 

Sample  Case  No.  12  is  identical  to  Sample  Case  No.  11  except  that 
the  left-running  characteristic  network  is  initiated  from  the  last  right- 
running  characteristic  from  the  initial-value  line  instead  of  from  the 
last  right-running  characteristic  from  the  initial  expansion  contour. 
Since  Sample  Case  No.  11  involved  a  basic  nozzle  only  (MODE  =  1),  this 
sample  case  also  involves  only  a  basic  nozzle.  The  analysis  of  a  basic 
only  by  emanating  left-running  characteristics  from  the  downstream 
extent  of  the  Initial -value  line  is  selected  by  specifying  H0DE®2.  Such 
an  analysis  in  a  basic  nozzle  is  analogous  to  H0DE=5  in  a  scarfed  nozzle. 
The  only  change  required  to  the  data  decx  for  Sample  Case  No.  11  to 
obtain  Sample  Case  No.  12  is  to  change  MODE  from  1  to  2.  The  data  deck 
for  Sample  Case  No,  12  is  presented  below. 

SAMPLE  CASE  NO.  12.  SAMPLE  CASE  NO.  11  WITH  M0DE«2 
SDATA  M0DE=2,  IWALL=1 ,  NI=11,  NT=15,  KWRITE=1, 

6=1.2,  RG=65.0,  PS=1000.0,  TS=50G0.0,  PA=14.696, 

YT°1 .0,  RTU=1 .0,  RTD=0.01,  AA-15.0,  XE=28. 51 233074  $ 

The  output  format  for  Sample  Case  No.  12  is  similar  to  the  output 
format  for  Sample  Case  No.  6  presented  in  Section  VIII. 7. 


14.  SAMPLE  CASE  NO.  13 

Sample  Case  No.  13  illustrates  the  quadratic  wall  option  for  the 
basic  nozzle  (IWALL=2).  The  envelope  of  the  quadratic  wall  is  identical 
to  the  envelope  of  the  conical  nozzle  considered  in  Sample  Case  No.  1. 

The  throat  attachment  angle  is  ea  =  25.0  deg,  25  points  are  to  be  specie 
fied  along  the  throat  downstream  circular  arc  contour,  the  nozzle  area 
ratio  e  =  10.0,  and  the  nozzle  length  is  xe  =  8.07104661  in.  (AA=25.0, 
NT=25.0,  EPS=10.0,  and  X£=8. 07 104661).  The  same  cylindrical  nozzle 
extension  considered  in  Sample  Case  No.  1  is  considered  here.  All  of  the 
remaining  parameters  in  namelist  DATA  are  the  same  as  for  Sample  Case  No. 
1.  The  data  deck  for  Sample  Case  No.  12  is  presented  below. 

SAMPLE  CASE  NO.  13.  QUADRATIC  WALL,  AA=25.0,  XE=8. 07104661 ,  EPS=10.0 
SDATA  H00E=4,  IWALL-2,  NI=U,  NT=25,  KWRITE=1, 

6=1.2,  RG=65.0,  PS=1000.0,  TS=500G.0,  PA=14.696, 

YT=1 .0,  RTtKI.0,  RTD=0.01,  AA=25.0,  XE=8. 071 04661 ,  EPS* 'l  0.0, 

AF=0.0,  6ETA=30.0  % 

The  output  for  Sample  Case  No.  13  is  presented  in  Append.  H. 
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15.  SAMPLE  CASE  NO,  14 

Sample  Case  No.  14  is  identical  to  Sample  Case  No.  13,  except  that 
the  quadratic  portion  of  the  basic  nozzle  is  specified  by  the  exit  lip 
radius  ye  instead  of  the  area  ratio  e.  The  identical  quadratic  contour 
is  obtained  for  ye  =  3.16227766  in.  Thus,  the  only  changes  to  the  input 
data  are  to  delate  EPS  =  10.0  and  to  add  YE  =  3.16227766.  The  data  deck 
for  Sample  Case  No.  14  is  presented  below. 

SAMPLE  CASE  NO.  14.  QUADRATIC  WALL,  AA=25.0,  XE=8. 071 04661 ,  YE=3. 16227766 
$DATAM0DE=4,  IWALL=2,  NI=11,  NT=25,  KWRITE=1, 

6=1.2,  RG=65.0,  PS=1000.0,  TS=5000.0,  PA=14.696, 

YT=1 .0,  RTU=1 .0,  RTD=0.01 ,  AA=25.0,  XE-d. 071 04661 ,  YE*3. 16227766, 

AF=0.0,  BETA=30.0  $ 


The  output  for  Sample  Case  No.  14  is  identical  to  the  output  for 


16.  SAMPLE  CASE  NO.  15. 

Sample  Case  No.  15  is  identical  to  Sample  Case  No.  13,  except  that 
the  quadratic  portion  of  the  basic  nozzle  is  specified  by  the  exit  lip 
angle  ee  instead  of  the  area  ratio  e.  The  identical  quadratic  contour 
is  obtained  for  6e  -  3.97861353  deg.  Thus,  the  only  changes  to  the  input 
data  are  to  delete  EPS  =  10.0  and  to  add  AE  =  3.97861353.  The  data  deck 
for  Sample  Case  No.  15  is  presented  below. 

SAMPLE  CASE  NO,  15.  QUADRATIC  WALL,  AA=25.0,  XE=8. 071 04662,  AE=3. 97861 353 
$DATA  M0DE=4,  IWALL=2,  NI=11,  NT=25,  KWRITE=1, 

6=1.2,  RG=65.0,  PS=1 000.0,  TS=5000.0,  PA=14.696, 

YT=1 .0,  RTU=1 .0,  RTD=0 . 01 ,  AA=25.0,  XE=8. 071 04661 ,  AE=3. 97861 353, 

AF=0.0,  BETA=30.0  $ 

The  output  for  Sample  Case  No.  15  is  identical  to  the  output  for 
Sample  Case  No.  13. 


SECTION  IX 
CONCLUSIONS 


An  analysis  is  presented  for  predicting  the  performance  of 


canted  scarfed  propulsive  nozzles.  A  computer  program  based  on  that 
analysis  was  developed.  Fifteen  sample  cases  are  presented  to 
illustrate  the  use  of  the  computer  program. 
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